Search NASA⌕ Search

SEARCH · Search NASA

Results for “DISCONTINUITY SURFACES”

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

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2022

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Investigation of Debonding Effect in Internal Replacement Pipe System Under Lateral Loading

The aging pipeline infrastructure around the world necessitates immediate rehabilitation. Internal replacement pipe (IRP) is a trenchless system offering a versatile and cost-effective solution across a variety of industries, including oil, natural gas, water, and wastewater. As a structural pipeline repair system, IRPs are subject to lateral deformation because of surface traffic loading. The present study evaluates the impact of adhesion between the host pipe and the IRP, with a focus on assessing the debonding effect on the behavior of the repair system under lateral deformation and bending. This was achieved using a comprehensive approach, including experimental, numerical, and analytical techniques. Varying levels of adhesive strength resulting from different methods of surface preparation were considered. The effectiveness of the IRP system on both discontinuous host pipes with various crack widths and continuous host pipes was also investigated. The results demonstrate that adhesive strength exerts a significant influence on the repair system, especially in the case of narrow circumferential cracks, while its impact on the continuous system is minimal. For optimal performance, it is essential to choose adhesives that possess sufficient shear strength while also accounting for the required debonding length. This approach ensures that minor discontinuities are effectively controlled, thereby enhancing the system′s fatigue life. The reliable determination of the maximum allowable shear strength for the adhesive or the debonding length can ensure that it does not negatively affect fatigue life. The findings presented in this study offer new insights into the development of trenchless repair techniques that can enhance system performance and extend service life.

Tien, Tri C. M.↗
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

Near-Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub-Arctic Watershed: Modeling Archive

This Modeling Archive is in support of a NGEE-Arctic publication: Shirley et al. (2022) “Near-Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub-Arctic Watershed". [DOI].The dataset contains outputs from the global sensitivity analysis (GSA) of the “ecosys” model as reported in Shirley et al. (2022). The study showed that discontinuous permafrost environments are characterized by complex feedback loops and strong spatial heterogeneity which is created by variability in near-surface hydrology and soil properties. Additionally, the study demonstrated that missing representation of sub-grid heterogeneity in terrestrial ecosystem models can lead to biased estimates of the high-latitude carbon budget. Included in this dataset are the factor values for each run in the GSA and the model outputs used in this study. Included are two *.csv data files and one *.pdf.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 OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2021

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

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

Chiral Topological Surface States on a Finite Square Photonic Crystal Bounded by Air

Chiral, topologically protected, photonic surface states can be found at the boundary between gyrotropic photonic crystals where a changing magnetic field induces different topology across the interface. Typically, photonic crystals with either a suitable band structure on both sides of the interface to provide a band gap and evanescent decay of the surface states away from the interface, or an outer layer with engineered material properties is required. In this paper, we show the emergence of topological, unidirectional surface states at the termination of finite gyrotropic photonic crystals with a simple square lattice and C 4 rotational symmetry bounded by a vacuum, eliminating the need for an outside layer to enable chiral surface modes. Here, we start from an infinite, time-reversal-symmetry-breaking photonic crystal with a band gap associated with bands with nonzero Chern numbers, different from all-zero Chern numbers in air. We then modify the photonic crystal to move this band gap below the light line, while maintaining the Chern-number discontinuities. Band-structure calculations for a supercell approximating a photonic crystal finite in the direction normal to the surface demonstrate the existence, dispersion, and chirality of the surface mode. Extensive direct scattering calculations for a point source and spatial Fourier analysis further reveal a unidirectional free-space topological surface state, which propagates counterclockwise around the surface of a finite photonic crystal, providing a nearly foolproof way to cross-check the surface-mode band structure unaffected by backscattering from local defects. Additionally, scattering simulations allow an independent characterization of the state dispersion and unveil the robustness of the topological plasmonic mode propagation around the 90° bends of the structure, being due to only radiation leakage. In contrast to buried topological surface states, the observed surface modes at the photonic crystal–air interface have the advantage of being accessible to the outside world, allowing one to take advantage of the defect-tolerant backscattering-free surface modes to engineer emission from photonic crystal surfaces into arbitrary free-space beam shapes and directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A closer look at how symmetry constraints and the spin–orbit coupling shape the electronic structure of Bi(111)

Abstract Fully relativistic density-functional-theory calculations of Bi(111) thin films are analyzed to revisit their two metallic surface-states branches. We first contrast these metallic branches with surface states arising at gaps in the valence band opened by the spin–orbit coupling (SOC). We find that the two metallic branches along Γ M ‾ do not overlap with the bulk band at the zone boundary,M. We show that the spin texture observed in such states cannot be traced to the lifting of Kramers’ degeneracy. Instead, we track them to the m j = ± 1 / 2 – m j = ± 3 / 2 SOC splitting, the potential anisotropy for in-plane and out-of-plane states, and the coupling between the opposite surfaces of a slab occurring nearM, which is driven by a spatial redistribution of the four metallic states composing the two metallic branches. Each of these branches appears to be non-degenerate at the tested surface, yet each is degenerate with another state of opposite spin at the other surface. Nevertheless, the four metallic states bear some contribution on both surfaces of the film because of their spatial redistribution nearM. The overlapping among these states nearM, afforded by their spatial redistribution on both surfaces, causes a hybridization that perpetuates the splitting between the two branches, makes the film’s electronic structure thickness dependent nearM, extinguishes the magnetic moment of the metallic states avoiding the magnetic-moment discontinuity atM, and denies the need or expectancy of the metallic branches becoming degenerate atM. We propose that theoppositespin polarization observed for the two metallic branches occurs because the surface atoms retain their covalent bonds and thus cannot afford magnetic polarization. We show that the Rashba-splitting of the metallic states for inversion-asymmetric films does not have a fixed magnitude but can be tuned by changing the perturbation breaking inversion symmetry.

Physics↗
DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2022

Rapid manufacturing of tailored preforms

A method for rapid manufacturing of three dimensional discontinuous fiber preforms is provided. The method includes the deposition of a polymeric material containing fibers on a surface to form a tailored charge for compression molding. The reinforced polymeric material may be a thermoplastic or a reactive polymer with viscosity low enough to allow flow through an orifice during deposition, yet high enough zero shear viscosity to retain the approximate shape of the deposited charge. The material can be deposited in a predetermined pattern to induce the desired mechanical properties through alignment of the fibers. This deposition can be performed in a single layer or in multiple layers. The alignment is achieved passively by shear alignment of the fibers or actively through fiber orientation control or mixing. The fibers can be of the desired material, length, and morphology, including short and long filaments.

Kunc, Vlastimil↗
DOE OSTI2024

Subsidence Measurements at the Teller 27 Field Site, Seward Peninsula, Alaska, 2021

In September 2021, differential GPS (dGPS) points were collected in areas of visual subsidence in the Teller watershed near mile marker 27 of the Bob Blodgett Nome-Teller Memorial Highway on the Seward Peninsula of Alaska. This study aimed to understand how structural permafrost loss is affecting meter-scale ground elevation. As ice thaws, heaving and slumping of the landscape creates new microtopographical features on the landscape that may alter surface hydrology and soil moisture and thus plant community composition and biogeochemical cycling. The Teller 27 field site is underlain with discontinuous permafrost, so the landscape features a variety of permafrost features. We collected ground elevations from eight subsidence areas in two forms: 1) Transects were sampled laterally through areas of ground subsidence to capture the high edges and low spots where the earth slumped. 2) The circumference of the subsidence area was sampled in order to measure the ground area affected by subsidence. These data were collected with an Emlid Reach RS2 dGPS and a base station. This dataset includes one *.csv of dGPS locations of the eight subsidence locations and one *.kml of measurement locations.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) 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 OSTI2023

Microstructure development and properties of micro-alloyed copper, Cu-0.3Zr-0.15Ag, produced by electron beam additive manufacturing

A micro-alloyed copper powder, Cu-0.3Zr-0.15Ag wt.%, was produced using gas atomization reaction synthesis. Zirconium was added to copper to sequester the oxygen present as copper oxide surface films on the powder particles. The as-received powders, as well as the intentionally oxidized powders were used to fabricate solid test articles by electron beam powder bed fusion additive manufacturing. Dense samples fabricated from as-received powder demonstrated nominal UTS, yield, and elongation values at 260 MPa, 150 MPa, and 34%, respectively. The average electrical conductivity of these samples was measured at 95% of the international annealed copper standard (IACS). Samples fabricated from the oxidized powder exhibited nominal UTS, yield, and elongation of 241 MPa, 146 MPa, and 43%, respectively, with an electrical conductivity of 95 % IACS. During characterization, it was observed that, rather than forming nano-scale dispersoids, the Zirconia (ZrO 2 ) appeared as discontinuous stringers in the metallographic cross-sections that crossed grain and melt pool boundaries. This was rationalized by tracing the presence of the micro-alloying addition of elemental zirconium, which was found to react with surface oxides dissociated in the melt pool to form ZrO 2 , which then solidified on the surface of the melt pool through an allotropic transformation to monoclinic ZrO 2 in discontinuous films and spheroids ranging in size from nanometers to microns. This was confirmed by microscopic analysis of the tops of the melt pools. On subsequent melt passes, these ZrO 2 structures were displaced and redistributed within the melt pool.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Graphic contrastive learning analyses of discontinuous molecular dynamics simulations: Study of protein folding upon adsorption

A comprehensive understanding of the interfacial behaviors of biomolecules holds great significance in the development of biomaterials and biosensing technologies. In this work, we used discontinuous molecular dynamics (DMD) simulations and graphic contrastive learning analysis to study the adsorption of ubiquitin protein on a graphene surface. Our high-throughput DMD simulations can explore the whole protein adsorption process including the protein structural evolution with sufficient accuracy. Contrastive learning was employed to train a protein contact map feature extractor aiming at generating contact map feature vectors. Subsequently, these features were grouped using the k-means clustering algorithm to identify the protein structural transition stages throughout the adsorption process. The machine learning analysis can illustrate the dynamics of protein structural changes, including the pathway and the rate-limiting step. Our study indicated that the protein–graphene surface hydrophobic interactions and the π–π stacking were crucial to the seven-stage adsorption process. Upon adsorption, the secondary structure and tertiary structure of ubiquitin disintegrated. The unfolding stages obtained by contrastive learning-based algorithm were not only consistent with the detailed analyses of protein structures but also provided more hidden information about the transition states and pathway of protein adsorption process and structural dynamics. Our combination of efficient DMD simulations and machine learning analysis could be a valuable approach to studying the interfacial behaviors of biomolecules.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Critical nuclei at hetero-phase interfaces

Two-step nucleation, in which a metastable intermediate phase acts as a precursor for nucleating a thermodynamically stable phase, has been widely observed in many materials systems and solid-state reactions. Among the advantages of two-step nucleation is that the stable phase may nucleate heterogeneously at the hetero-phase interface between the original and the precursory phases. Although heterogeneous nucleation (HN) theories for homo-phase grain boundaries and inert surfaces are well established, our understanding of HN at reactive hetero-phase interfaces remains incomplete. This deficiency stems from the discontinuity of the chemical potential driving force across the hetero-phase interface, which profoundly affects the fundamental properties of the nucleus in a way that is not properly accounted for in existing models. In this work, we incorporate these effects to extend the classical nucleation theory to HN at hetero-phase interfaces. Our extended model demonstrates that the nucleus shape along the minimum energy path is strongly size-dependent, and this additional degree of freedom can result in the reduction of the critical nucleus volume and associated activation energy barrier by orders of magnitude relative to conventional predictions. The simulation results are used to construct a sensitivity map in the parameter space of interfacial energy and bulk driving force ratios, which quantifies the difference in nucleation barriers predicted by different models.

36 MATERIALS SCIENCE↗