Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural dynamics”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

On the Use of Vortex-Fitting in the Numerical Simulation of Blade-Vortex Interaction

The usefulness of vortex-fitting in the computational fluid dynamics (CFD) methods to preserve the vortex strength and structure while convecting in a uniform free stream is demonstrated through the numerical simulations of two- and three-dimensional blade-vortex interactions. The fundamental premise of the formulation is the velocity and pressure field of the interacting vortex are unaltered either in the presence of an airfoil or a rotor blade or by the resulting nonlinear interactional flowfield. Although, the governing Euler and Navier-Stokes equations are nonlinear and independent solutions cannot be superposed, the interactional flowfield can be accurately captured by adding and subtracting the flowfield of the convecting vortex at each instant. The aerodynamics and aeroacoustics of two- and three-dimensional blade-vortex interactions have been calculated in Refs. 1-6 using this concept. Some of the results from these publications and similar other published material will be summarized in this paper.

Srinivasan, G. R.↗
NASA NTRS2006

Applying Parallel Adaptive Methods with GeoFEST/PYRAMID to Simulate Earth Surface Crustal Dynamics

This viewgraph presentation reviews the use Adaptive Mesh Refinement (AMR) in simulating the Crustal Dynamics of Earth's Surface. AMR simultaneously improves solution quality, time to solution, and computer memory requirements when compared to generating/running on a globally fine mesh. The use of AMR in simulating the dynamics of the Earth's Surface is spurred by future proposed NASA missions, such as InSAR for Earth surface deformation and other measurements. These missions will require support for large-scale adaptive numerical methods using AMR to model observations. AMR was chosen because it has been successful in computation fluid dynamics for predictive simulation of complex flows around complex structures.

GeoFEST↗
NASA NTRS2010

Boundary Layer Protuberance Simulations in Channel Nozzle Arc-Jet

Two protuberance designs were modeled in the channel nozzle of the NASA Johnson Space Center Atmospheric Reentry Materials and Structures Facility with the Data-Parallel Line Relaxation computational fluid dynamics code. The heating on the protuberance was compared to nominal baseline heating at a single fixed arc-jet condition in order to obtain heating augmentation factors for flight traceability in the Boundary Layer Transition Flight Experiment on Space Shuttle Orbiter flights STS-119 and STS-128. The arc-jet simulations were performed in conjunction with the actual ground tests performed on the protuberances. The arc-jet simulations included non-uniform inflow conditions based on the current best practices methodology and used variable enthalpy and constant mass flow rate across the throat. Channel walls were modeled as fully catalytic isothermal surfaces, while the test section (consisting of Reaction Cured Glass tiles) was modeled as a partially catalytic radiative equilibrium wall. The results of the protuberance and baseline simulations were compared to the applicable ground test results, and the effects of the protuberance shock on the opposite channel wall were investigated.

Marichalar, J. J.↗
NASA NTRS2016

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗
NASA NTRS2016

SUPIN: A Computational Tool for Supersonic Inlet Design

A computational tool named SUPIN is being developed to design and analyze the aerodynamic performance of supersonic inlets. The inlet types available include the axisymmetric pitot, three-dimensional pitot, axisymmetric outward-turning, two-dimensional single-duct, two-dimensional bifurcated-duct, and streamline-traced inlets. The aerodynamic performance is characterized by the flow rates, total pressure recovery, and drag. The inlet flow-field is divided into parts to provide a framework for the geometry and aerodynamic modeling. Each part of the inlet is defined in terms of geometric factors. The low-fidelity aerodynamic analysis and design methods are based on analytic, empirical, and numerical methods which provide for quick design and analysis. SUPIN provides inlet geometry in the form of coordinates, surface angles, and cross-sectional areas. SUPIN can generate inlet surface grids and three-dimensional, structured volume grids for use with higher-fidelity computational fluid dynamics (CFD) analysis. Capabilities highlighted in this paper include the design and analysis of streamline-traced external-compression inlets, modeling of porous bleed, and the design and analysis of mixed-compression inlets. CFD analyses are used to verify the SUPIN results.

computational fluid dynamics↗
NASA NTRS2017

Thermosphere Global Time Response to Geomagnetic Storms Caused by Coronal Mass Ejections

We investigate, for the first time with a spatial superposed epoch analysis study, the thermosphere global time response to 159 geomagnetic storms caused by coronal mass ejections (CMEs) observed in the solar wind at Earth's orbit during the period of September 2001 to September 2011. The thermosphere neutral mass density is obtained from the CHAMP (CHAllenge Mini-Satellite Payload) and GRACE (Gravity Recovery Climate Experiment) spacecraft. All density measurements are intercalibrated against densities computed by the Jacchia-Bowman 2008 empirical model under the regime of very low geomagnetic activity. We explore both the effects of the pre-CME shock impact on the thermosphere and of the storm main phase onset by taking their times of occurrence as zero epoch times (CME impact and interplanetary magnetic field Bz southward turning) for each storm. We find that the shock impact produces quick and transient responses at the two high-latitude regions with minimal propagation toward lower latitudes. In both cases, thermosphere is heated in very high latitude regions within several minutes. The Bz southward turning of the storm onset has a fast heating manifestation at the two high-latitude regions, and it takes approximately 3 h for that heating to propagate down to equatorial latitudes and to globalize in the thermosphere. This heating propagation is presumably accomplished, at least in part, with traveling atmospheric disturbances and complex meridional wind structures. Current models use longer lag times in computing thermosphere density dynamics during storms. Our results suggest that the thermosphere response time scales are shorter and should be accordingly adjusted in thermospheric empirical models.

mass ejections↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

Science & Technology - Other Topics↗
DOE OSTI2025

Selective mass scaling for single-layer thick shell elements in DYNA3D

Hexahedral elements can be adapted to model thin and moderately thick structures by neglecting the coupling of through-thickness stress, resulting in a fully three-dimensional, but simplified, state of stress. These specialized elements, often referred to as “thick” or “solid” shells, are generally employed to model thin-walled structures using continuum mechanics-based material models. In explicit dynamics simulations, where computational speed is important, these elements are integrated with a single quadrature point and a set of anti-hourglassing (stabilizing) forces. Thick shells, by definition, have a thickness dimension smaller than their in-plane dimensions, and this small thickness often determines the stable time step size in simulations, despite the mechanics being approximated. To alleviate this limitation while retaining the relevant dynamics of thin-walled structures, selective mass scaling (SMS), or selective mass “augmentation,” has been proposed in the literature. In this technical report, we explore the application of SMS to single-layer thick shells in the simulation software DYNA3D.

42 ENGINEERING↗
NASA NTRS2011

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗
DOE OSTI2021

Exceptional Radiation Absorption in a Pentagon-Based Si Allotrope

Excellent photovoltaic performance is predicted in a pentagonal covalent network of Si in a hollow structure exhibiting both thermal and dynamical stability. Consisting of a combination of sp 2 and sp 3 hybridized Si atomic orbitals, the GW0 computed band structure shows an indirect band gap near the zone edge and also a manifold of directly absorbing transitions at frequencies in the window of visible light, in distinction with conventional Si. Hydrogenation of a single sp 2 site is predicted to lead to a robust local magnetic moment. In this study, we find a low formation energy at low pressure that is compatible with other experimentally known phases, suggesting that a stable phase might be obtained.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Computational insights into the lattice dynamics of Pu(IV) oxalates

Despite its use in PuO 2 production, the structure of anhydrous Pu(C 2 O 4 ) 2 is still not completely understood. Recently, two candidate structures for Pu(C 2 O 4 ) 2 were proposed via density functional theory (DFT), after which the first experimental optical vibrational spectra were reported. Here, in this study, we calculated the lattice dynamics of the candidate structures using DFT and found that the primary difference between them is the presence of a vibrational mode near 1380 cm -1 in one structure. The frequency and optical activity of this mode agree well with the published experimental results, providing strong support for this calculated structure as that of anhydrous Pu(C 2 O 4 ) 2 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Structural and Dynamical Design Principles for Iron-Sulfur Clusters in Sequence Defined Polymers

This integrated synthetic and computational effort uses sequence defined polymers such as peptoids, and triazines, capable of mimicking a flexible protein environment to control the positioning of side chains around redox cofactors. This strategy will allow modulation of redox properties through both direct and indirect mechanisms. To this end the aim is defined to test the underlying hypothesis: that principles and critical functionalities governing the binding and tuning of naturally occurring electron carrier co-factors can be identified and used to predictably design and synthesize maquettes based on non-natural polymers resulting in new modular biomimetic scaffolds for Electron Transfer Materials. This will be achieved through design and synthesis of ligands that can stabilize iron-sulfur clusters and tune the redox potential based on sequence.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Dynamic response analysis of geometrically non-linear structures subjected to high impact.

Description of an efficient digital computer method for the determination of the propagation of elastic stresses and deformations in certain geometrically nonlinear structures subjected to high impact loading. The finite-element matrix displacement approach utilizing curved quadrilateral shell elements in conjunction with a nodewise predictor-corrector method employing Runge-Kutta extrapolation techniques has been adopted for the present solution. The related computer program written in FORTRAN V for the UNIVAC 1108 computer has proved to be effective for the solution of a range of practical problems including rectangular and cylindrical panels. Numerical results are presented for a relevant structure, the cell container, and the negative electrode of an impact-resistant battery subjected to high impact, simulating its free landing on a planetary surface.

Gupta, K. K.↗
NASA NTRS2001

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
NASA NTRS2017

Realistic Modeling of Interaction of Quiet-Sun Magnetic Fields with the Chromosphere

High­‐resolution observations and 3D MHD simulations reveal intense interaction between the convection zone dynamics and the solar atmosphere on subarcsecond scales. To investigate processes of the dynamical coupling and energy exchange between the subsurface layers and the chromosphere we perform 3D radiative MHD modeling for a computational domain that includes the upper convection zone and the chromosphere, and investigate the structure and dynamics for different intensity of the photospheric magnetic flux. For comparison with observations, the simulation models have been used to calculate synthetic Stokes profiles of various spectral lines. The results show intense energy exchange through small­‐scale magnetized vortex tubes rooted below the photosphere, which provide extra heating of the chromosphere, initiate shock waves, and small­‐scale eruptions.

quiet-sun↗