Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shallow”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

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

Chang, Hilary↗
DOE OSTI2021

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗
DOE OSTI2022

Report on Depth-Dependent Q from Frequency-Dependent Lg Q

As a phase that travels through the crustal waveguide, Lg samples the velocity and attenuation structure of the crust. We have developed methods (Pasyanos et al., 2009a; Pasyanos et al., 2009b) to estimate Qp and Qs of the crust and upper mantle from the amplitudes of regional phases, and have applied it to a number of regions around the world, and over a broad frequency band (0.5 - 10 Hz). In prior attenuation modeling, Lg phase amplitudes are able to effectively distinguish the high Q (low attenuation) of old stable cratonic crust from the low Q (high attenuation) of more recently active tectonic regions. For a variety of purposes (e.g. waveform modeling, seismic hazard, etc.), it is useful to have more information on the depth-dependence of the Q structure in the crust. One obvious candidate for doing this is using the frequency-dependence information of Lg Q, in the manner of surface waves. Unfortunately, however, Lg samples the whole of the crust without much differentiation with respect to frequency that is characteristic of fundamental mode surface waves. Lg is generally thought to be the sum of supercritically reflected S-waves trapped in the crustal waveguide (e.g. Xie and Lay, 1994). The phase can also be modeled as a superposition of surface waves. At lower frequencies, it is primarily composed of fundamental and low-order surface waves, while at higher frequencies, it is comprised purely of many high-order surface waves. The low-order surface waves that comprise Lg in low frequency bands (< 1 Hz) are preferentially sensitive to shallow earth structure, such as sedimentary basins and upper crust, while high-order surface waves have a more homogenized sensitivity to the whole crust. Furthermore, at the lower frequencies, we might expect to see some differentiation in the sensitivity between early Lg (herein referred to as Lg1) and late Lg (referred to as Lg2) where the latter is comprised of lower-order surface waves and, hence, shallower structure. Therefore, we might expect to see the signature of shallow features, primarily the lower Q of sediments as compared to crystalline crust, at lower frequencies and in later arriving Lg amplitudes.

58 GEOSCIENCES↗
DOE OSTI2022

200-IA-1 Operable Unit Ecological Risk Screen

The purpose of this environmental calculation file (ECF) is to document the assumptions, inputs, methodology, equations, and results of a screening level ecological risk assessment (SLERA) conducted for eight of 12 representative waste sites within the newly formed 200 IA-1 Operable Unit (OU); the remaining four waste sites lacked soil data within the shallow vadose zone and were not further evaluated. The SLERA was conducted to evaluate ecological receptor exposure and potential risk using soil data for radionuclides and nonradionuclides occurring within the shallow vadose zone (0 to 15 ft) soil of the OU. The representative Hanford ecological receptors evaluated are: plants, soil invertebrates, California quail (Callipepla californica), Meadowlark (Sturnella neglecta), Killdeer (Charadrius vociferus), Red-tailed hawk (Buteo jamaicensis), Great Basin pocket mouse (Perognathus parvus), Deer mouse (Peromyscus maniculatus), Grasshopper mouse (Onychomys leucogaster), and badger (Taxidea taxus). During the SLERA, detected soil sample results from the shallow vadose zone were used to calculate Hazard Quotients (HQs) for nonradionuclides and Sum of Fractions (SOF) for radionuclides. The results of this SLERA support the preparation of DOE/RL-2020-51, 200-IA-1 Operable Unit Focused Feasibility Study.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗
DOE OSTI2023

Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake

An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.

58 GEOSCIENCES↗
DOE OSTI2023

Automated Quantification of Wind Turbine Blade Leading Edge Erosion from Field Images

Wind turbine blade leading edge erosion is a major source of power production loss and early detection benefits optimization of repair strategies. Two machine learning (ML) models are developed and evaluated for automated quantification of the areal extent, morphology and nature (deep, shallow) of damage from field images. The supervised ML model employs convolutional neural networks (CNN) and learns features (specific types of damage) present in an annotated set of training images. The unsupervised approach aggregates pixel intensity thresholding with calculation of pixel-by-pixel shadow ratio (PTS) to independently identify features within images. The models are developed and tested using a dataset of 140 field images. The images sample across a range of blade orientation, aspect ratio, lighting and resolution. Each model (CNN v PTS) is applied to quantify the percent area of the visible blade that is damaged and classifies the damage into deep or shallow using only the images as input. Both models successfully identify approximately 65% of total damage area in the independent images, and both perform better at quantifying deep damage. The CNN is more successful at identifying shallow damage and exhibits better performance when applied to the images after they are preprocessed to a common blade orientation.

Aird, Jeanie A.↗
DOE OSTI2022

Pitting Propagation Behavior on Low Alloy AISI 4130 (UNS G41300) Steel Exposed to Various Alkali and Alkaline Earth Metal Chlorides

Pit propagation studies were conducted to elucidate whether alkali and alkali earth metal cations such as Na+ and Rb+ present in the form of metal chloride salts such as RbCl affect pitting behavior in distinctly different manner than NaCl. Pit propagation studies were conducted on a low alloy steel using one-dimensional (1-D) pit method over pit depths from 300-1000 µm. LSV and EIS of planar electrodes of 4130 in a range of Cl- solutions were conducted and revealed no detected differences in impedance, open circuit, corrosion potential (Ecorr), passive current density (ipass), and pitting potential (Epit) as a function of salt type. In the case of one-dimensional pits during fast downward scan rates, the saturation potential (Esat) varied as a function of cation identity when pit depths were shallow. Mass transported limited current density also differed with cations in shallow pits when various alkali metal and alkaline metal cations were present. The limiting current density was increased for RbCl relative to NaCl The pit surface potential (Esurf) of activated pit surfaces reached Ecorr prior to establishing a condition where the pit electrolyte surface concentration (Csurf) was less than the critical concentration for active acidified pitting (i.e. Csurf<Ccrit) in this marginally passivating steel. For various Esurf and pit current density (ipit) combinations at constant Csurf where Ccrit< Csurf < Csat, E-log(i) plots were constructed using the method of Tianshu to create IR ohmic voltage corrected Tafel plots for various fixed pit solution concentrations. Under these conditions, the influence of cation identity on charge transfer controlled kinetics indicated slight differences in Tafel behavior where RbCl was slightly more aggressive than NaCl for low alloy steel oxidation. Differences in metal cation identity exert no effect on passive or breakdown on planar electrodes and only effect pit propagation stage in shallow pits.

Demarest, Chalres↗
DOE OSTI2024

Environmental controls on isolated convection during the Amazonian wet season

The Amazon rainforest is a vital component of the global climate system, influencing the hydrological cycle and tropical circulation. However, understanding and modeling the evolution of convection in this region remain a scientific challenge. Here, we assess the environmental conditions associated with shallow, congestus, and isolated deep convection days during the wet season (December to April), employing measurements from the Green Ocean Amazon 2014–2015 (GoAmazon2014/5) experiment and large-scale wind fields from the constrained variational analysis. Composites of deep days show moister than average conditions below 3 km early in the morning. Analyzing the water budget at the surface through observations only, we estimated the water vapor convergence term as a residual of the water balance closure. Convergence remains nearly zero during the deep days until early afternoon (13:00 LST), when it becomes a dominant factor in the water budget. At 14:00 LST, the deep days experience a robust upward large-scale vertical velocity, especially above 4 km, which supports the shallow-to-deep convective transition occurring around 16:00–17:00 LST. In contrast, shallow and congestus days exhibit drier pre-convective conditions, along with diurnal water vapor divergence and large-scale subsidence that extend from the surface to the lower free troposphere. Moreover, afternoon precipitation exhibits the strongest linear correlation (0.6) with large-scale vertical velocity, nearly double the magnitude observed for other environmental factors, even moisture, at different levels and periods of the day. Precipitation also exhibits a moderate increase with low-level wind shear, while upper-level shear has a relatively minor negative impact on convection.

Environmental Sciences & Ecology↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Viscous stratification of the earth and convection

The shallow model of the earth's mantle is discussed along with a variety of geophysical arguments for its correctness and against the existence of deep convection. The main agrument is summarized in the proposal that the astheno sphere is less viscous (by a factor of 10 to 100) than has generally been assumed. In this shallow model, the return flow is essentially through the asthenosphere. The dynamical agent is the steep temperature gradient in the upper mantle. Speculations as to the historical variation of this gradient are advanced. The effects on the model of a nonuniform earth aggregation are considered and shown to favor shallow convection as well as a top convective layer (lithosphere plus asthenosphere) whose depth increases slowly over the earth's life, leading to a tectonic activity that increases gradually with time.

Elsasser, W. M.↗
NASA NTRS1974

Shell deformation studies using holographic interferometry

The buckling of shallow spherical shells under pressure has been the subject of many theoretical and experimental papers. Experimental data above the theoretical buckling load of Huang have given rise to speculation that shallow shell theory may not adequately predict the stability of nonsymmetric modes in higher-rise shells which are normally classified as shallow by the Reissner criterion. This article considers holographic interferometry as a noncontact, high-resolution method of measuring prebuckling deformations. Prebuckling deformations of a lambda = 9, h/b = 0.038 shell are Fourier-analyzed. Buckling is found to occur in an N = 5 mode as predicted by Huang's theory. The N = 4 mode was unusually stable, suggesting that even at this low value of h/b, stabilizing effects may be at work.

Parmerter, R. R.↗
NASA NTRS1976

Formation of ferrobasalt at east Pacific midocean spreading centers

A numerical modeling using a least squares petrologic mixing program demonstrates that ridge ferrobasalts from three east Pacific regions are generated from normal midocean ridge balsalt by crystal fractionation of plagioclase, clinopyroxene, and minor olivine at shallow mantle depths. Clinopyroxene is an important fractionating phase in this model in contrast to most previously proposed fractionation models. Calculations demonstrate that the ferrobasalts represent a much smaller percentage of residual liquid than has previously been proposed. Implications of extensive shallow crystal fractionation for the determination of depths of magma generation are discussed. The different fields on an anorthite-olivine-hypersthene ternary diagram for basalts from east Pacific fast and moderate spreading ridges can be explained by variable degrees of shallow crystal fractionation rather than by different depths of magma generation as suggested by Bass (1971) and Scheidegger (1973).

Clague, D. A.↗
NASA NTRS1976

Martian crater depth/diameter relationships - Comparison with the moon and Mercury

New Martian topographic data from Mariner 9 ultraviolet spectrometer (UVS) profiles provide depth data for 139 Martian craters of all degrees of degradation, between the diameters of 15 and 201 km. The population of Martian craters, including morphologically fresh examples, is shallower than both lunar and Mercurian fresh crater populations. Because the surface gravities of Mercury and Mars are identical within 5%, these differences in fresh crater depths suggest that factors other than gravity may play important roles in determining initial crater depths (e.g., differences in impact velocity, substrate variations, and Martian atmospheric effects during the crater-forming event). Degraded Martian craters are, on the average, no shallower than lunar pre-Imbrian craters of similar sizes. If the early bombardment of Mars was as significant a degradational agent as it was on the moon, then major levels of crater degradation and crater shallowing on Mars were associated with this mechanism. Continued eolian infilling, although locally significant, may be a less significant cause of morphometric degradation of large old Martian craters.

Cintala, M. J.↗
NASA NTRS1977

New processes and materials for ultraviolet detection with solid state devices

The three major effects that degrade external responsivity of silicon from the 1/lambda theoretical curve for a quantum detector are: surface reflectance, surface recombination, and junction depth. Since the p-n junction must be very shallow, problems relating to surface are further enhanced. MOS type of processing is necessary. HCl oxides and numerous acid clean-ups are utilized in order to obtain a contamination free surface with low Qss levels. Stringent process controls such as CV shifts, spreading resistance measurements, thickness monitoring etc., are used to analyze the surface contaminations, surface mobile charges, surface concentrations, junction depth, oxide thickness etc. Low surface concentrations of 10 to the 18th atoms/cu cm are achieved by low temperature boron nitride depositions. Shallow junction depths of the order of a few tenths of a micron are achieved by low temperature controlled diffusions. In order to improve breakdown characteristics of these shallow junction devices, field plate and deep diffused p(+) ring geometries are used.

Chopra, D.↗
NASA NTRS1981

Seismic energy release of the moon

Lunar seismicity is investigated by calculating various source parameters for a number of shallow and deep-focus moonquakes. The seismic moment, seismic energy release, annual seismic energy release, stress drop, and body-wave magnitude are determined for the largest shallow moonquakes and for large deep-focus events. It is found that the shallow events dominate the lunar seismic energy release, that tidal dissipation may account for the energy release by the deep-focus events, and that the stress drops for the deep-focus events are comparable to or smaller than the calculated tidal stresses. A comparison of the results with terrestrial data indicates that the seismic characteristics of a planet are controlled more by tectonic style and state than by the relative magnitude of the driving forces.

Goins, N. R.↗