Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural mechanics”

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 19 records

Modern chemical graph theory

Abstract Graph theory has a long history in chemistry. Yet as the breadth and variety of chemical data is rapidly changing, so too do graph encoding methods and analyses that yield qualitative and quantitative insights. Using illustrative cases within a basic mathematical framework, we showcase modern chemical graph theory's utility in Chemists' analysis and model development toolkit. The encoding of both experimental and simulation data is discussed at various levels of granularity of information. This is followed by a discussion of the two major classes of graph theoretical analyses: identifying connectivity patterns and partitioning methods. Measures, metrics, descriptors, and topological indices are then introduced with an emphasis upon enhancing interpretability and incorporation into physical models. Challenging data cases are described that include strategies for studying time dependence. Throughout, we incorporate recent advancements in computer science and applied mathematics that are propelling chemical graph theory into new domains of chemical study. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Structure and Mechanism > Computational Materials Science Structure and Mechanism > Molecular Structures

Leite, Leonardo S. G.↗

UPDATES FROM THE INVOLUTE WORKING GROUP

The HFIR, RHF, and FRM II reactors represent a particular class of Research and Test Reactors that provide some of the most intense and continuous neutron fluxes for science, industry, and medical applications. These high-performance reactors have achieved compact cores by operating with Highly Enriched Uranium fuel (HEU, 235U/U ≥ 20 wt. %) and utilizing fuel plates curved as an involute. Due to the proliferation risks, the international community aims to reduce or eliminate, when possible, the use of HEU fuel in civilian facilities by converting them to a Low-Enriched Uranium fuel (LEU, 235U/U < 20 wt. %). Conversion of these reactors without significantly compromising their performance or safety is a challenging endeavor that can tremendously benefit from advanced computational tools and thus, eliminate unnecessary conservatism to ensure sufficient thermal margins. Therefore, models are being developed using modern Computational Fluid Dynamics (CFD) and Computational Structural Mechanics (CSM) software to evaluate the steady-state safety margins of various LEU designs instead of being reliant on the more traditional, conservative methods. To gain the confidence and acceptance of high-fidelity modeling by the nuclear regulators, Argonne National Laboratory (ANL) and the involute reactors have formed an informal scientific group, the Involute Working Group (IWG). The IWG facilitates inter-organizational collaboration on experimental benchmarking, code-to-code comparisons, and Verification and Validation (V&V). This paper describes some of the recent IWG efforts in validating software against the existing experimental data, as well as code-to-code comparisons of different software used by the IWG members.

Bergeron, Aurelien↗

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗

Understanding flavin electronic structure and spectra

Abstract Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low‐grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light‐driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Electronic Structure Theory > Combined QM/MM Methods Theoretical and Physical Chemistry > Spectroscopy

59 BASIC BIOLOGICAL SCIENCES↗

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗

First Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the sensor was initially investigated with analytical derivations and with COMSOL RF Module computer simulations of resonant frequency spectral shift due to uniform load. We also investigated the mechanical integrity of the flowmeter’s membrane through analytical modelling and COMSOL Structural Mechanics Module computer simulations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the material’s yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure. Using results from models, we have developed an initial design for a microwave K-band sensor. A cylindrical resonator prototype was fabricated from brass for the initial tests. The external dimensions of the cavity are matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed consisting of a piping Tee with a bulkhead WR-42 microwave waveguide installed in a leak-proof assembly. A microwave waveguide circulator was installed in the setup to suppress the effect of reflections at the cavity entrance by increasing the isolation between the input and the output port. Preliminary spectral characterization of cavity spectral response was performed with a portable PXIe chassis microwave VNA with a custom GUI. Preliminary dry tests of the transducer response were conducted with a set of calibrated weights. Transducer frequency shift was shown to be monotonically increasing with increasing pressure. The next steps will involve investigation of the transducer performance for water flow sensing.

42 ENGINEERING↗

Chemical transformations and transport phenomena at interfaces

Abstract Interfaces, the boundary that separates two or more chemical compositions and/or phases of matter, alters basic chemical and physical properties including the thermodynamics of selectivity, transition states, and pathways of chemical reactions, nucleation events and phase growth, and kinetic barriers and mechanisms for mass transport and heat transport. While progress has been made in advancing more interface‐sensitive experimental approaches, their interpretation requires new theoretical methods and models that in turn can further elaborate on the microscopic physics that make interfacial chemistry so unique compared to the bulk phase. In this review, we describe some of the most recent theoretical efforts in modeling interfaces, and what has been learned about the transport and chemical transformations that occur at the air–liquid and solid–liquid interfaces. This article is categorized under: Structure and Mechanism > Reaction Mechanisms and Catalysis Structure and Mechanism > Computational Materials Science Software > Quantum Chemistry Software > Simulation Methods

Hao, Hongxia↗

Wind Turbine Rotor Design Using High-Fidelity Aerostructural Optimization

Large wind turbines yield more energy but demand careful aeroelastic blade design. Coupled multiphysics design strategies can reduce wind energy costs by exploiting fluid-structure interactions. This work presents the first high-fidelity aerostructural optimization study of a large wind turbine rotor. We use blade-resolved fluid dynamics and structural solvers in a monolithic gradient-based optimization framework to explore steady-state torque and blade mass tradeoffs. The coupled-adjoint approach computes gradients efficiently, enabling the optimization of over 100 structural and geometric parameters simultaneously. Our optimization study modifies a DTU 10 MW benchmark with a simplified structure and isotropic material properties. The tightly coupled optimizations increase torque by 14% while reducing rotor mass by 9% or reduce blade mass by 27% while maintaining torque. Blade-resolved models provide greater design freedom, enabling 5% higher mass reductions than conventional parameterizations at equal torque. This framework paves the way for more detailed high-fidelity optimization studies to complement conventional design approaches.

17 WIND ENERGY↗

Mechanical and Thermal Properties of Low-Density Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y Alloys

Refractory high-entropy alloys (RHEAs) Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y ((x, y) = (0, 0), (0, 10), and (10, 15)) were computationally studied to obtain a low density and a better mechanical property. The density functional theory (DFT) method was employed to compute the structural and mechanical properties of the alloys, based on a large unit cell model of randomly distributed elements. Debye–Grüneisen theory was used to study the thermal properties of Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y . The phase diagram calculation shows that all three RHEAs have a single body-centered cubic (BCC) structure at high temperatures ranging from 1000 K to 2000 K. The RHEA Al 30 Cr 10 Mo 5 Ti 20 V 35 has shown a low density of 5.16 g/cm 3 and a hardness of 5.56 GPa. The studied RHEAs could be potential candidates for high-temperature application materials where high hardness, ductility, and low density are required.

36 MATERIALS SCIENCE↗

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

36 MATERIALS SCIENCE↗

Insight into the Catalytic Mechanism of GH11 Xylanase: Computational Analysis of Substrate Distortion Based on a Neutron Structure

The reaction mechanism of biomass decomposition by xylanases remains the subject of debate. Here, to clarify the mechanism, we investigated the glycosylation step of GH11 xylanase, an enzyme that catalyzes the hydrolysis of lignocellulosic hemicellulose (xylan). Making use of a recent neutron crystal structure, which revealed the protonation states of relevant residues, we used ab initio quantum mechanics/molecular mechanics (QM/MM) calculations to determine the detailed reaction mechanism of the glycosylation step. In particular, our focus is on the controversial question of whether or not an oxocarbenium ion intermediate is formed on the reaction pathway. The calculations support the validity of a basic retaining mechanism within a double-displacement scheme. The estimated free energy barrier of this reaction is ~18 kcal/mol with QM/MM-CCSD(T)/6-31(+)G**//MP2/6-31+G**/AMBER calculations, and the rate-determining step of the glycosylation is scission of the glycosidic bond after proton transfer from the acidic Glu177. The estimated lifetime of the oxocarbenium ion intermediate (on the order of tens of ps) and the secondary kinetic isotope effect suggest that there is no accumulation of this intermediate on the reaction path, although the intermediate can be transiently formed. In the enzyme–substrate (ES) complex, the carbohydrate structure of the xylose residue at the -1 subsite has a rather distorted (skewed) geometry, and this xylose unit at the active site has an apparent half-chair conformation when the oxocarbenium ion intermediate is formed. The major catalytic role of the protein environment is to orient residues that take part in the initial proton transfer. Because of a fine alignment of catalytic residues, the enzyme can accelerate the glycosylation reaction without paying a reorganization energy penalty.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fractal dimensions of jammed packings with power-law particle size distributions in two and three dimensions

Static structure factors are computed for large-scale, mechanically stable, jammed packings of frictionless spheres (three dimensions) and disks (two dimensions) with broad, power-law size dispersity characterized by the exponent -β. Here, the static structure factor exhibits diverging power-law behavior for small wave numbers, allowing us to identify a structural fractal dimension d f . In three dimensions, d f ≈ 2.0 for 2.5≤ β ≤ 3.8, such that each of the structure factors can be collapsed onto a universal curve. In two dimensions, we instead find 1.0≲ d f ≲1.34 for 2.1≤ β ≤ 2.9. Furthermore, we show that the fractal behavior persists when rattler particles are removed, indicating that the long-wavelength structural properties of the packings are controlled by the large particle backbone conferring mechanical rigidity to the system. A numerical scheme for computing structure factors for triclinic unit cells is presented and employed to analyze the jammed packings.

97 MATHEMATICS AND COMPUTING↗

Integrated Computational Materials and Mechanical Modeling for Additive Manufacturing of Alloys with Graded Structure Used in Fossil Fuel Power Plants

Wire-arc additive manufacturing (WAAM) has demonstrated its unique capability of producing large-size alloy components with a significantly reduced fabrication time and enhanced geometry design freedom. In this project, the team has developed an ICME (Integrated Computational Materials Engineering) modeling framework, which supports the WAAM of the AUSC (Advanced Ultra-Supercritical) power plant components. The manufacturing design has been applied to Inconel 740H, steel P91, as well as the dissimilar alloy components between steel P91 and Inconel 740H. The ICME model framework is developed by considering two types of modeling. First, mechanistic modeling has been applied to control the printing quality and understand the sequence of the dissimilar printing of the wall structure. The following models have been included in the developed ICME framework: finite element thermal model, grain structure model, residual stress simulation, crystal plasticity model, CALPHAD-based precipitation kinetic model, phase stability prediction, thermal expansion predictive model, and heuristic creep model. Secondary, a physics-based machine learning model has also been developed based on the ICME model structure. The machine learning model development is based on the ICME model prediction with calibration of the experiments. In addition, the WAAM has been utilized as a high-throughput experimental tool rapidly generating a gradient of alloy composition to facilitate experimental database generation for process-structure-property relationships. Such a database directly supported the ICME-enhanced machine learning, which further assisted in intermediate composition block design between P91 and 740H. A high-throughput screening study of the oxidation resistance has been performed based on such high-throughput experimentation. Based on the computational design, several dissimilar alloy manufacturing with post-heat treatment have been performed with a comprehensive evaluation of mechanical performance, including hardness mapping, yield strength, creep resistance. In this project, the single component of P91 and 740H processed by WAAM after heat treatment designed by ICME has demonstrated higher performance in yield strength and creep resistance than the wrought materials. The P91 sample prepared by WAAM with ICME-designed heat treatment performs better than P92 in creep resistance. The designed graded alloy printing with intermediate block shows a promising performance that exceeds the traditional welding. Moreover, the current research indicates the high need for location-specific design analysis with uncertainty quantification, an important topic that deserves more dedicated research. The achievement of this project demonstrated the promising future of WAAM in structural alloy manufacturing for energy power plant development. Successful printing requires synergetic efforts made by manufacturing, mechanical, and materials sciences.

20 FOSSIL-FUELED POWER PLANTS↗

Mapping 3D grain and precipitate structure during in situ mechanical testing of open-cell metal foam using micro-computed tomography and high-energy X-ray diffraction microscopy

Open-cell metal foams are ultra-low-density cellular metals with complex hierarchical structures that span bulk, cell, ligament, and sub-ligament scales and give rise to desirable properties such as high strength-to-weight ratio and excellent energy absorption. Although literature suggests that intrinsic material structures at sub-ligament length scales (e.g., grains and precipitates) play an important role in mechanical behavior of open-cell metal foams, there are very few experimental measurements of such structures in three dimensions and for meaningful volumes of foam. This study seeks to map and track the three-dimensional (3D) grain and precipitate structures of an intact volume of open-cell aluminum foam by advancing microstructural characterization techniques that leverage X-ray micro-computed tomography (μCT) and far-field high-energy X-ray diffraction microscopy (FF-HEDM). A 6%-relative-density aluminum foam sample was mechanically tested in compression while μCT and FF-HEDM measurements were collected at interrupted loading states at beamline 1-ID of the Advanced Photon Source. Further, a new scanning strategy and reconstruction algorithm were established to enable characterization of a foam volume with diameter approximately four times wider than the nominal width of the X-ray beam. The result is a set of maps that detail both the 3D grain and precipitate structures throughout the foam volume at successive strain steps. A novel grain tracking procedure was developed to track individual grains within the foam volume by accounting for the large rigid-body motions that individual ligaments can undergo during mechanical loading. The ability to track grains and precipitate structures in three dimensions throughout large bulk deformation of ultra-low-density polycrystalline materials enables new possibilities for validating numerical models and investigating local failure mechanisms. Furthermore, the methods and procedures developed in this study could be applied to other ultra-low-density structures, such as additively manufactured lattices.

36 MATERIALS SCIENCE↗

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗

Four-Coordinate Fe N 2 and Imido Complexes Supported by a Hemilabile NNC Heteroscorpionate Ligand

Inspired by mechanistic proposals for N 2 reduction at the nitrogenase FeMo cofactor, we report herein a new, strongly σ-donating heteroscorpionate ligand featuring two weak-field pyrazoles and an alkyl donor. This ligand supports four-coordinate Fe(I)-N 2 , Fe(II)-Cl, and Fe(III)-imido complexes, which we have characterized using a variety of spectroscopic and computational methods. Structural and quantum mechanical analysis reveal the nature of the Fe–C bonds to be essentially invariant between the complexes, with conversion between the (formally) low-valent Fe-N 2 and high-valent Fe-imido complexes mediated by pyrazole hemilability. This presents a useful strategy for substrate reduction at such low-coordinate centers and suggests a mechanism by which FeMoco might accommodate the binding of both π-acidic and π-basic nitrogenous substrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strategies for simultaneous strengthening and toughening via nanoscopic intracrystalline defects in a biogenic ceramic

While many organisms synthesize robust skeletal composites consisting of spatially discrete organic and mineral (ceramic) phases, the intrinsic mechanical properties of the mineral phases are poorly understood. Using the shell of the marine bivalve Atrina rigida as a model system, and through a combination of multiscale structural and mechanical characterization in conjunction with theoretical and computational modeling, we uncover the underlying mechanical roles of a ubiquitous structural motif in biogenic calcite, their nanoscopic intracrystalline defects. These nanoscopic defects not only suppress the soft yielding of pure calcite through the classical precipitation strengthening mechanism, but also enhance energy dissipation through controlled nano- and micro-fracture, where the defects’ size, geometry, orientation, and distribution facilitate and guide crack initialization and propagation. These nano- and micro-scale cracks are further confined by larger scale intercrystalline organic interfaces, enabling further improved damage tolerance.

36 MATERIALS SCIENCE↗