Search NASA⌕ Search

Engineering topics

Gao, Ke

Publications and source records attributed to Gao, Ke.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗
DOE OSTI2022

The Crystalline Behavior and Device Function of Nonfullerene Acceptors in Organic Solar Cells

The current research investigates the structure features and intermolecular interactions of nonfullerene acceptors (NFAs) in single crystal and thin films, as well as their solar cell applications. Guiding parameters and key intermolecular forces that lead to 2D brickwork or 3D web packing are identified. The atomic modification is shown as the key to induce hydrogen bonding or π–π stacking column, which results in different crystalline packing. The molecular assembly in thin film is initiated by hydrogen bonding and completed by π–π stacking reorganization. The packing energy is seen as a guiding parameter that dictates the NFA crystalline morphology in blended thin films. The crystalline packing motif is not directly related with device efficiency. However, the crystalline morphology is the key parameter to influence exciton/carrier dynamics and device performance. In conclusion, a broader picture on the scaling behavior of organic semiconductor crystals ranging from oligoacenes to NFAs is established.

14 SOLAR ENERGY↗
DOE OSTI2022

The Molecular Ordering and Double-Channel Carrier Generation of Nonfullerene Photovoltaics within Multi-Length-Scale Morphology

The success of nonfullerene acceptor (NFA) solar cells lies in their unique physical properties beyond the extended absorption and suitable energy levels. The current study investigates the morphology and photophysical behavior of PBDB-T donor blending with ITIC, 4TIC, and 6TIC acceptors. Single-crystal study shows that the π–π stacking and side-chain interaction dictate molecular assembly, which can be carried to blended films, forming a multi-length-scale morphology. Spontaneous carrier generation is seen in ITIC, 4TIC, and 6TIC neat films and their blended thin films using the PBDB-T donor, providing a new avenue of zero-energy-loss carrier formation. The molecular packing associated with specific contacts and geometry is key in influencing the photophysics, as demonstrated by the charge transfer and carrier lifetime results. The 2D layer of 6TIC facilitates the exciton-to-polaron conversion, and the largest photogenerated polaron yield is obtained. Finally, the new mechanism, together with the highly efficient blending region carrier generation, has the prospect of the fundamental advantage for NFA solar cells, from molecular assembly to thin-film morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The coupling and competition of crystallization and phase separation, correlating thermodynamics and kinetics in OPV morphology and performances

The active layer morphology transition of organic photovoltaics under non-equilibrium conditions are of vital importance in determining the device power conversion efficiency and stability; however, a general and unified picture on this issue has not been well addressed. Using combined in situ and ex situ morphology characterizations, morphological parameters relating to kinetics and thermodynamics of morphology evolution are extracted and studied in model systems under thermal annealing. The coupling and competition of crystallization and demixing are found to be critical in morphology evolution, phase purification and interfacial orientation. A unified model summarizing different phase diagrams and all possible kinetic routes is proposed. The current observations address the fundamental issues underlying the formation of the complex multi-length scale morphology in bulk heterojunction blends and provide useful morphology optimization guidelines for processing devices with higher efficiency and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2020

IC Project: w19_hossfault “Modelling of stick-slip behavior in sheared granular fault gouge & nonlinear elasticity behavior in cracked solid”

Low frequency earthquakes, non-volcanic tremor, and acoustic emissions are examples of weak seismic signals that may help detect major seismic events, i.e., earthquakes. In a laboratory setting, in which stick-slip events are simulated, acoustic emissions are detected far from the stick-slip events. In both, field and laboratory scale cases, the acoustic emission signals are sourced or detected in a volume remote from the volume that spawns the earthquake. Therefore, it is relevant to establish the causal relationship between signals detected on passive, remote monitors and the dynamics of the elastic structures that launch important seismic events. The work conducted under this IC allocation allowed us to utilize a numerical model that let us follow this causality, i.e., examine and connect the dynamics in a granular system (fault gouge), to signals detected on passive remote monitors. It was demonstrated that stress chains are key in the dynamics of granular systems and that their evolution is the source of the acoustic emission.

58 GEOSCIENCES↗
DOE OSTI2020

As‐Cast Ternary Organic Solar Cells Based on an Asymmetric Side‐Chains Featured Acceptor with Reduced Voltage Loss and 14.0% Efficiency

Abstract A new small‐molecule nonfullerene acceptor based on the benzo[1,2‐ b :4,5‐ b ′]dithiophene (BDT) fused central core with asymmetrical alkoxy and thienyl side chains, namely TOBDT , is designed and synthesized. The alkoxy unit helps narrow the bandgap, and thienyl side chain helps enhance the intermolecular interaction. As a result, TOBDT is suitable to match the deep‐lying highest occupied molecular orbital (HOMO) of polymer donor PM6 . Then, a strong crystalline acceptor IDIC is introduced as the third component to fabricate as‐cast nonfullerene ternary devices to achieve absorption and morphology control. Addition of IDIC not only mixes well with TOBDT but modulates the morphology of the blend film, which helps to balance the charge transport properties and reduce the photovoltage loss of ternary devices. All these contribute to synergetic improvement of J sc , V oc , and fill factor parameters, leading to a power conversion efficiency of 14.0% for the as‐cast fullerene‐free ternary device.

Chen, Xuebin↗