Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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

Advancing Multi-Hazard Risk and Safety Considerations for Aging Nuclear Facilities (Revision 1)

This project will demonstrate a multi-hazard time-dependent probabilistic risk assessment (PRA) approach for nuclear facilities considering aging-related deterioration of structures. A generic pressurized water reactor (PWR) reactor subjected to seismic mainshock-aftershock sequences considering the aging of the containment structure will be used as a case study to demonstrate the multi-hazard PRA approach. Using advanced modeling and simulation, seismic mainshock-aftershock fragility functions will be simulated for the containment structure considering aging effects. A multi-hazard PRA model for a generic PWR reactor will be built to quantify the multi-hazard core damage frequency (CDF) and large early release frequency (LERF) with explicit time-dependent modeling of event sequences. To date, the cascading impacts of multi-hazards are not adequately accounted for in the PRA models for nuclear facilities. In addition, for both initial and periodic evaluation of facilities to withstand natural phenomena hazards (NPH), deterioration of the SSCs due to aging and other effects is not adequately considered. By advancing multi-hazard considerations accounting for aging effects, this project will contribute to improved understanding of the safety of aging nuclear facilities. Project outcomes such as the multi-hazard CDF and LERF will also allow facility owners to optimize upgrade and retrofit protocols. This project is divided into two components: (1) advanced modeling and simulations; and (2) multi-hazard PRA. For the first component, Idaho National Laboratory’s (INL) Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) and BlackBear codes will be used to simulate the seismic response and damage of a containment structure under cascading mainshock-aftershock sequences with aging effects. Uncertainties related to the seismic inputs, material parameters, and environmental factors such as temperature and humidity will be identified and propagated to the fragility functions. These fragility functions, which consider aging effects, will be time dependent. The second component will take the fragility information from the first component to build a multi-hazard time-dependent PRA model starting from a generic PWR model to evaluate the multi-hazard CDF and LERF and characterize the associated consequences. For this component, OpenPRA Web Application and SAPHIRE code will be used. Events such as the Fukushima Daiichi accident have highlighted the importance of considering the cascading impacts of multi-hazards for PRA. Moreover, many of the reactors in the current nuclear fleet in the United States (US) are already operating well beyond their initially planned design life, and applications for further extensions to operating licenses are being considered. Therefore, considering multi-hazard effects and aging deterioration in the NPH risk assessment process will contribute to the safety of both existing and future nuclear facilities. Outcomes of this project will thus directly benefit the standards DOE-STD-1020-16 and DOE-HDBK-1224-18.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Development and Verification of the Soil-Pile Interaction Extension for SubDyn

SubDyn is the substructure structural-dynamics module for the aero-hydro-servo-elastic tool FAST v8. SubDyn uses a finite-element model (FEM) to simulate complex multimember lattice structures connected to conventional turbines and towers, and it can make use of the Craig-Bampton model reduction. Here we describe the newly added capability to handle soil-pile stiffness and compare results for monopile and jacket-based offshore wind turbines as obtained with FAST v8, SACS, and EDP (the latter two are modeling software packages commonly used in the offshore oil and gas industry). The level of agreement in terms of modal properties and loads for the entire offshore wind turbine components is excellent, thus allowing SubDyn and FAST v8 to accurately simulate offshore wind turbines on fixed-bottom structures and accounting for the effect of soil dynamics, thus reducing risk to the project.

17 WIND ENERGY↗
DOE OSTI2018

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO 2 -16.0B 2 O 3 -12.6Na 2 O-3.8Al 2 O 3 -5.7CaO-1.7ZrO 2 . Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additional short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO 4 ]-[SiO 4 ] and [BO 4 ]-[SiO 4 ] connections are the most abundant while [ZrO 6 ] octohedra having higher probability to connect to another [SiO 4 ], [BO 4 ] and [BO 3 ]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Structural features of sodium silicate glasses from reactive force field-based molecular dynamics simulations

Atomistic computer simulations can offer insights into silicate glass-environment interactions with the recent development of reactive potentials. Yet, the accuracy of generated glass structures with these potential was usually not fully examined. In this paper, the capability of the reactive force field (ReaxFF) to describe the short and medium range structure features of sodium silicate glasses in molecular dynamics simulations is investigated by comparing a widely used partial charge pairwise potential and available experimental data. Glass structure information such as pair distribution function (PDF), coordination number, Qn species, neutron broadened structure factor, and X-ray broadened structure factor of the glass structures from ReaxFF simulations were calculated and compared to evaluate the generated glass structure. Advantages and limitations of the potentials and glass forming procedures, as well as areas of further improvement, were discussed. The results show that the recently refined ReaxFF parameters through the proposed procedure enable the simulations of sodium silicate glass structures with minimal defects, which paves the way to investigate water-glass interaction mechanisms with the reactive enabled potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Resonant interaction between phonons and PbTe/PbSe (001) misfit dislocation networks

Here this work aims at a quantitative and mechanistic understanding of the dynamic process of the phonon-dislocation interaction in PbTe/PbSe (001) heterostructures using the Concurrent Atomistic-Continuum (CAC) method as the simulation tool. The misfit dislocation network and the atomic-scale dislocation core structure obtained in the simulations are found to agree reasonably well with the experimental observations of the PbTe/PbSe (001) interface. Through visualizing the dynamic interaction between phonons and dislocations, as well as quantifying the dislocation vibration amplitude, the phonon energy transmission, and the thermal resistance of the misfit interfaces, this work has illustrated and quantified two mechanisms for phonon-dislocation interaction: (1) phonon scattering by the strain field of dislocations, and (2) phonon scattering by dislocations that vibrate via the local modes of a dislocation network; the latter, leads to resonant phonon-dislocation interaction, which is manifested as local maxima of out-of-phase vibration of the atoms on the two sides of the slip plane, leading to local minima of the energy transmission in the heterostructure that contains one interface. The local vibrational modes are found to be excited only by shear stress induced by transverse phonons. Among various resonant modes, the one with the lowest frequency has the strongest effect. This work has also demonstrated the collective motion of dislocations under ultrafast phonon pulses. In addition, the dynamic properties of the misfit dislocation network localized within one interface are found to be significantly altered by the presence of misfit dislocations at other interfaces, thus further confirming the cooperative dynamic nature of the motion of dislocations and phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2018

One-Pot Pyrolysis Method to Fabricate Carbon Nanotube Supported Ni Single-Atom Catalysts with Ultrahigh Loading

The practical application of single atom catalysts (SACs) is constrained by the low achievable loading of single metal atoms. In this work, nickel SACs stabilized on a nitrogen-doped carbon nanotube structure (NiSA-N-CNT) with ultrahigh Ni atomic loading up to 20.3 wt % have been successfully synthesized using a new one-pot pyrolysis method employing Ni acetylacetonate (Ni(acac) 2 ) and dicyandiamide (DCD) as precursors. The yield and formation of NiSA-N-CNT depends strongly on the Ni(acac) 2 /DCD ratio and annealing temperature. Pyrolysis at 350 and 650 °C led to the formation of Ni single atom dispersed melem and graphitic carbon nitride (Ni-melem and Ni-g-C 3 N 4 ). Transition from a stacked and layered Ni-g-C 3 N 4 structure to a bamboo-shaped tubular NiSA-N-CNT structure most likely occurs via a solid-to-solid curling or rolling-up mechanism, thermally activated at temperatures of 700–900 °C. Extended X-ray absorption fine structure (EXAFS) experiments and simulations show that Ni single atoms are stabilized in the N-CNT structure through nitrogen coordination, forming a structure with four nearest N coordination shell surrounded by two carbon shells, Ni–N 4 . The NiSA-N-CNT catalysts show an excellent activity and selectivity for the electrochemical reduction of CO 2 , achieving a turnover frequency (TOF) of 11.7 s–1 at -0.55 V (vs RHE), but a low activity for the O 2 reduction and O 2 evolution reactions, as compared to Ni nanoparticles supported on N-CNTs.

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

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Influence of Blood Rheology and Turbulence Models in the Numerical Simulation of Aneurysms

An aneurysm is a vascular malformation that can be classified according to its location (cerebral, aortic) or shape (saccular, fusiform, and mycotic). Recently, the study of blood flow interaction with aneurysms has gained attention from physicians and engineers. Shear stresses, oscillatory shear index (OSI), gradient oscillatory number (GON), and residence time have been used as variables to describe the hemodynamics as well as the origin and evolution of aneurysms. However, the causes and hemodynamic conditions that promote their growth are still under debate. The present work presents numerical simulations of three types of aneurysms: two aortic and one cerebral. Simulation results showed that the blood rheology is not relevant for aortic aneurysms. However, for the cerebral aneurysm case, blood rheology could play a relevant role in the hemodynamics. The evaluated turbulence models showed equivalent results in both cases. Lastly, a simulation considering the fluid–structure interaction (FSI) showed that this phenomenon is the dominant factor for aneurysm simulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Molecular-Scale Structure of Electrode–Electrolyte Interfaces: The Case of Platinum in Aqueous Sulfuric Acid

Knowledge of the molecular composition and electronic structure of electrified solid-liquid interfaces is key to understanding elemental processes in heterogeneous reactions. Using X-ray absorption spectroscopy in the interface-sensitive electron yield mode (EY-XAS), first-principles electronic structure calculations, and multiscale simulations, we determined the chemical composition of the interfacial region of a polycrystalline platinum electrode in contact with aqueous sulfuric acid solution at potentials between the hydrogen and oxygen evolution reactions. We found that between 0.7 and 1.3 V vs Ag/AgCl the electrical double layer (EDL) region comprises adsorbed sulfate ions with hydrated hydronium ions in the next layer. No evidence was found for bisulfate or Pt-O/Pt-OH species, which have very distinctive spectral signatures. In addition to resolving the long-standing issue of the EDL structure, our work establishes interface- and element-sensitive EY-XAS as a powerful spectroscopic tool for studying condensed phase, buried solid-liquid interfaces relevant to various electrochemical processes and devices.

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

Increasing thermostability of the key photorespiratory enzyme glycerate 3‐kinase by structure‐based recombination

As global temperatures rise, improving crop yields will require enhancing the thermotolerance of crops. One approach for improving thermotolerance is using bioengineering to increase the thermostability of enzymes catalysing essential biological processes. Photorespiration is an essential recycling process in plants that is integral to photosynthesis and crop growth. The enzymes of photorespiration are targets for enhancing plant thermotolerance as this pathway limits carbon fixation at elevated temperatures. We explored the effects of temperature on the activity of the photorespiratory enzyme glycerate kinase (GLYK) from various organisms and the homologue from the thermophilic alga Cyanidioschyzon merolae was more thermotolerant than those from mesophilic plants, including Arabidopsis thaliana. To understand enzyme features underlying the thermotolerance of C. merolae GLYK (CmGLYK), we performed molecular dynamics simulations using AlphaFold-predicted structures, which revealed greater movement of loop regions of mesophilic plant GLYKs at higher temperatures compared to CmGLYK. Based on these simulations, hybrid proteins were produced and analysed. These hybrid enzymes contained loop regions from CmGLYK replacing the most mobile corresponding loops of AtGLYK. Two of these hybrid enzymes had enhanced thermostability, with melting temperatures increased by 6 °C. One hybrid with three grafted loops maintained higher activity at elevated temperatures. Whilst this hybrid enzyme exhibited enhanced thermostability and a similar Km for ATP compared to AtGLYK, its Km for glycerate increased threefold. This study demonstrates that molecular dynamics simulation-guided structure-based recombination offers a promising strategy for enhancing the thermostability of other plant enzymes with possible application to increasing the thermotolerance of plants under warming climates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal↗
DOE OSTI2018

Simulation Analysis and Performance Study of CoCrMo Porous Structure Manufactured by Selective Laser Melting

To fabricate porous implants with improved biocompatibility and mechanical properties that are matched to their application using selective laser melting (SLM), flow within the mold and compressive properties and performance of the porous structures must be comprehensively studied. Parametric modeling was used to build 3D models of octahedron and hexahedron structures. Finite element analysis was used to evaluate the mold flow and compressive properties of the parametric porous structures. A DiMetal-100 SLM molding apparatus was used to manufacture the porous structures and the results evaluated by light microscopy. The results showed that parametric modeling can produce robust models. Square structures caused higher blood cell adhesion than cylindrical structures. “Vortex” flow in square structures resulted in chaotic distribution of blood elements, whereas they were mostly distributed around the connecting parts in the cylindrical structures. No significant difference in elastic moduli or compressive strength was observed in square and cylindrical porous structures of identical characteristics. Hexahedron, square and cylindrical porous structures had the same stress–strain properties. For octahedron porous structures, cylindrical structures had higher stress–strain properties. Using these modeling and molding results, an important basis for designing and the direct manufacture of fixed biological implants is provided.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2018

Cold Commissioning Simulant Selection for the Hanford Tank Waste Treatment Plant LAW Facility - 18465

Commissioning the Low-Activity Waste (LAW) Facility of the Hanford Tank Waste Treatment and Immobilization Plant (WTP) will require large volumes (∼2,000,000 liters) of a non-radioactive, chemically analogous waste simulant. Cold commissioning will provide for operator training in an integrated facility, demonstration of mechanical process equipment capacity, environmental performance demonstration of emissions control systems, and process control systems shakedown and tuning. Selection of the cold commissioning simulant is following a structured four-phase process to arrive at the preferred composition and recipe. Phases I, II and III are complete. Phase I laid the framework for decision making with preparation of a Project Management Plan, an update to the Simulant Development, Approval, Validation, and Documentation procedure, and a kickoff meeting to gather stakeholder input into the simulant definition and selection criteria. Phase II consisted of a comprehensive literature search and compilation of LAW simulant recipes, followed by initial screening against the selection criteria identified in Phase I. This process identified 72 candidate simulant recipes and down-selected 9 for further study. Phase III studies involved a) statistical comparisons of the similarity of candidate simulant compositions to the tank waste feed vector, b) process modelling of the commissioning sequence with each simulant to estimate glass former and simulant material consumption, offgas and effluent volumes and compositions, and product glass volumes and compositions, and c) evaluation of simulants against the immobilized LAW glass composition algorithm to determine allo able waste loading to ensure compliance with product specifications and process constraints. The results of these studies identified two of the candidate recipes to be carried forward to full-scale testing. Phase IV consists of testing these two recipes at vendor facilities to determine scalability, fabricability, and stability of the simulants when produced at industrial scale. Scalability is the ability to take a recipe from bench scale (1 to 40 liters) batch size to full scale (>35,000 liters) batch size. Fabricability is the ability of a recipe to be prepared in large batches without the need for complex or expensive mechanical equipment, heating, or cooling systems. Stability is the ability of a simulant to withstand long term storage and variations in environmental conditions without unacceptable changes in key characteristics (the set of chemical, physical, and rheological parameters that the simulant is intended to exhibit during cold commissioning tests). Simulants will also be tested for compatibility with other startup material (i.e., frit), tuning simulants, and environmental performance demonstration test spiking solutions to ensure against unexpected situations. Using a structured process is essential when making high impact decisions such as waste simulant selection. Proper planning and a clear explanation of how the decision process will unfold is key to obtaining stakeholder buy-in as more information is gathered. The documentation resulting from this structured selection process will be compiled and summarized in a Cold Commissioning Simulant Definition deliverable due to the Department of Energy 24 months before the start of Cold Commissioning. (authors)

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

Connecting Structural and Transport Properties of Ionic Liquids with Cationic Oligoether Chains

We used X-ray diffraction and molecular dynamics simulations to probe the structures of two families of ionic liquids containing oligoether tails on the cations. Imidazolium and pyrrolidinium bis(trifluoromethylsulfonyl)amide ILs with side chains ranging from 4 to 10 atoms in length, including both linear alkyl and oligo-ethylene oxide tails, were prepared. Furthermore, their physical properties, such as viscosity, conductivity and thermal profile, were measured and compared for systematic trends. Consistent with earlier literature, a single ether substituent substantially decreases the viscosity of pyrrolidinium and imidazolium ILs compared to their alkyl congeners. Remarkably, as the number of ether units in the pyrrolidinium ILs increases there is hardly any increase in the viscosity, in contrast to alkylpyrrolidinium ILs where the viscosity increases steadily with chain length. Viscosities of imidazolium ether ILs increase with chain length but always remain well below their alkyl congeners. To complement the experimentally determined properties, molecular dynamics simulations were run on the two ILs with the longest ether chains. Our results point to specific aspects that could be useful for researchers designing ILs for specific applications.

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

Machine learning inversion from small-angle scattering for charged polymers

We develop Monte Carlo simulations for uniformly charged polymers and a machine learning algorithm to interpret the intra-polymer structure factor of the charged polymer system, which can be obtained from small-angle scattering experiments. The polymer is modeled as a chain of fixed-length bonds, where the connected bonds are subject to bending energy, and there is also a screened Coulomb potential for charge interaction between all joints. The bending energy is determined by the intrinsic bending stiffness, and the charge interaction depends on the interaction strength and screening length. All three contribute to the stiffness of the polymer chain and lead to longer and larger polymer conformations. The screening length also introduces a second length scale for the polymer besides the bending persistence length. To obtain the inverse mapping from the structure factor to these polymer conformation and energy-related parameters, we generate a large data set of structure factors by running simulations for a wide range of polymer energy parameters. We use principal component analysis to investigate the intra-polymer structure factors and determine the feasibility of the inversion using the nearest neighbor distance. We employ Gaussian process regression to achieve the inverse mapping and extract the characteristic parameters of polymers from the structure factor with low relative error.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗
DOE OSTI2017

Influence of Asymmetric Cyclic Loading on Structural Evolution and Deformation Behavior of Cu-5 at.% Zr Alloy: An Atomistic Simulation-Based Study

Molecular dynamics (MD) simulation-based studies of tensile test and structural evolution of Cu-5 at.% Zr alloy under asymmetric cyclic loading (i.e., ratcheting behavior) considering various stress ratios such as − 0.2, − 0.4 and − 0.6 for different temperatures, viz.≈ 100, 300 and 500 K have been performed using embedded atom model Finnis–Sinclair potential. According to obtained stress–strain response from MD calculation, Cu-5 at.% Zr alloy specimen is pristine in nature as sudden drop in stress just after yield stress and subsequent elastic type deformation are observed for this alloy. Predicted ratcheting strain by MD simulation for Cu-5 at.% Zr alloy varies from 4.5 to 5%. Significant increase in ratcheting strain has been observed with the increase in temperature. Slight reduction in crystallinity is identified at the middle of the each loading cycle from the performed radial distribution function analysis and cluster analysis.

36 MATERIALS SCIENCE↗