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 1,207 records · Page 67

Performance Evaluation of the Gravity Probe B Design

This report documents the simulation of the Lockheed Martin designed Gravity Probe B (GPB) spacecraft developed tool by bd Systems Inc using the TREETOPS simulation. This study quantifies the effects of flexibility and liquid helium slosh on GPB spacecraft control performance. The TREETOPS simulation tool permits the simulation of flexible structures given that a flexible body model of the structure is available. For purposes of this study, a flexible model of the GPB spacecraft was obtained from Lockheed Martin. To model the liquid helium slosh effects, computational fluid dynamics (CFD) results' were obtained, and used to develop a dynamic model of the slosh effects. The flexible body and slosh effects were incorporated separately into the TREETOPS simulation, which places the vehicle in a 650 km circular polar orbit and subjects the spacecraft to realistic environmental disturbances and sensor error quantities. In all of the analysis conducted in this study the spacecraft is pointed at an inertially fixed guide star (GS) and is rotating at a constant rate about this line of sight.

Francis, Ronnie↗

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗

3D Bacterial flagella as both synthetic biotemplates and ultrathin spacers for enhanced inter-particle coupling and solar energy harvesting

Linear light-absorbing nanomaterials are ideal for film-based solar harvesting applications as they form porous structures that can maximize the absorption and minimize the reflection of the solar light. Conventional 1D nanochains of plasmonic nanoparticle assemblies can achieve significantly broadened optical absorption through surface plasmon coupling, but their optical bands are still not broad enough to absorb through the solar spectrum and thus are not efficient solar absorbers. Here we discovered first by simulation that 3D structured nanochains of plasmonic nanoparticles presented a remarkably increased optical broadening effect and much longer redshift of the optical peaks due to the enhanced inter-particle coupling effect. Then we fabricated 3D nanochains by assembling gold nanoparticles (AuNPs) around 14 nm ultrathin bionanofibers, the bacterial flagella. Additionally, the ultrathin biotemplates enabled the 3D arrangement of 50 nm AuNPs along the nanofiber with a very small inter-particle gap, allowing the strong coupling of surface plasmons in a 3D manner. Consistent with the theoretical prediction, the 3D nanochains, when assembled into films, could effectively convert nearly the full spectrum of solar energy into heat, which was further efficiently converted into electricity through a thermoelectric generation unit. Our work represents a nanobiomaterial approach to highly efficient solar thermal power generation.

36 MATERIALS SCIENCE↗

Large-eddy simulation of flow around an airfoil on a structured mesh

The diversity of flow characteristics encountered in a flow over an airfoil near maximum lift taxes the presently available statistical turbulence models. This work describes our first attempt to apply the technique of large-eddy simulation to a flow of aeronautical interest. The challenge for this simulation comes from the high Reynolds number of the flow as well as the variety of flow regimes encountered, including a thin laminar boundary layer at the nose, transition, boundary layer growth under adverse pressure gradient, incipient separation near the trailing edge, and merging of two shear layers at the trailing edge. The flow configuration chosen is a NACA 4412 airfoil near maximum lift. The corresponding angle of attack was determined independently by Wadcock (1987) and Hastings & Williams (1984, 1987) to be close to 12 deg. The simulation matches the chord Reynolds number U(sub infinity)c/v = 1.64 x 10(exp 6) of Wadcock's experiment.

Kaltenbach, Hans-Jakob↗

Systematic Validation of Ensemble Cloud-Process Simulations using Polarimetric Radar Observations and Simulator over the NASA Wallops Flight Facility

The BiLateral Operational Storm-Scale Observation and Modeling (BLOSSOM) project was initiated to establish a long-term supersite to improve understanding of cloud physical states and processes as well as to support satellite and climate model programs over the WFF site via a bilateral approach of storm-scale observations and process modeling. This study highlights a noble systematic validation framework of the BLOSSOM ensemble cloud-process simulations through mixed-phase, light-rain, and deep-convective precipitation cases. The framework consists of creating a domain-shifted ensemble of large-scale forcing datasets, and configuring and performing cloud-process simulations with three different bulk microphysics schemes. Validation uses NASA S-band dual-POLarimetric radar (NPOL) observations in the form of statistical composites and skill scores via a polarimetric radar simulator and newly developed CfRad Data tool (CfRAD). While the simulations capture the overall structures of the reflectivity composites, polarimetric signals are still poorly simulated, mainly due to a lack of representation of ice microphysics diversity in shapes, orientation distributions, and their complex mixtures. Despite the limitation, this new ensemble-based validation framework demonstrates that 1) no particular forcing or microphysics scheme outperforms the rest and 2) the skill scores of coarse- and fine-resolution ensemble simulations with different domain-shifted forcing and microphysics schemes are highly correlated with each other with no clear improvement. On the other hand, this suggests that coarse-resolution ensemble simulations are relevant for selecting the best meteorological forcing and microphysics scheme before conducting computationally demanding large eddy simulations (LESs) in support of aircraft and satellite instrument development as well as cloud-precipitation-convection parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Trapping of spermine, Kukoamine A, and polyamine toxin blockers in GluK2 kainate receptor channels

Abstract Kainate receptors (KARs) are a subtype of ionotropic glutamate receptor (iGluR) channels, a superfamily of ligand-gated ion channels which mediate the majority of excitatory neurotransmission in the central nervous system. KARs modulate neuronal circuits and plasticity during development and are implicated in neurological disorders, including epilepsy, depression, schizophrenia, anxiety, and autism. Calcium-permeable KARs undergo ion channel block, but the therapeutic potential of channel blockers remains underdeveloped, mainly due to limited structural knowledge. Here, we present closed-state structures of GluK2 KAR homotetramers in complex with ion channel blockers NpTx-8, PhTx-74, Kukoamine A, and spermine. We find that blockers reside inside the GluK2 ion channel pore, intracellular to the closed M3 helix bundle-crossing gate, with their hydrophobic heads filling the central cavity and positively charged polyamine tails spanning the selectivity filter. Molecular dynamics (MD) simulations of our structures illuminate interactions responsible for different affinity and binding poses of the blockers. Our structures elucidate the trapping mechanism of KAR channel block and provide a template for designing new blockers that can selectively target calcium-permeable KARs in neuropathologies.

Science & Technology - Other Topics↗

Anisotropic X-Ray Scattering of Transiently Oriented Water

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ~ 0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗

Pulsar Simulations for the Fermi Large Area Telescope

Pulsars are among the prime targets for the Large Area Telescope (LAT) aboard the recently launched Fermi observatory. The LAT will study the gamma-ray Universe between 20 MeV and 300 GeV with unprecedented detail. Increasing numbers of gamma-ray pulsars are being firmly identified, yet their emission mechanisms are far from being understood. To better investigate and exploit the tAT capabilities for pulsar science. a set of new detailed pulsar simulation tools have been developed within the LAT collaboration. The structure of the pulsar simulator package (PulsarSpeccrum) is presented here. Starting from photon distributions in energy and phase obtained from theoretical calculations or phenomenological considerations, gamma-rays are generated and their arrival times at the spacecraft are determined by taking Into account effects such as barycentric effects and timing noise. Pulsars in binary systems also can be simulated given orbital parameters. We present how simulations can be used for generating a realistic set of gamma rays as observed by the LAT, focusing on some case studies that show the performance of the LAT for pulsar observations.

Razzano, M.↗

Phase transitions and 2D spherical convection in a large icy satellite

Three ice phase transitions, based on experimental and theoretical results, were incorporated into a 2D spherical convection model. These phase transitions are the exothermic Ice 1-2, the endothermic Ice 2-6, and the exothermic Ice 6-8 phase transitions. The fluid is assumed isoviscous and the anelastic liquid formulation was used. The object is an icy satellite whose physical characteristics resemble those of Ganymede or Callisto. Structural models, similar to those of Mueller and McKinnon, with different core sizes, i.e. with different degree of differentiation were studied. The shells, i.e. the mantle of the satellite, is heated both internally and from below to account for the decaying radiogenic heating and the heat flow from the solid core. The lower boundary of the mantle is rigid and isothermal, the upper boundary is isothermal. Calculations with rigid and shear stress free upper boundary are carried out in order to assess the role of the different boundary conditions. Two different Rayleigh numbers, depending on the assumed value of the viscosity, were used in the calculations and the thermal evolution of the satellite was studied. The suit of calculations presented demonstrates that phase transition cannot be ignored when the thermal evolution of a large icy satellite has to be studied. The importance of the construction of the realistic phase diagram for ice was also documented. On the other hand the degree of differentiation, simulated with the structures with different core's radius does not seem to drastically influence the overall behavior of the convective activity. Within the limits of our simulations, the post-accretional evolution seems to be independent on the degree of primordial differentiation.

Forni, O.↗

Divergence in land surface modeling: linking spread to structure

Divergence in land carbon cycle simulation is persistent and widespread. Regardless of model intercomparison project, results from individual models diverge significantly from each other and, in consequence, from reference datasets. Here we link model spread to structure using a 15-member ensemble of land surface models from the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) as a test case. Our analysis uses functional benchmarks and model structure as predicted by model skill in a machine learning framework to isolate discrete aspects of model structure associated with divergence. We also quantify how initial conditions prejudice present-day model outcomes after centennial-scale transient simulations. Overall, the functional benchmark and machine learning exercises emphasize the importance of ecosystem structure in correctly simulating carbon and water cycling, highlight uncertainties in the structure of carbon pools, and advise against hard parametric limits on ecosystem function. We also find that initial conditions explain 90% of variation in global satellite-era values—initial conditions largely predetermine transient endpoints, historical environmental change notwithstanding. As MsTMIP prescribes forcing data and spin-up protocol, the range in initial conditions and high levels of predetermination are also structural. Our results suggest that methodological tools linking divergence to discrete aspects of model structure would complement current community best practices in model development.

58 GEOSCIENCES↗

A New Standard DNA Damage (SDD) Data Format

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called “indirect” damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates intermodel comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to significantly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Simulating seawater intrusion in a complex coastal karst aquifer using an improved variable-density flow and solute transport–conduit flow process model

VDFST-CFP (variable-density flow and solute transport–conduit flow process) is a density-dependent discrete-continuum numerical model for simulating seawater intrusion in a dual-permeability coastal karst aquifer. A previous study (Xu and Hu 2017) simulates variable-density flow only in a single conduit, and studies the parameter sensitivities only in the horizontal case (2D domain as horizontal section) by the VDFST-CFP model. This paper focuses on the density-dependent vertical case (2D domain as vertical section) with two major improvements: 1) when implementing double-conduit networks in the domain, simulated intruded plumes in the porous medium are extended in the double-conduit scenario, compared to the single-conduit system; 2) by quantifying micro-textures on the conduit wall by the Goudar-Sonnad equation and considering macro-structures as local head loss. Sensitivity analysis shows that medium hydraulic conductivity, conduit diameter and effective porosity are important parameters for simulating seawater intrusion in the discrete-continuum system. On the other hand, rougher micro-structures and additional macro-structure components on the conduit wall would reduce the distance of seawater intrusion to the conduit system, but, rarely affect salinity distribution in the matrix. Compared to the equivalent mean roughness height, the new method (with more detailed description of structure) simulates seawater intrusion slightly landward in the conduit system. Here, the macro-structure measured by local head loss is more reasonable but needs further study on conduit flow. Xu and Hu (2017) Development of a discrete-continuum VDFST-CFP numerical model for simulating seawater intrusion to a coastal karst aquifer with a conduit system.

58 GEOSCIENCES↗

Wall-Modeled Large-Eddy Simulation of Supersonic Parachute Inflation

Supersonic parachutes have been used in nearly every robotic mission to another planetary body with an atmosphere because they are one of the most mass-efficient ways to decelerate a payload to land on the surface. Short of performing flight tests in the upper Earth atmosphere, we currently cannot reliably predict a novel parachute system’s performance or potential failure modes, or even confidently explain it after the fact, as was the case with the Low Density Supersonic Decelerator flight tests and subsequent investigations. Fluid-structure interaction (FSI) simulations have the potential to bridge this gap. To this end, we present several improvements to the state-of-the-art for simulating supersonic parachutes using FSI: a higher effective resolution convective flux, an immersed boundary turbulent wall layer modeling approach to capture viscous effects, and a novel method to obtain a more realistic initial parachute shape. A recent supersonic parachute flight test is simulated and compared to measurements for the purposes of model validation.

ESM↗

Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: a benchmark Quantum Monte Carlo study

The accurate description of non-covalent interactions is critical for understanding the structure, dynamics, and eventual function of biomolecules. The adenine dimer serves as a benchmark system for computational methods due to its role in nucleic acid structures and its rich conformational landscape. In this study, we employ benchmark diffusion quantum Monte Carlo (DMC) methods to investigate the relative energies and role of electron correlation on a set of adenine dimer conformations generated via a search of the potential energy landscape using the global optimizer algorithm. Relative DMC energies are compared against a wide range of density functional theory (DFT) approximation results. We find that although most of the DFT functionals perform well for low-energy structures, their accuracy varies significantly for higher-energy conformations, including stacked and T-shaped structures. A large fraction of the variation is due to the treatment of the van der Waals interaction. BLYP, B3LYP, and PBE0 significantly improve with added D4 dispersion, while the recent r2SCAN-D4 and ωB97M-V functionals show the least scatter and closest agreement with the DMC. These findings highlight the delicate nature of these interactions in biomolecular systems and provide guidance for simulations of their structure and dynamics and for the development of machine learned interatomic potentials.

Washburn, Laurel [ORNL] (ORCID:0000000324179335)↗

Sensitivity of High-Resolution Simulations of Hurricane Bob (1991) to Planetary Boundary Layer Parameterizations

The MM5 mesoscale model is used to simulate Hurricane Bob (1991) using grids nested to high resolution (4 km). Tests are conducted to determine the sensitivity of the simulation to the available planetary boundary layer parameterizations, including the bulk-aerodynamic, Blackadar, Medium-RanGe Forecast (MRF) model, and Burk-Thompson boundary-layer schemes. Significant sensitivity is seen, with minimum central pressures varying by up to 17 mb. The Burk-Thompson and bulk-aerodynamic boundary-layer schemes produced the strongest storms while the MRF scheme produced the weakest storm. Precipitation structure of the simulated hurricanes also varied substantially with the boundary layer parameterizations. Diagnostics of boundary-layer variables indicated that the intensity of the simulated hurricanes generally increased as the ratio of the surface exchange coefficients for heat and momentum, C(sub h)/C(sub M), although the manner in which the vertical mixing takes place was also important. Findings specific to the boundary-layer schemes include: 1) the MRF scheme produces mixing that is too deep and causes drying of the lower boundary layer in the inner-core region of the hurricane; 2) the bulk-aerodynamic scheme produces mixing that is probably too shallow, but results in a strong hurricane because of a large value of C(sub h)/C(sub M) (approximately 1.3); 3) the MRF and Blackadar schemes are weak partly because of smaller surface moisture fluxes that result in a reduced value of C(sub h)/C(sub M) (approximately 0.7); 4) the Burk-Thompson scheme produces a strong storm with C(sub h)/C(sub M) approximately 1; and 5) the formulation of the wind-speed dependence of the surface roughness parameter, z(sub 0), is important for getting appropriate values of the surface exchange coefficients in hurricanes based upon current estimates of these parameters.

Braun, Scott A.↗

Multi-Task Learning of Scanning Electron Microscopy and Synthetic Thermal Tomography Images for Detection of Defects in Additively Manufactured Metals

One of the key challenges in laser powder bed fusion (LPBF) additive manufacturing of metals is the appearance of microscopic pores in 3D-printed metallic structures. Quality control in LPBF can be accomplished with non-destructive imaging of the actual 3D-printed structures. Thermal tomography (TT) is a promising non-contact, non-destructive imaging method, which allows for the visualization of subsurface defects in arbitrary-sized metallic structures. However, because imaging is based on heat diffusion, TT images suffer from blurring, which increases with depth. We have been investigating the enhancement of TT imaging capability using machine learning. In this work, we introduce a novel multi-task learning (MTL) approach, which simultaneously performs the classification of synthetic TT images, and segmentation of experimental scanning electron microscopy (SEM) images. Synthetic TT images are obtained from computer simulations of metallic structures with subsurface elliptical-shaped defects, while experimental SEM images are obtained from imaging of LPBF-printed stainless-steel coupons. MTL network is implemented as a shared U-net encoder between the classification and the segmentation tasks. Results of this study show that the MTL network performs better in both the classification of synthetic TT images and the segmentation of SEM images tasks, as compared to the conventional approach when the individual tasks are performed independently of each other.

36 MATERIALS SCIENCE↗