Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface modeling”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Bayesian Learning of Adatom Interactions from Atomically Resolved Imaging Data

Atomic structures and adatom geometries of surfaces encode information about the thermodynamics and kinetics of the processes that lead to their formation, and which can be captured by a generative physical model. In this work, we develop a workflow based on a machine-learning-based analysis of scanning tunneling microscopy images to reconstruct the atomic and adatom positions, and a Bayesian optimization procedure to minimize statistical distance between the chosen physical models and experimental observations. We optimize the parameters of a 2- and 3-parameter Ising model describing surface ordering and use the derived generative model to make predictions across the parameter space. For concentration dependence, we compare the predicted morphologies at different adatom concentrations with the dissimilar regions on the sample surfaces that serendipitously had different adatom concentrations. The proposed workflow can be used to reconstruct the thermodynamic models and associated uncertainties from the experimental observations of materials microstructures. The code used in the manuscript is available at https://github.com/saimani5/Adatom_interactions.

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

Insights into the Interaction of Redox Active Organic Molecules and Solvents with the Pristine and Defective Graphene Surfaces from Density Functional Theory

A systematic study is reported on the interaction of two representative redox active organic molecules and two solvent molecules with pristine and defective graphene surfaces as a model of an electrode surface of a redox flow battery (RFB). The redox active molecules include a catholyte, 2,5-di-tert-butyl-1,4-dimethoxybenzene (DDB), and an anolyte molecule, benzothiadiazole (BTZ), and the solvent molecules include acetonitrile (MeCN) and ethylene carbonate (EC). The graphene defects investigated include a single vacancy, double vacancy, zigzag step edge, and armchair step edge. Computations suggest that the interactions of all molecules with a pristine graphene surface are relatively weak (0.2 to 0.8 eV) and dominated by van der Waals effects; therefore, these molecules are chemically stable upon interacting with pristine nondefective graphene. The BTZ, DDB, MeCN, and EC molecules interact strongly (1.5 to 5.5 eV) with the single vacancy and zigzag step edges of graphene that leads to the possible decomposition of the molecules with strength of interaction in the order of MeCN > BTZ > DDB > EC and MeCN > EC > BTZ > DDB, respectively. Calculations show that the BTZ, DDB, MeCN, and EC molecules interact less strongly (0.2 to 1.4 eV) to the double vacancy and armchair step edge than to the single vacancy and zigzag step edge. The binding energies of the molecules were significantly reduced when interacting with the passivated defects, suggesting that the passivation of the defects could help prevent unwanted chemical interactions between the neutral molecules and the electrode surface. In all organic RFBs, one of the crucial bottlenecks is the stability of the constituent molecules with the electrode material, and this study provides insights into the chemical interaction of selected candidate species with a model carbon electrode.

Howard, Jason D.↗
DOE OSTI2021

Simulating Oxidation of Carbon Surfaces by Atomic Oxygen Coupled with a Finite Rate Oxidation Model

Oxidation is a prominent factor in degradation of carbon based ablative materials at high temperatures, resulting in mass loss and material recession through chemical and mechanical erosion. In this paper, a finite rate carbon oxidation model coupled with surface recession is implemented in the Porous Ablation and Thermal Response (Panther) software to investigate carbon oxidation by atomic oxygen over a range of temperatures and pressures. Panther is utilized to investigate the finite rate oxidation model proposed by Swaminathan-Gopalan et al. at conditions of higher pressure and flux relative to the experiments it was built on.

Kirmanoglu, Kaan K.↗
DOE OSTI2021

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Size effect in stainless steel thin wires under tension

Miniaturization of modern devices drives the exploration into how materials behave in small dimensions. However, the mechanical behaviors of micro-scale specimens under varied loading conditions are still challenging to characterize and predict. In this work, specimen size-dependent weakening effect was experimentally revealed in 316L stainless steel thin wires under tension with wire diameter decreasing from 200 μm down to 61 μm. The underlying mechanisms were primarily associated with surface grain softening and strain localization, which were largely dependent on the number of grains across wire diameter. Through an analog composite model with surface and interior grains that follows the Hall-Petch relation, the extent of surface grain softening was found to coincide with the flow stress reduction due to lack of boundary strengthening. Further using elasto-plastic self-consistent polycrystal modeling, the stress-strain behaviors of interior and surface grains were well described in analogy to those of bulk materials with similar-size grains and extremely large grains, respectively. Subsequently, the analog composite model well simulated the specimen size-dependent weakening effect under tension. This article provided more insights into the size effect exploration and proposed a practical model for describing size-dependent mechanical behaviors of small-scale materials to be applied to miniaturized devices.

316L stainless steel↗
DOE OSTI2019

Inertial Effects During the Process of Supercritical CO 2 Displacing Brine in a Sandstone: Lattice Boltzmann Simulations Based on the Continuum-Surface-Force and Geometrical Wetting Models

Inertial effects during the process of supercritical CO 2 displacing brine in porous media may not be negligible according to recent studies. Capturing the inertial effects of the physical CO 2 -brine system imposes a requirement on the grid resolution and viscosity to surface tension ratio for pore-scale simulations, which some commonly used simulators may not be able to meet. To fulfill the parameter requirement, we combine the continuum-surface-force based color-gradient lattice Boltzmann (LB) multiphase model and the geometrical wetting model and extend the model to 3D under the multiple-relaxation-time framework. We validate the model via simple benchmarks which show significant improvement over the traditional models. Furthermore, we perform 3D drainage simulations in a heterogeneous micromodel where the simulation result agrees well with experimental data, while our previous work fails to reproduce certain displacement patterns in the experiment due to the use of a traditional LB model that cannot fulfill the parameter requirement. Finally, we perform high-fidelity 3D drainage simulations to study the inertial effects in a Bentheimer sandstone sample. These findings show that stronger inertial effects generally help develop more CO 2 flow pathways for the same capillary number which results in higher CO 2 saturation, consistent with the micromodel results. The phenomena can be found in both low and high capillary number cases, indicating that the inertial effects are not dependent on the mean velocity. In addition, the change of the invasion patterns is not proportional to the change of inertial effects, thus exhibiting threshold behavior.

58 GEOSCIENCES↗
DOE OSTI2025

Continuous snow depth and temperature measurements from dense network of above-ground distributed temperature profiling systems from 2021-09-23 to 2024-08-23, Seward Peninsula, Alaska

The dataset contains temperature measurements from distributed temperature profiling (DTP) systems (Dafflon et al., 2022; Wielandt et al., 2022; Wang et al., 2024a; Fiolleau et al., 2024) deployed vertically above the ground surface at a large number of locations from 2021 to 2024. The research is designed to improve understanding of the local heterogeneity in snow depth and snow thermal insulation dynamics, as well as their interactions in a discontinuous permafrost region (Wang et al., 2025). The DTP systems were deployed at 96 locations in a watershed along the Nome-Teller road at mile marker 27 (T27) and at 54 locations on a hillslope along the Kougarok road at mile marker 64 (K64) in the Seward Peninsula, Alaska. The probe location information is stored in Probe_locations_T27.csv and Probe_locations_K64.csv. Temperature measurements were recorded at 15-minute intervals using high-precision digital sensors (accuracy: ±0.1°C, resolution: 0.0078°C). The temperature probes, either 1.4 m or 1.6 m long, contain sensors spaced every 5 cm or 10 cm along their length. The temperature data are stored in compressed files following the format: DTP_snow_air_temperature_(site)_(start)_(end).zip, where site is either T27 or K64, and start and end represent the time series period. Within each ZIP file, individual CSV files are named by probe ID and contain temperature records at different heights above the ground surface.This dataset also includes derived snow depth time series over three snow seasons, estimated from temperature measurements. Snow depth was estimated by identifying the consecutive sensor pair that exhibited the largest drop in high-frequency temperature fluctuations (detailed in the methods). These data are stored in: Snow_depths_flags_(site)_(start)_(end).csv, which includes snow depth time series and corresponding quality flags (defined in the methods) from different probes. Additionally, the dataset includes derived metrics and supporting measurements at selected locations over two snow seasons, contributing to the manuscript of Wang et al., 2025. These locations were chosen based on the availability of high-quality snow depth time series during both seasons. The additional data include: (1) Air temperature proxies measured from the top sensors on the pole when they were not buried by snow, stored in Air_temperature_proxies_(site)_(start)_(end).csv (2) Ground interface temperature, recorded at 3 cm above the ground, stored in Ground_interface_temperature_(site)_(start)_(end).csv (3) Site characteristics, including vegetation height, elevation, and the topographic position index (TPI) within a 50 m radius, stored in Selected_probe_locations_gps_vegheight_tpi_elevation_(site).csv. These metrics were derived from 1 m resolution summer LiDAR-based digital elevation models and digital surface models from Singhania et al., 2023, DOI:10.5440/1832016. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv.This dataset is an updated version of a previous archive (Wang et al., 2024b, DOI: 10.15485/2475020), incorporating multiple seasons and improved snow depth estimation. Please note that due to large amount of information present in this dataset, many specificities associated with the acquisition of snow temperature, air temperature proxy and estimation of snow depth, and the future archiving of additional datasets on the soil temperature, thaw depth and soil characteristics at these locations, the author would welcome being contacted by people planning to use this dataset.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗
DOE OSTI2024

Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States

Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.

Geochemistry & Geophysics↗
DOE OSTI2021

Effect of electrode rinse solutions on the electrodialysis of concentrated salts

The data presented indicate that modest attention to the chemistry of the electrode rinse solution can yield as much as 56% improvement in process efficiency. When operating electrodialysis (ED) with salt concentrations greater than 10,000 mg/l (as NaCl) the amperage can be relatively high even at low voltages. To maximize the overall rate of ion transport from diluate to concentrate, there is a need to minimize resistances in the electrode cells since these can represent as much as 30% of the total resistances in the entire process and in this particular ED unit, represented resistance roughly equivalent to that within the membrane stack. This work describes the performance of a 10-cell pair, 200 cm 2 (0.02 m 2 ) per membrane, pilot ED unit operated in batch mode at 5 V potential. The feedstock to the stack was varied from 0.5% to 6% NaCl. Results from Volt-amp profiles (2.0–15 V) were used as the rationale for choosing the preferred electrode rinse solution, 90 g/l (kg/m 3 ) disodium sulfate at pH 12.5. The recommended solution of 30 g/l (kg/m 3 ) disodium sulfate at neutral pH was tested against stronger solutions (60, 90, and 120 g/l (kg/m 3 ) disodium sulfate) at neutral pH and with 1 g/l (kg/m 3 ) sodium hydroxide added to yield solutions at pH 12.5. There were stark differences in the performance of the ED unit between the solutions at pH 7 and those at pH 12.5. The slopes of the Volt-amp profiles improved indicating less resistance to ion flow at pH 12.5. Most of the difference can be attributed to the greater conductivity of sodium hydroxide compared to disodium sulfate, where even a modest addition of 1 g/l (kg/m 3 ) sodium hydroxide increased the conductivity of the solution between 8 and 20% depending on the concentration of disodium sulfate. Most notably, the initiation voltage was 2.39 ± 0.08 V with solutions at pH 7 and 1.99 V ± 0.04 V when the rinse solutions were around pH 12.5, a difference of around 0.42 V. Since the full batch runs were performed at 5.0 V, this 0.42 V differential represented a loss of 8% in process efficiency. A series of tests showed a shift in initiation voltage occurred between pH 11.5 and 12.5. Nernst equations for water electrolysis are coupled with a simple flux model that estimates the concentration of hydroxide and hydronium ions at the electrode surfaces. This model indicates that the potential demanded by the electrodes is directly linked to the current. Furthermore, there is a steep drop in the potential demanded by the anode when there is sufficient hydroxide available to generate oxygen from hydroxide as opposed to generation from water. The voltage demand is minimal at pH 12.5 and reflective of the reference standard voltage (just over 1.23 V predicted at 0.1 amps). However, the flux model predicts that the voltage demanded by the electrodes at pH 7 is 1.61 V at 0.1 amps (the lowest amperage measurable with the power pack). In conclusion, comparing this result to the reference standard (E 0 = -1.23 V) yields a satisfactory explanation for the observed difference below reference standard of 0.42 V.

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

RiverPIXELS: paired Landsat images and expert-labeled sediment and water pixels for a selection of rivers v1.0

RiverPIXELS contains GeoTIFFs of hand-labeled water and sediment pixels from Landsat images containing rivers. Each of the 104 labeled patches contains 256 x 256 Landsat pixels (30 meter resolution). Our aim in releasing RiverPIXELS is to provide an "off-the-shelf" training and testing dataset for building machine-learned models to automatically identify rivers from multispectral imagery. While a number of trained models and/or surface water products already exist, RiverPIXELS aims for pixel-level accuracy in order to precisely identify river boundaries in particular. Our selection of rivers is heavily Arctic, but we include tropical and temperate rivers as well. Patches are provided for the Colville (7), Indigirka (6), Kolyma (4), Ucayali (54), Waitaki (21), and Yana (12) Rivers. For each patch, all surface water pixels are labeled (1) and all in-channel sediment pixels are labeled (2). Sediments not in-channel are considered part of the land (0) class. RiverPIXELS also includes paired surface water data from the Global Surface Water dataset that may be useful as additional features in machine learning models. Each patch therefore contains four aligned GeoTIFFs: labeled, landsat, gswmo, and gswocc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Generalized Hall conductivities in local commuting projector models: Generalized symmetries and protected surface modes

Hall conductivities are important characterizations of phases of matter. It is known that nonzero Hall conductivities are difficult to realize in local commuting projector lattice models due to no-go theorems in (2+1)⁢D. In this work we construct local commuting projector models in (2+1)⁢D and (3+1)⁢D with nonzero generalized Hall conductivities for ordinary and higher-form continuous symmetries on tensor product Hilbert space of finite local dimension. The model is given by a standard ℤ 𝑁 toric code, but the symmetries do not admit expression in terms of on-site charge operators. The symmetry do not have local charges or currents on the lattice in the absence of boundaries, but there is still a notion of Hall conductivities that coincide with the continuum field theories. We construct protected gapless boundaries of the lattice models using modified Villain formalism. The generalized Hall conductivities are computed by surface currents as well as bulk flux insertion and many-body Chern number.

Anomalies↗
DOE OSTI2023

Amazon drought resilience – emerging results point to new empirical needs

Will climate change turn tropical forests from a carbon sink to a source (Pan et al., 2011)? Tropical forests cycle more carbon, water, and energy than any other biome (Bonan, 2008). The Amazon is the largest tropical forest, spanning one-third of South America, and any perturbation in its function has repercussions on the global climate. Mild and severe droughts in the Amazon are predicted to double and triple in area, respectively, by 2100, while the area under wet extremes will increase (Duffy et al., 2015). Our understanding of the role of plants' hydrological environments in determining their drought response, however, is limited (Chitra-Tarak et al., 2018, 2021). An important review by Costa et al. (2023; pp. 714–733) published in this issue of New Phytologist highlights that shallow-water table (WT) forests constitute c. 50% of the Amazon and may act as hydrological refugia during droughts, yet they have been neglected in Amazon forest research to date. Emerging results suggest that the shallow-WT forests that dominate the Amazon basin may increase in productivity under moderate droughts due to relief from hypoxia, suggesting a potential buffer from drought relative to deeper WT forests (Fig. 1). During severe droughts, however, shallow-WT forests may be vulnerable to collapse due to drought-intolerant traits. Addressing the underrepresentation of relatively drought-resilient shallow-WT forests in the Amazon's forest inventories may help resolve the much-debated incongruence in recent studies of the Amazon's drought resilience. Data from forest inventories (biased toward deep-WT forests) indicate that the Amazon's carbon sink is declining, partially because of drought impacts (Phillips et al., 2009; Brienen et al., 2015), whereas basin-wide satellite-based measures of gross primary productivity, which include shallow-WT forests, indicate varied regional responses to droughts (Saleska et al., 2007; Brando et al., 2010). Costa et al.'s (2023) descriptive insights of the intricate balance of geology, topography, hydrology, vegetation, and drought on ecosystem function, as well as their conceptual predictive framework, are useful for developing new empirical research in these understudied ecosystems and improving Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Review of Wind–Wave Coupling Models for Large-Eddy Simulation of the Marine Atmospheric Boundary Layer

Abstract We present a review of existing wind–wave coupling models and parameterizations used for large-eddy simulation of the marine atmospheric boundary layer. The models are classified into two main categories: (i) the wave-phase-averaged, sea surface–roughness models and (ii) the wave-phase-resolved models. Both categories are discussed from their implementation, validity, and computational efficiency viewpoints, with emphasis given on their applicability in offshore wind energy problems. In addition to the various models discussed, a review of laboratory-scale and field-measurement databases is presented thereafter. The majority of the presented data have been gathered over many decades of studying air–sea interaction phenomena, with the most recent ones compiled to reflect an offshore wind energy perspective. Both provide valuable data for model validation. We also discuss the modeling knowledge gaps and computational challenges ahead.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impact of the Earthquake Rupture on Ground-Motion Variability of the 24 August 2016 M w 6.2 Amatrice, Italy, Earthquake

Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.

3D Ground Motions Simulations↗
DOE OSTI2019

Surface Engineering by Simultaneous Action of Multiple External Fields (Final Technical Report)

This research project aims at enabling surface engineering strategies based on the fundamental understanding and predictive modeling of the surface morphological response of solid materials subjected to the combined action of multiple external forces. Surface morphological response to mismatch strain in strained-layer heterostructures also is analyzed. Our research focuses on the surface morphological evolution and stability of electrically conducting (metallic, such as Cu) and semiconducting (such as Si) solids under the simultaneous application of mechanical stresses, electric fields, and temperature gradients. We emphasize on identifying the conditions under which the multiply driven surface morphology is stable and on exploring the complexity of the corresponding various morphologically stable surface patterns. We also emphasize on the driven evolution of nanoscale features aiming at stabilizing and controlling nanoscale patterns on surfaces by their manipulation through simultaneously applied multiple external forces. Specifically, we address systematically the morphological response to the combined action of multiple external forces of bulk solid surfaces, surfaces of thin films grown epitaxially on thick or thin substrates, as well as nanoscale surface features such as coherently strained islands grown epitaxially on substrate surfaces. The research is based on a modeling approach that combines theoretical analyses of surface morphological stability with self-consistent dynamical simulations of surface evolution based on properly parameterized continuum and multiscale surface transport models.

36 MATERIALS SCIENCE↗