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,279 records · Page 71

Secondary-Source Core Reload Modeling with VERA

The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary-source assemblies, use of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous-energy Monte Carlo solver Shift for final detector response evaluation using the CADIS methodology for variance reduction. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data for the first few fuel loading steps and were found to demonstrate excellent agreement with the measured data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

In situ characterization of laser-generated melt pools using synchronized ultrasound and high-speed X-ray imaging

Metal additive manufacturing is a fabrication method that forms a part by fusing layers of powder to one another. An energy source, such as a laser, is commonly used to heat the metal powder sufficiently to cause a molten pool to form, which is known as the melt pool. The melt pool can exist in the conduction or the keyhole mode where the material begins to rapidly evaporate. The interaction between the laser and the material is physically complex and difficult to predict or measure. Here, high-speed X-ray imaging was combined with immersion ultrasound to obtain synchronized measurements of stationary laser-generated melt pools. Furthermore, two-dimensional and three-dimensional finite-element simulations were conducted to help explain the ultrasonic response in the experiments. In particular, the time-of-flight and amplitude in pulse-echo configuration were observed to have a linear relationship to the depth of the melt pool. These results are promising for the use of ultrasound to characterize the melt pool behavior and for finite-element simulations to aid in interpretation.

42 ENGINEERING↗

Reducing Uncertainty in High-Resolution Sea Ice Models

Arctic sea ice is an important component of the global climate system, reflecting a significant amount of solar radiation, insulating the ocean from the atmosphere and influencing ocean circulation by modifying the salinity of the upper ocean. The thickness and extent of Arctic sea ice have shown a significant decline in recent decades with implications for global climate as well as regional geopolitics. Increasing interest in exploration as well as climate feedback effects make predictive mathematical modeling of sea ice a task of tremendous practical import. Satellite data obtained over the last few decades have provided a wealth of information on sea ice motion and deformation. The data clearly show that ice deformation is focused along narrow linear features and this type of deformation is not well-represented in existing models. To improve sea ice dynamics we have incorporated an anisotropic rheology into the Los Alamos National Laboratory global sea ice model, CICE. Sensitivity analyses were performed using the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) to determine the impact of material parameters on sea ice response functions. Two material strength parameters that exhibited the most significant impact on responses were further analyzed to evaluate their influence on quantitative comparisons between model output and data. The sensitivity analysis along with ten year model runs indicate that while the anisotropic rheology provides some benefit in velocity predictions, additional improvements are required to make this material model a viable alternative for global sea ice simulations.

Peterson, Kara J.↗

Buckling and postbuckling analyses of laminated anisotropic structures

A review is given of recent advances in two aspects of the numerical simulation of the buckling and postbuckling responses of composite structures. The first aspect is exploiting non-traditional symmetries exhibited by composite structures; and strategies for reducing the size of the model and the cost of the buckling and postbuckling analyses in the presence of symmetry-breaking conditions (e.g., asymmetry of the material, geometry, and/or loading). The second aspect pertains to the prediction of onset of local delamination in the postbuckling range and accurate determination of transverse shear stresses in the structure. The accuracy and effectiveness of the strategies developed are demonstrated by means of a numerical example.

Noor, Ahmed K.↗

Advanced Porous Coating for Low-Density Ceramic Insulation Materials

The need for improved coatings on low-density reusable surface insulation (RSI) materials used on the space shuttle has stimulated research into developing tougher coatings. The processing of a new porous composite "coating" for RST called toughened unipiece fibrous insulation Is discussed. Characteristics including performance in a simulated high-speed atmospheric entry, morphological structure before and after this exposure, resistance to Impact, and thermal response to a typical heat pulse are described. It is shown that this coating has improved impact resistance while maintaining optical and thermal properties comparable to the previously available reaction-cured glass coating.

Leiser, Daniel B.↗

Modeling Dynamic Crush Behavior of Carbon Fiber Reinforced Polymer Composite Structures Using MAT213

Modeling crushing of carbon fiber reinforced polymer (CFRP) composites is challenging, and current simulation methodologies involve tuning of non-physical parameters. MAT213, a next-generation material model, has advanced functionality to better simulate dynamic impact loading. The objective of the present investigation is to evaluate the potential for using MAT213 to simulate dynamic crushing of CFRPs. Two sets of simulations were performed: one for calibration based on a set of coupon-level experiments and another for prediction of the response of structural elements. Simulations involving dynamic crushing of flat specimens were iteratively run to calibrate model parameters. The calibration demonstrated that MAT213 could produce a simulated force-displacement response within experimental scatter. The simulated failure morphology was also comparable to the experiments. After successful calibration, predictive simulations of dynamic crushing of C-channel shaped specimens were completed using a simulated crash sled test rig and two pairs of impactor mass/velocity conditions. The simulated force-displacement curve in the crash sled simulations for the lower-velocity condition fell within the experimental scatter, but the stable crush force was underpredicted by 27% in the higher-velocity simulations. Better correlation in the lower-velocity test condition likely results from the calibration condition being a similar velocity to the lower-velocity crash sled condition.

Dynamic Crush↗

Performance Study of an aSi Flat Panel Detector for Fast Neutron Imaging of Nuclear Waste

Radioactive waste must be characterized to check its conformance for intermediate storage and final disposal according to national regulations. For the determination of radio-toxic and chemo-toxic contents of radioactive waste packages non-destructive analytical techniques are preferentially used. Fast neutron imaging is a promising technique to assay large and dense items providing, in complementarity to photon imaging, additional information on the presence of structures in radioactive waste packages. Therefore the feasibility of a compact Neutron Imaging System for Radioactive waste Analysis (NISRA) using 14 MeV neutrons is studied in a cooperation framework of Forschungszentrum Juelich GmbH, RWTH Aachen University and Siemens AG. However due to the low neutron emission of neutron generators in comparison to research reactors the challenging task resides in the development of an imaging detector with a high efficiency, a low sensitivity to gamma radiation and a resolution sufficient for the purpose. The setup is composed of a commercial D-T neutron generator (Genie16GT, Sodern) with a surrounding shielding made of polyethylene, which acts as a collimator and an amorphous silicon flat panel detector (aSi, 40 x 40 cm{sup 2}, XRD-1642, Perkin Elmer). Neutron detection is achieved using a general propose plastic scintillator (EJ-260, Eljen Technology) linked to the detector. The thermal noise of the photodiodes is reduced by employing an entrance window made of aluminium. Optimal gain and integration time for data acquisition are set by measuring the response of the detector to the radiation of a 500 MBq {sup 241}Am-source. Detector performance was studied by recording neutron radiography images of materials with various, but well known, chemical compositions, densities and dimensions (Al, C, Fe, Pb, W, concrete, polyethylene, 5 x 8 x 10 cm{sup 3}). To simulate gamma-ray emitting waste radiographs in presence of a gamma-ray sources ({sup 60}Co, {sup 137}Cs, {sup 241}Am) were performed. A homemade algorithm was developed to determine a value which is related to the neutron absorption of the sample with the analysis of the raw detector data. The detector was placed 42 cm away from the neutron source. Distance between detector and the samples was 0.5 cm. At the sample position the fast neutron flux was estimated to 9x10{sup 3} n cm{sup -2} s{sup -1} for a neutron emission of 10{sup 8} n s{sup -1}. The acquisition time was 15 minutes. First neutron radiographs were successfully recorded despite the low detector efficiency and low neutron emission. Analysis of the data shows a correlation between the measured signal and determined neutron absorption. Thus discrimination between different materials of same thicknesses may be achieved. The measurements and results will be presented and discussed in details.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Human Factors Guidelines for UAS in the National Airspace System

The ground control stations (GCS) of some UAS have been characterized by less-than-adequate human-system interfaces. In some cases this may reflect a failure to apply an existing regulation or human factors standard. In other cases, the problem may indicate a lack of suitable guidance material. NASA is leading a community effort to develop recommendations for human factors guidelines for GCS to support routine beyond-line-of-sight UAS operations in the national airspace system (NAS). In contrast to regulations, guidelines are not mandatory requirements. However, by encapsulating solutions to identified problems or areas of risk, guidelines can provide assistance to system developers, users and regulatory agencies. To be effective, guidelines must be relevant to a wide range of systems, must not be overly prescriptive, and must not impose premature standardization on evolving technologies. By assuming that a pilot will be responsible for each UAS operating in the NAS, and that the aircraft will be required to operate in a manner comparable to conventionally piloted aircraft, it is possible to identify a generic set of pilot tasks and the information, control and communication requirements needed to support these tasks. Areas where guidelines will be useful can then be identified, utilizing information from simulations, operational experience and the human factors literature. In developing guidelines, we recognize that existing regulatory and guidance material will, at times, provide adequate coverage of an area. In other cases suitable guidelines may be found in existing military or industry human factors standards. In cases where appropriate existing standards cannot be identified, original guidelines will be proposed.

Guidelines↗

Symmetry Breaking in Tetrahedral Chiral Plasmonic Nanoparticle Assemblies

Self-assembled plasmonic structures combine the specificity and tunability of chemical synthesis with collective plasmonic properties. In this work we systematically explore the effects of symmetry breaking on the chiroptical response of an assembly of plasmonic nanoparticles using simulation. The design is based on a tetrahedral nanoparticle frame with two different types of nanoparticles, where chirality is induced by targeted stimuli that change the distance along one edge of the assembly. Here, we show that the intensity, spectral position, and handedness of the CD response are tunable with small structural changes, making it usable as a nanoscale plasmonic ruler. We then build upon this initial design to show that the symmetry breaking principle may also be used to design a chiral pyramid using a mixture of different nanoparticle materials, which affords tunability over a broad spectral range, and retrieves nanoscale conformational changes over a range of length scales.

36 MATERIALS SCIENCE↗

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗

Derivatization and diffusive motion of molecular fullerenes: Ab initio and atomistic simulations

Using first principles density functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of derivatization on the electronic and transport properties of C{sub 60} fullerene. As a typical example, we consider [6,6]-phenyl-C{sub 61}-butyric acid methyl ester (PCBM), which forms one of the most efficient organic photovoltaic materials in combination with electron donating polymers. Extra peaks are observed in the density of states (DOS) due to the formation of new electronic states localized at/near the attached molecule. Despite such peculiar behavior in the DOS of an isolated molecule, derivatization does not have a pronounced effect on the electronic transport properties of the fullerene molecular junctions. Both C{sub 60} and PCBM show the same response to finite voltage biasing with new features in the transmission spectrum due to voltage induced delocalization of some electronic states. We also study the diffusive motion of molecular fullerenes in ethanol solvent and inside poly(3-hexylthiophene) lamella using reactive molecular dynamics simulations. We found that the mobility of the fullerene reduces considerably due to derivatization; the diffusion coefficient of C{sub 60} is an order of magnitude larger than the one for PCBM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The effects of Si-doped prelayers on the optical properties of InGaN/GaN single quantum well structures

In this paper, we report on the effects of including Si-doped (In)GaN prelayers on the low temperature optical properties of a blue-light emitting InGaN/GaN single quantum well. We observed a large blue shift of the photoluminescence peak emission energy and significant increases in the radiative recombination rate for the quantum well structures that incorporated Si-doped prelayers. Simulations of the variation of the conduction and valence band energies show that a strong modification of the band profile occurs for the quantum wells on Si-doped prelayers due to an increase in strength of the surface polarization field. The enhanced surface polarization field opposes the built-in field across the quantum well and thus reduces this built-in electric field. This reduction of the electric field across the quantum well reduces the Quantum Confined Stark Effect and is responsible for the observed blue shift and the change in the recombination dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗

TU-F-18A-07: To Explore the More Realistic Energy Responses of the In-Depth Photon Counting Detectors

Purpose: We study the effect of the secondary photon events on modeling the energy response functions of the In-depth photon counting X-ray detectors (PCXD) and the potential impact of the spectral distortion on material decompositions. Methods: Square-shape wafers of three potential PCXD materials were constructed (5-by-20-by-30 mm{sup 3} for Si, 4-by-20-by-5 mm{sup 3} for GaAs and 4-by-20-by-3 mm{sup 3} for CdTe), with pixel size of 5-by-4 mm{sup 2} for Si and 5-by-5 mm{sup 2} for GaAs and CdTe. The depth direction (z-direction) was segmented into 5 layers with exponentially increasing thicknesses of each layer. X-rays from 10keV to 120keV with 20000 photons per keV bin was simulated to characterize the energy response function of each PCXD using Geant4. Secondary photons events were recorded and we omitted the photons exiting the detector. The Energy Response Functions (ERFs) from the Monte Carlo (MC) simulations were compared with those from a semi-ideal model developed earlier. Results: For Si, detection of secondary events in the center detector were minimal due to the long aspect ratio of the detector, which results in the agreement between the theoretical prediction and the MC simulation with and without the secondary photons. For CdTe, the secondary photons captured by the center pixel were important, leading to obvious disagreement between the analytical and the simulated ERF. After correction for secondary events, the two curves were more similar except for the escape peaks which are not correctly portrayed by the semi-ideal model. For GaAs, the behavior is in between Si and CdTe. Conclusion: Given the complexity of the In-Depth PCXD's geometry, the uniform semi-ideal model does not fully characterize the ERF at each layer. Therefore, more realistic models need to be explored for better modeling of the spectral distortion.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development and Operation of a Material Identification and Discrimination Imaging Spectroradiometer

Many imaging applications require quantitative determination of a scene's spectral radiance. This paper describes a new system capable of real-time spectroradiometric imagery. Operating at a full-spectrum update rate of 30Hz, this imager is capable of collecting a 30 point spectrum from each of three imaging heads: the first operates from 400 nm to 950 nm, with a 2% bandwidth; the second operates from 1.5 micro-m to 5.5 micro-m with a 1.5% bandwidth; the third operates from 5 micro-m to 12 micro-m, also at a 1.5% bandwidth. Standard image format is 256 x 256, with 512 x 512 possible in the VIS/NIR head. Spectra of up to 256 points are available at proportionately lower frame rates. In order to make such a tremendous amount of data more manageable, internal processing electronics perform four important operations on the spectral imagery data in real-time. First, all data in the spatial/spectral cube of data is spectro-radiometrically calibrated as it is collected. Second, to allow the imager to simulate sensors with arbitrary spectral response, any set of three spectral response functions may be loaded into the imager including delta functions to allow single wavelength viewing; the instrument then evaluates the integral of the product of the scene spectral radiances and the response function. Third, more powerful exploitation of the gathered spectral radiances can be effected by application of various spectral-matched filtering algorithms to identify pixels whose relative spectral radiance distribution matches a sought-after spectral radiance distribution, allowing materials-based identification and discrimination. Fourth, the instrument allows determination of spectral reflectance, surface temperature, and spectral emissivity, also in real-time. The spectral imaging technique used in the instrument allows tailoring of the frame rate and/or the spectral bandwidth to suit the scene radiance levels, i.e., frame rate can be reduced, or bandwidth increased to improve SNR when viewing low radiance scenes. The unique challenges of design and calibration are described. Pixel readout rates of 160 MHz, for full frame readout rates of 1000 Hz (512 x 512 image) present the first challenge; processing rates of nearly 600 million integer operations per second for sensor emulation, or over 2 billion per second for matched filtering, present the second. Spatial and spectral calibration of 66,536 pixels (262,144 for the 512 x 512 version) and up to 1,000 spectral positions mandate novel decoupling methods to keep the required calibration memory to a reasonable size. Large radiometric dynamic range also requires care to maintain precision operation with minimum memory size.

Dombrowski, Mark↗

Electron transport and electron energy distributions within the wurtzite and zinc-blende phases of indium nitride: Response to the application of a constant and uniform electric field

Within the framework of an ensemble semi-classical three-valley Monte Carlo electron transport simulation approach, we critically contrast the nature of the electron transport that occurs within the wurtzite and zinc-blende phases of indium nitride in response to the application of a constant and uniform electric field. We use the electron energy distribution and its relationship with the electron transport characteristics in order to pursue this analysis. For the case of zinc-blende indium nitride, only a peak corresponding to the electrons within the lowest energy conduction band valley is observed, this peak being seen to broaden and shift to higher energies in response to increases in the applied electric field strength, negligible amounts of upper energy conduction band valley occupancy being observed. In contrast, for the case of wurtzite indium nitride, in addition to the aforementioned lowest energy conduction band valley peak in the electron energy distribution, and its broadening and shifting to higher energies in response to increases in the applied electric field strength, beyond a certain critical electric field strength, 30 kV/cm for the case of this particular material, upper energy conduction band valley occupancy is observed, this occupancy being further enhanced in response to further increases in the applied electric field strength. Reasons for these results are provided. The potential for device consequences is then commented upon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗