Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface capturing”

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 343 records · Page 19

Visual Style Elements in Human–System Interface Design for Nuclear Power Operations: Does Style Affect Performance and Preference?

As the world increasingly adopts renewable and sustainable energy systems, transitionary solutions include nuclear power, which currently provides 20% of the United States’ electricity and is the largest single source of carbon-free electricity generation. Advanced reactors are a critical component of a carbon-free mixed energy portfolio that require careful design of first-of-a-kind control rooms. The application of Human Factors Engineering (HFE) is essential for scientific and iterative testing of novel human–system interface (HSI) concepts to ensure effective, efficient, and safe plant operations. Microworlds are simulators that use simplified physics models and control systems to distill nuclear power operations into essential functions. HFE scientists used the Rancor Microworld Simulator to obtain preference and performance metrics for novel and traditional static HSI design styles. Participants comprised advanced reactor company employees and nuclear industry consultants. A mixture of quantitative and qualitative data was captured. There was a preference for the basic graphical style that included high contrast and traditional color scheme elements. No single HSI design outperformed the others, and the participants did not perform better using their preferred HSI style. We report this experiment is the first in a series of HFE testing for HSIs in advanced reactor control room development. Clear user preferences emerged for elements within static displays. The cutting-edge neumorphic style was the least preferred. Future directions include tests of dynamic displays. HFE is used in evaluating and designing HSI devices that will improve the efficiency and safety of advanced nuclear power operations.

99 GENERAL AND MISCELLANEOUS↗

Microscopic and Macroscopic Characterization of Grain Boundary Energy and Strength in Silicon Carbide via Machine-Learning Techniques

Predicting the properties of grain boundaries poses a challenge because of the complex relationships between structural and chemical attributes both at the atomic and continuum scales. Grain boundary systems are typically characterized by parameters used to classify local atomic arrangements in order to extract features such as grain boundary energy or grain boundary strength. The present work utilizes a combination of high-throughput atomistic simulations, macroscopic and microscopic descriptors, and machine-learning techniques to characterize the energy and strength of silicon carbide grain boundaries. Additionally, a diverse data set of symmetric tilt and twist grain boundaries are described using macroscopic metrics such as misorientation, the alignment of critical low-index planes, and the Schmid factor, but also in terms of microscopic metrics, by quantifying the local atomic structure and chemistry at the interface. These descriptors are used to create random-forest regression models, allowing for their relative importance to the grain boundary energy and decohesion stress to be better understood. Results show that while the energetics of the grain boundary were best described using the microscopic descriptors, the ability of the macroscopic descriptors to reasonably predict grain boundaries with low energy suggests a link between the crystallographic orientation and the resultant atomic structure that forms at the grain boundary within this regime. For grain boundary strength, neither microscopic nor macroscopic descriptors were able to fully capture the response individually. However, when both descriptor sets were utilized, the decohesion stress of the grain boundary could be accurately predicted. These results highlight the importance of considering both macroscopic and microscopic factors when constructing constitutive models for grain boundary systems, which has significant implications for both understanding the fundamental mechanisms at work and the ability to bridge length scales.

36 MATERIALS SCIENCE↗

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

Materials data science using CRADLE: A distributed, data-centric approach

Abstract There is a paradigm shift towards data-centric AI, where model efficacy relies on quality, unified data. The common research analytics and data lifecycle environment (CRADLE™) is an infrastructure and framework that supports a data-centric paradigm and materials data science at scale through heterogeneous data management, elastic scaling, and accessible interfaces. We demonstrate CRADLE’s capabilities through five materials science studies: phase identification in X-ray diffraction, defect segmentation in X-ray computed tomography, polymer crystallization analysis in atomic force microscopy, feature extraction from additive manufacturing, and geospatial data fusion. CRADLE catalyzes scalable, reproducible insights to transform how data is captured, stored, and analyzed. Graphical abstract

97 MATHEMATICS AND COMPUTING↗

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing Electrolyte Influence on CO 2 Reduction in Aprotic Solvents

Selective CO 2 capture and electrochemical conversion are important tools in the fight against climate change. Industrially, CO 2 is captured using a variety of aprotic solvents due to their high CO 2 solubility. However, most research efforts on electrochemical CO 2 conversion use aqueous media and are plagued by competing hydrogen evolution reaction (HER) from water breakdown. Fortunately, aprotic solvents can circumvent HER, making it important to develop strategies that enable integrated CO 2 capture and conversion. However, the influence of ion solvation and solvent selection within nonaqueous electrolytes for efficient and selective CO 2 reduction is unclear. In this work, we show that the bulk solvation behavior within the nonaqueous electrolyte can control the CO2 reduction reaction and product distribution occurring at the catalyst-electrolyte interface. We study different tetrabutylammonium (TBA) salts in two electrolyte systems with glyme ethers (e.g., 1,2 dimethoxyethane or DME) and dimethyl sulfoxide (DMSO) as a low and high dielectric constant medium, respectively. Using spectroscopic tools, we quantify the fraction of ion pairs that forms within the electrolyte. Also, we show how ion pair formation is prevalent in DME and is dependent on the anion type. More importantly, we show that as ion pair formation decreases within the electrolyte, CO 2 current densities increase, and a higher CO Faradaic efficiency is observed at low overpotentials. Meanwhile, in an electrolyte medium where the ion pair fraction does not change with the anion type (such as in DMSO), a smaller influence of solvation is observed on CO 2 current densities and product distribution. By directly coupling bulk solvation to interfacial reactions and product distribution, we showcase the importance and utility of controlling the reaction microenvironment in tuning the electrocatalytic reaction pathways. Insights gained from this work will enable novel electrolyte designs for efficient and selective CO 2 conversion to desired fuels and chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Facile hermetic TEM grid preparation for molecular imaging of hydrated biological samples at room temperature

Abstract Although structures of vitrified supramolecular complexes have been determined at near-atomic resolution, elucidating in situ molecular structure in living cells remains a challenge. Here, we report a straightforward liquid cell technique, originally developed for real-time visualization of dynamics at a liquid-gas interface using transmission electron microscopy, to image wet biological samples. Due to the scattering effects from the liquid phase, the micrographs display an amplitude contrast comparable to that observed in negatively stained samples. We succeed in resolving subunits within the protein complex GroEL imaged in a buffer solution at room temperature. Additionally, we capture various stages of virus cell entry, a process for which only sparse structural data exists due to their transient nature. To scrutinize the morphological details further, we used individual particle electron tomography for 3D reconstruction of each virus. These findings showcase this approach potential as an efficient, cost-effective complement to other microscopy technique in addressing biological questions at the molecular level.

59 BASIC BIOLOGICAL SCIENCES↗

Improving energy yield in photovoltaic modules with photonic structures

When installed outdoors, Si solar cells typically operate 20 – 30 K above ambient conditions. These elevated operating temperatures lead to diminished performance, reducing efficiency by ~0.4%/K for crystalline Si cells and decreasing module lifetime. Since much of the elevated temperature of the modules derives from parasitic absorption of sub-bandgap sunlight by the contacts, encapsulants, or other materials, reflectors that remove the sub-bandgap radiation from the module before it is absorbed promise to reduce the operating temperature. However, there are many competing strategies to reject this light: selective reflectors could be integrated into the outer module glass, at the textured Si-encapsulant interface, or on the rear. Furthermore, many of these approaches are complex and expensive to fabricate. The issue of heating can be even more critical in bifacial modules: high fractions of rear irradiance can lead to higher operating temperatures. Prior to this project, relatively little was known about the optical strategies for temperature reduction in bifacial modules. The goal of this project was to quantitatively assess optical strategies for improving the energy yield via a combination of improved anti-reflection and sub-bandgap light rejection. We studied two complementary concepts: the performance of ideal structures to determine the upper limits to performance, and low-complexity structures made from standard materials as cost-effective implementations. We also studied the performance of these reflectors in different types of modules, including Al BSF, PERC, and bifacial. Here we build from our prior accomplishments in developing mirror optimization techniques, and make use of simulation methods that accurately capture the optical properties of the modules throughout the visible and near-infrared spectrum. We connect these optical properties to energy yield under realistic, outdoor conditions. We also experimentally fabricated mirrors, incorporated them into mini-modules, and performed outdoor testing of performance at NREL. We found definitively that spectrally-selective mirrors on the outside of the module cover glass offer the best performance, both under ideal conditions and for realistic, low- complexity mirrors. Mirrors at the Si cell/encapsulant interface are limited by the multiple reflections needed to remove sub-bandgap light. Mirrors at the cell surface and mirrors on the rear also cannot remove all the sub-bandgap light, since parasitic absorption also occurs in the cover glass, encapsulant, and front contacts, depending on the module configuration. Also, cell texturing places higher demands on any cell-surface mirror because multiple interactions with the mirror are required to reject sub bandgap light which compounds losses. For bifacial modules, mirrors on the back do provide additional cooling and energy yield improvements, and the mirrors can have the same structure. In the monofacial case up to 3.3K of cooling is possible for mirrors on the front, compared to 2.2K at the cell surface and 1.2K at the rear. For bifacial modules, up to 2.4K of cooling is possible when mirrors are included at both the front and back. Our realistic designs utilize 4 – 6 layers of common photovoltaic materials such as SiO 2 , TiO 2 , and SiN x , making them more cost effective. Finally, we used our models to study thermal management in bifacial modules with different back lamination materials, finding that glass-glass laminated modules operated hotter than equivalent glass-polymer backsheet modules, due solely to the increased absorption of light incident from the rear of the module. Improving the energy yield by 2%, which is achievable with these realistic designs, would lead to LCOE reductions of $0.03/kWh by 2030. These benefits make photovoltaic panels more accessible to the public and may reduce associated maintenance costs.

14 SOLAR ENERGY↗

Developing Flexible, Networked Lighting Control Systems That Reliably Save Energy in California Buildings

An important strategy to meet California's ambitious energy efficiency goals is to use innovative wireless communications, embedded sensors, data analytics and controls to significantly reduce lighting energy use in commercial buildings. This project developed a suite of networked lighting solutions to further this goal. The technologies include a platform for low-cost sensing, distributed intelligence and communications, the “PermaMote,” which is a self-powered sensor and controller for lighting applications. The project team also developed a task ambient daylighting system that integrates sensors with data-driven daylighting control using an open communication interface, called the “Readings-At-Desk” (RAD) system. To address the problem of building occupants being confused about how to operate traditional lighting control systems, the research team created content that could be the basis for a user interface standard for lighting controls. Finally, to address the difficulty of ensuring that advanced lighting control systems actually deliver their promised energy savings, the project team developed a new method for evaluating and specifying lighting systems’ performance. The research team validated these technologies in the laboratory, showing significant lighting energy savings, up to 73% for the PermaMote sensor system from occupancy control and daylight dimming features, compared to the same light source (LED replacement lamps) operated via simple on/off scheduling. The project team also developed a proposed standard lighting data model and user interface elements, which were contributed to the ANSI Lighting Systems Committee (C137) for standardization. Existing data models are incomplete and inconsistent, whereas the lighting-specific data model developed here is clear and comprehensive, to serve as a starting point for creating common, universally agreed upon semantic definitions of key lighting parameters, to promote interoperability. For the task on verifiable performance of lighting systems, the project team developed a more effective metric for capturing the actual energy impact of a lighting system over time — the energy usage intensity (kWh/ft2/year). Three commercial lighting systems were tested in FLEXLAB® using this new metric, and the tests show a wide range in the accuracy of the self-reported energy-use metric, from 0.5% to 28% error compared to direct measurement of lighting energy using dedicated submeters. Overall, the project team estimates that these advanced technologies can reduce California office lighting energy use by 20% (above and beyond normal advanced lighting controls mandated by Title 24), resulting in about 1,600 GWh/year in savings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Single-Molecule Detection of the Encounter and Productive Electron Transfer Complexes of a Photosynthetic Reaction Center

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanoconfined Interfaces for Highly Selective Separation of Critical Rare Earth Elements

Industrial demand for rare earth elements (REEs) has surged over the past three decades due to their unique properties that support sustainable energy and new technologies. Separating individual REEs is challenging and hazardous, typically done through liquid-liquid extraction. There is an urgent need for environmentally friendly and efficient separation technologies for REEs. Porous materials offer promising advances for sustainable REE separation via ion-selective capture. We hypothesize that REE separation can be efficiently achieved in reactive nanopores, such as Zr(IV) and Cr(III) metal-organic frameworks (MOFs), through surface functionalization. By integrating material synthesis, interfacial chemistry experiments, theory, computation, and machine learning, we gained insights into the chemical factors controlling REE speciation and their competitive adsorption on MOFs. Our findings show that these materials’ selectivity can be tuned by surface functionalization. The machine learning component addressed ion-specific diffusion based on MOF topology and chemistry.

36 MATERIALS SCIENCE↗

Spin–Phonon Coupling and Magnetic Transition in an Organic Molecule Intercalated Cr 2 Ge 2 Te 6

The manipulation of spin–phonon coupling in both formations and explorations of magnetism in two-dimensional van der Waals ferromagnetic semiconductors facilitates unprecedented prospects for spintronic devices. The interlayer engineering with spin–phonon coupling promises controllable magnetism via organic cation intercalation. Here, in this study, spectroscopic evidence reveals the intercalation effect on the intrinsic magnetic and electronic transitions in quasi-two-dimensional Cr 2 Ge 2 Te 6 using tetrabutyl ammonium (TBA + ) as the intercalant. The temperature evolution of Raman modes, E g 3 and A g 1 , along with the magnetization measurements, unambiguously captures the enhancement of the ferromagnetic Curie temperature in the intercalated heterostructure. Moreover, the E g 4 mode highlights the increased effect of spin–phonon interaction in magnetic-order-induced lattice distortion. Combined with the first-principle calculations, we observed a substantial number of electrons transferred from TBA + to Cr through the interface. The interplay between spin–phonon coupling and magnetic ordering in van der Waals magnets appeals for further understanding of the manipulation of magnetism in layered heterostructures.

2D magnet↗

Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie

Ecosystems at coastal terrestrial–aquatic interfaces play a significant role in global biogeochemical cycles. In this study, we aimed to characterize coastal wetlands with particular focus on the co-variability between plant dynamics, topography, soil, and other environmental factors. We proposed a functional zonation approach based on machine learning clustering to identify the spatial regions, i.e., zones that capture these co-varied properties. This approach was applied to publicly available datasets along Lake Erie, in the Great Lakes Region. We investigated the heterogeneity of coastal ecosystem structures as a function of along-shore distance and transverse distance, based on the spatial data layers, including topography, wetland vegetation cover, and the time series of Landsat’s enhanced vegetation index (EVI) between 1990 and 2020. Results showed that the topographic metrics (elevation and slope), soil texture, and plant productivity influence the spatial distribution of wetland land-covers (emergent and phragmites). These results highlight a natural organization along the transverse axis, where the elevation and the EVI increase further away from the coastline. In addition, the clustering analysis allowed us to identify regions with distinct environmental characteristics, as well as the ones that are more sensitive to interannual lake-level variations.

47 OTHER INSTRUMENTATION↗

Probabilistic Modeling of Commercial Building Occupancy Patterns Using Location-Based Map Data: Preprint

Considering occupancy patterns is crucial to simulate buildings' energy use. Current energy models use inputs that simplify the actual diversity in occupancy into static occupancy patterns and are not able to represent the numerous variations in occupancy patterns between buildings and across different locations. Recently, inferring occupancy schedules from metered electricity consumption data was used to model occupancy in commercial buildings. However, the translation from metered data to occupancy schedules requires many assumptions that might not capture the reality, and the process is hindered by the availability of data from advanced metering infrastructure. With the development of information technologies, occupancy modeling should not be limited to traditional approaches. The prevalence of social networks and location services with real-time user feedback provides publicly accessible data via Maps Application Programming Interfaces (APIs) such as Google Maps, SafeGraph, Mapbox, Foursquare, etc. This paper presents an automated framework for modeling parametric occupancy patterns using such APIs to calibrate commercial district buildings' energy models. This process includes three main steps: data extraction and processing, parametric schedules generation, and schedules integration. We demonstrated this framework in districts where we used maps API to generate more accurate behavioral patterns for operations and electric vehicle charging events. We used these patterns to determine differences in energy use across key sociodemographic and spatial parameters. The presented method has the potential for worldwide applications. Users can utilize this framework to extract data for selected locations of interest to create more realistic behavioral patterns for commercial facilities across different districts.

building energy modeling↗

Synchrotron X-ray based particle image velocimetry to measure multiphase streamflow and densitometry

Synchrotron-based radiation techniques are increasingly used in the field of flow condition measurements. As a complicated multiphase flow, cavitation has been investigated in a wide variety of industrial fields, to prevent damages on spillways at large dams or intensify the bacteria eradication rate during water treatments. To control the detrimental effects of cavitating flows, the dynamics of and mechanisms affecting cavitation development shall be visualized and identified in real-time. Due to the large reflection and scattering effects of multiphase flow, standard visible-light imaging technologies cannot penetrate vapor bubbles and provide extremely limited information. In this study, the synchrotron source emitted high-flux X-ray pulses are used to capture flow motion and visualize internal structures with sufficiently high spatial and temporal resolutions. The image processing procedures on raw x-ray images are developed using 2D Fourier transform (FT) and wavelet transform (WT) to implement the contrast enhancement and de-noising. Through the edge detection algorithm, tracking either seeded particles or phase interfaces inside the opaque multiphase flow can be achieved to perform the particle image velocimetry. Further, with decomposing particle-only and particle-eliminated subplots from the X-ray phase-contrast image, the time-resolved velocity and void fraction fields are obtained simultaneously, paving a way for further flow condition analysis such as densitometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Short-Term Groundwater Level Fluctuations Drive Subsurface Redox Variability

As global change processes modify the extent and functions of terrestrial-aquatic interfaces, the variability of critical and dynamic transitional zones between wetlands and uplands increases. However, it is still unclear how fluctuating water levels at these dynamic boundaries alter groundwater biogeochemical cycling. Here, we used high-temporal resolution data along gradients from wetlands to uplands and during fluctuating water levels at freshwater coastal areas to capture spatiotemporal patterns of groundwater redox potential (E h ). We observed that topography influences groundwater E h that is higher in uplands than wetlands; however, the high variability within TAI zones challenged the establishment of distinct redox zonation. Declining water levels generally decreased E h , but most locations exhibited significant E h variability, which is associated with rare instances of short-term water level fluctuations introducing oxygen. The E h -oxygen relationship showed distinct hysteresis patterns, reflecting redox poising capacity at higher E h , maintaining more oxidizing states longer than dissolved oxygen presence. Surprisingly, we observed more frequent oxidizing states in transitional areas and wetlands than in uplands. We infer that occasional oxygen entering specific wetland-upland boundaries acts as critical biogeochemical control points. Finally, high-resolution data can capture such rare yet significant biogeochemical instances, supporting redox-informed models and advancing the predictability of climate change feedback.

54 ENVIRONMENTAL SCIENCES↗

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗