Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response 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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Tunable nano-wrinkling of chiral surfaces: Structure and diffraction optics

Periodic surface nano-wrinkling is found throughout biological liquid crystalline materials, such as collagen films, spider silk gland ducts, exoskeleton of beetles, and flower petals. These surface ultrastructures are responsible for structural colors observed in some beetles and plants that can dynamically respond to external conditions, such as humidity and temperature. In this paper, the formation of the surface undulations is investigated through the interaction of anisotropic interfacial tension, swelling through hydration, and capillarity at free surfaces. Focusing on the cellulosic cholesteric liquid crystal (CCLC) material model, the generalized shape equation for anisotropic interfaces using the Cahn-Hoffman capillarity vector and the Rapini-Papoular anchoring energy are applied to analyze periodic nano-wrinkling in plant-based plywood free surfaces with water-induced cholesteric pitch gradients. Scaling is used to derive the explicit relations between the undulations’ amplitude expressed as a function of the anchoring strength and the spatially varying pitch. The optical responses of the periodic nano-structured surfaces are studied through finite difference time domain simulations indicating that CCLC surfaces with spatially varying pitch reflect light in a wavelength higher than that of a CCLC’s surface with constant pitch. This structural color change is controlled by the pitch gradient through hydration. All these findings provide a foundation to understand structural color phenomena in nature and for the design of optical sensor devices.

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

Light-induced dilation in nanosheets of charge-transfer complexes

We observe a sizable photostrictive effect of 5.7% with fast, submillisecond response times, arising from a light-induced lattice dilation of a molecular nanosheet, composed of the molecular charge-transfer compound dibenzotetrathiafulvalene (DBTTF) and C 60 . An interfacial self-assembly approach is introduced for the thickness-controlled growth of the thin films. From photoabsorption measurements, molecular simulations, and electronic structure calculations, we suggest that photostriction within these films arises from a transformation in the molecular structure of constituent molecules upon photoinduced charge transfer, as well as the accommodation of free charge carriers within the material. Additionally, we find that the photostrictive properties of the nanosheets are thickness-dependent, a phenomenon that we suggest arises from surface-induced conformational disorder in the molecular components of the film. Moreover, because of the molecular structure in the films, which results largely from interactions between the constituent π-systems and the sulfur atoms of DBTTF, the optoelectronic properties are found to be anisotropic. This work enables the fabrication of 2D molecular charge-transfer nanosheets with tunable thicknesses and properties, suitable for a wide range of applications in flexible electronic technologies.

2D nanosheets↗
DOE OSTI2015

Modeling the Mesoscale Transport of Lithium-Magnetite Electrodes Using Insight from Discharge and Voltage Recovery Experiments

A multi-scale mathematical model, which accounts for mass transport on the crystal and agglomerate length-scales, is used to investigate the electrochemical performance of lithium-magnetite electrochemical cells. Experimental discharge and voltage recovery data are compared to three sets of simulations, which incorporate crystal-only, agglomerate-only, or multi-scale transport effects. Mass transport diffusion coefficients are determined by fitting the simulated voltage recovery times to experimental data. In addition, a further extension of the multi-scale model is proposed which accounts for the impact of agglomerate size distributions on electrochemical performance. The results of the study indicate that, depending on the crystal size, the low utilization of the active material is caused by transport limitations on the agglomerate and/or crystal length-scales. For electrodes composed of small crystals (6 and 8 nm diameters), it is concluded that the transport limitations in the agglomerate are primarily responsible for the long voltage recovery times and low utilization of the active mass. In the electrodes composed of large crystals (32 nm diameter), the slow voltage recovery is attributed to transport limitations on both the agglomerate and crystal length-scales.

25 ENERGY STORAGE↗
NASA NTRS1985

Occurrence of giant impacts during the growth of the terrestrial planets

Three dimensional Monte Carlo simulations of the accumulation of the terrestrial planets in the absence of gas drag produced results that are in general agreement with the number and distribution of the present planets. The accumulation process appears to be characterized by impact of bodies as large as three times the mass of Mars at velocities of about 9 kilometers per second. These giant impacts on earth may have supplied the material and angular momentum that formed the moon, should have heated earth to the melting point, and may have been responsible for the differences in the content of inert gases of the atmospheres of earth and Venus.

Wetherill, G. W.↗
NASA NTRS1989

Computational procedures for postbuckling of composite shells

A recently developed finite-element capability for general nonlinear shell analysis, featuring the use of three-dimensional constitutive equations within an efficient resultant-oriented framework, is employed to simulate the postbuckling response of an axially compressed composite cylindrical panel with a circular cutout. The problem is a generic example of modern composite aircraft components for which postbuckling strength (i.e., fail-safety) is desired in the presence of local discontinuities such as holes and cracked stiffeners. While the computational software does a reasonable job of predicting both the buckling load and the qualitative aspects of postbuckling (compared both with experiment and another code) there are some discrepancies due to: (1) uncertainties in the nominal layer material properties, (2) structural sensitivity to initial imperfections, and (3) the neglect of dynamic and local material delamination effects in the numerical model. Corresponding refinements are suggested for the realistic continuation of this type of analysis.

Stanley, G. M.↗
DOE OSTI2023

Microstructure characterization and elastic-plastic self-consistent simulation studies of anisotropic deformation of β -tin

Uniaxial compression tests of 99.9% pure polycrystalline $β$-tin (Sn) were conducted at various strain rates (10 -3 /s, 10 -1 /s and 1/s) and temperatures (294 K and 193 K) to understand the effects of these variables on the stress–strain response, microstructure evolution, and stress relaxation behavior. Multiple Sn specimens were subjected to complex compressive loading/unloading/reloading paths at different strain rates. Specimens initially compressed at higher strain rates showed strain rate-dependent texture evolution and more pronounced relaxation upon reloading, as compared to those pre-strained at lower rates at room temperature (RT). Compression tests conducted at low temperature (193 K) revealed increased strength, similar to the enhanced strength observed when the strain rate was increased at RT. Neutron diffraction was employed to characterize the initial and final bulk textures of the RT specimens. Electron backscatter diffraction was utilized to examine the crystallographic grain orientation and morphology, thereby identifying the signatures of dislocation-mediated deformation, grain refinement, recrystallization, and twinning behaviors. Elastic-plastic self-consistent simulations were performed to investigate the deformation modes responsible for the strain rate-dependent macroscopic stress–strain response and texture evolution. Effects of crystallographic orientation on stress relaxation behavior was also examined. In conclusion, the model predictions are in reasonable agreement with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Dependence of the rate of LiF ion pairing on the description of molecular interaction

We present an analysis of the dynamics of ion-pairing of Lithium Fluoride (LiF) in aqueous solvent using both detailed molecular simulation as well as reduced models within a Gener- alized Langevin Equation (GLE) framework. We explored the sensitivity of the ion-pairing phenomena to the details of descriptions of molecular interaction, comparing two empirical potentials to explicit quantum based density functional theory. We find quantitative differences in the potentials of mean force for ion-pairing as well as time dependent frictions that lead to variations in the rate constant and reactive flux correlation functions. These details reflect differences in solvent response to ion-pairing between different representations of molecular interaction and influence anharmonicity of the dynamic response. We find that the short time anharmonic response is recovered with a GLE parameterization. Recovery of the details of long time response may require extensions to the reduced model. We show that the utility of using a reduced model leads to a straight forward application of variational transition state the- ory concepts to the condensed phase system. The significance of this is reflected in the analysis of committor distributions and the variation of planar hypersurfaces, leading to an improved understanding of factors that determine the rate of LiF ion-pairing. CJM and GKS are supported by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest Na- tional Laboratory (PNNL) is operated for the Department of Energy by Battelle. MDB is grateful for the support of Laboratory Directed Research and Development funding under the auspices of PNNL’s Laboratory Initiative Materials Synthesis and Simulation across Scales (MS3). Additional computing resources were generously allocated by PNNL’s Institutional Computing program. EP acknowledges support from PNNL’s Alternate Sponsored Fellowship program and IMPRS Dres- den.Support to P.J. from the Czech Science Foundation (grant P208/12/G016) and the Academy of Sciences (Praemium Academie award) is gratefully acknowledged.

Pluharova, Eva↗
DOE OSTI2023

Intermetallics Based on Sodium Chalcogenides Promote Stable Electrodeposition–Electrodissolution of Sodium Metal Anodes

Sodiophilic micro-composite films of sodium-chalcogenide intermetallics (Na 2 Te and Na 2 S) and Cu particles are fabricated onto commercial copper foam current collectors (Na 2 Te@CF and Na 2 S@CF). For the first time a controllable capacity thermal infusion process is demonstrated. Enhanced wetting by the metal electrodeposition leads to state-of-the-art electrochemical performance. For example, Na 2 Te@CF-based half-cells demonstrate stable cycling at 6 mA cm -2 and 6 mAh cm -2 , corresponding to 54 µm of Na electrodeposited/electrodissolved by geometric area. Sodium metal batteries with Na 3 V 2 (PO 4 ) 3 cathodes are stable at 30C (7 mA cm -2 ) and for 10 000 cycles at 5C and 10C. Cross-sectional cryogenic focused ion beam (cryo-FIB) microscopy details deposited and remnant dissolved microstructures. Sodium metal electrodeposition onto Na 2 Te@CF is dense, smooth, and free of dendrites or pores. On unmodified copper foam, sodium grows in a filament-like manner, not requiring cycling to achieve this geometry. Substrate–metal interaction critically affects the metal–electrolyte interface, namely the thickness and morphology of the solid electrolyte interphase. Density functional theory and mesoscale simulations provide insight into support-adatom energetics, nucleation response, and early-stage morphological evolution. On Na 2 Te sodium atomic dispersion is thermodynamically more stable than isolated clusters, leading to conformal adatom coverage of the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Nanoantennas for nanowire photovoltaics

We consider the use of plasmonic nanoantenna elements, hemispherical and cylindrical, for application in semiconductor nanowire (NW) vertical arrays. Using Mie theory and a finite element method, scattering and absorption efficiencies are obtained for the desired enhancement of interaction with light in the NWs. We find an optimal mixture of nanoantennae for efficient scattering of solar spectrum in the NW array. Spectral radiation patterns of scattered light are computed, and, for representing the total response of the nanoantenna-equipped NWs to the solar AM1.5G spectrum, the weighted average of scattering patterns for unpolarized normal incidence is obtained showing an advantageous overall directivity toward the NWs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Multi-scale modeling of inter-granular fracture in UO2

A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Effects of geometry and composition of soft polymer films embedded with nanoparticles on rates for optothermal heat dissipation

Embedding soft matter with nanoparticles (NPs) can provide electromagnetic tunability at sub-micron scales for growing number of uses in healthcare, sustainable energy, and chemical processing. But use of NP-embedded soft material in temperature-sensitive applications has been constrained by difficulty in validating prediction of rates for energy dissipation across thermally insulating to conducting behavior. This work improved embedment of monodisperse NP to stably decrease inter-NP spacings in polydimethylsiloxane (PDMS) to nano-scales. Lumped-parameter and finite element analyses were refined to apportion effects of structure and composition of NP-embedded soft polymer to rates for conductive, convective, and radiative heat dissipation. These advances allowed rational selection of PDMS size and NP composition to optimize measured rates of internal (conductive) and external (convective and radiative) heat dissipation. Stably reducing distance between monodisperse NP to nano-scale intervals increased overall heat dissipation rate by up to 29%. Refined fabrication of NP-embedded polymer enabled tunability of dynamic thermal response, the ratio of internal to external dissipation rate, by a factor of 3.1 to achieve a value of 0.091, the largest reported to date. Heat dissipation rates simulated a priori were consistent with 130 μm resolution thermal images across 2- to 15-fold changes in NP-PDMS geometry and composition. Nusselt number was observed to increase with the fourth root of Rayleigh number across thermally insulative and conductive regimes, further validating the approach. These developments support model-informed design of soft media embedded with nano-scale spaced NP to optimize heat dissipation rates for evolving temperature-sensitive diagnostic and therapeutic modalities as well as emerging uses in flexible bioelectronics, cell and tissue culture, and solar-thermal heating.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Photoluminescence Intensity Fluctuations and Temperature-Dependent Decay Dynamics of Individual Carbon Nanotube sp 3 Defects

Recent demonstration of room temperature, telecommunication wavelength single photon generation by sp 3 defects of single wall carbon nanotubes established these defects as a new class of quantum materials. However, their practical utilization in development of quantum light sources calls for a significant improvement in their imperfect quantum yield (QY~10–30%). PL intensity fluctuations observed with some defects also need to be eliminated. Aiming toward attaining fundamental understanding necessary for addressing these critical issues, we investigate PL intensity fluctuation and PL decay dynamics of aryl sp 3 defects of (6,5), (7,5), and (10,3) single wall carbon nanotubes (SWCNTs) at temperatures ranging from 300 to 4 K. By correlating defect-state PL intensity fluctuations with change (or lack of change) in PL decay dynamics, we identified random variations in the trapping efficiency of E 11 band-edge excitons (likely resulting from the existence of a fluctuating potential barrier in the vicinity of the defect) as the mechanism mainly responsible for the defect PL intensity fluctuations. Furthermore, by analyzing the temperature dependence of PL intensity and decay dynamics of individual defects based on a kinetic model involving the trapping and detrapping of excitons by optically allowed and forbidden (bright and dark) defect states, we estimate the height of the potential barrier to be in the 3–22 meV range. Our analysis also provides further confirmation of recent DFT simulation results that the emissive sp 3 defect state is accompanied by an energetically higher-lying optically forbidden (dark) exciton state.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides

Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Creep fronts and complexity in laboratory earthquake sequences illuminate delayed earthquake triggering

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.

58 GEOSCIENCES↗
DOE OSTI2022

Requirements for Soldering Fluxes Research using the B-53 Test Board

IPC J-STD-004B standard prescribes general requirements for the classification and testing of soldering flux for high qualify interconnections. This standard defines the classification of soldering materials through specifications of test methods and inspection criteria. The materials include liquid flux, paste flux, solderpaste flux, solder preform flux, and flux-cored solder. This research will use the proposed IPC-53 Surface Insulation Resistance (SIR) test patterns by means of an open comb (2D) and closed comb (3D to simulate a component over the comb pattern). The 2D open comb has uniformity of conductor spacing, sheet resistance, and flux outgassing. The 3D-closed comb simulates the effect of leadless or bottom-terminated components, which have non-uniform sheet resistance and flux outgassing. The response variables will include SIR testing and visual imaging. The objective is to investigate IPC test method improvements for characterizing soldering fluxes when using leadless components with narrow pad-to-pad spacing.

Diamond, Louis↗
DOE OSTI2014

Influence of crystallographic orientation on the response of copper crystallites to nanoindentation

Molecular dynamics simulation was performed to study the features of nucleation and development of plastic deformation in copper crystallites in nanoindentation with different crystallographic orientations of their loaded surface: (011), (001), and (111). Atomic interaction was described by a potential constructed in terms of the embedded atom method. It is shown that behavior of the crystallite reaction force correlates well with a change in the fraction of atoms involved in local structural rearrangements. The generation of local structural changes decreases the slope of the crystallite reaction force curve or results in an extremum due to internal stress relaxation. Analysis of structural changes in the material being indented demonstrates that the orientation of its loaded surface greatly affects the features of nucleation and development of plastic deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
NASA NTRS1998

Effect of Apex Flap Deflection on Vertical Tail Buffeting

A computational study of the effect of vortex breakdown location on vertical tail buffeting is conducted. The position of the breakdown is modified by employing an apex flap deflected by an experimentally determined optimal angle. The delayed breakdown flow and buffeting response is then compared to the nominal undeflected case. This multidisciplinary problem is solved sequentially for the fluid flow, the elastic tail deformations and the grid displacements. The fluid flow is simulated by time accurately solving the unsteady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, upwind, flux-difference splitting finite volume scheme. The elastic vibrations of the tails are modeled by uncoupled bending and torsion beam equations. These equations are solved accurately in time using the Galerkin method and a five-stage Runge-Kutta-Verner scheme. The grid for the fluid dynamics calculations is continuously deformed using interpolation functions to disperse the displacements smoothly throughout the computational domain. An angle-of-attack of 35 deg.is chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on tile vertical tails. The dimensions and material properties of the vertical tails are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. Results are presented for a baseline uncontrolled buffeting case and a delayed breakdown case in which the apex flap has been deflected 15 deg. The flap was found to be very effective in delaying the breakdown, increasing the location from 50%c to 94%c, which resulted in a 6% increase in lift coefficient and pitching moment. However, the integrated buffet loads and tip responses were roughly equivalent for the two cases.

Massey, Steven J.↗

Refine your search to explore more results.