SEARCH · Engineering Papers
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.
A model for the intermittent fine structure of high Reynolds number turbulence is proposed. The model consists of slender axially strained spiral vortex solutions of the Navier-Stokes equation. The tightening of the spiral turns by the differential rotation of the induced swirling velocity produces a cascade of velocity fluctuations to smaller scale. The Kolmogorov energy spectrum is a result of this model.
Since the 14th Critical Assessment of Techniques for Protein Structure Prediction (CASP14), AlphaFold2 has become the standard method for protein tertiary structure prediction. One remaining challenge is to further improve its prediction. We developed a new version of the MULTICOM system to sample diverse multiple sequence alignments (MSAs) and structural templates to improve the input for AlphaFold2 to generate structural models. The models are then ranked by both the pairwise model similarity and AlphaFold2 self-reported model quality score. The top ranked models are refined by a novel structure alignment-based refinement method powered by Foldseek. Moreover, for a monomer target that is a subunit of a protein assembly (complex), MULTICOM integrates tertiary and quaternary structure predictions to account for tertiary structural changes induced by protein-protein interaction. The system participated in the tertiary structure prediction in 2022 CASP15 experiment. Our server predictor MULTICOM_refine ranked 3rd among 47 CASP15 server predictors and our human predictor MULTICOM ranked 7th among all 132 human and server predictors. The average GDT-TS score and TM-score of the first structural models that MULTICOM_refine predicted for 94 CASP15 domains are ~0.80 and ~0.92, 9.6% and 8.2% higher than ~0.73 and 0.85 of the standard AlphaFold2 predictor respectively.
The structural models presently defined for Ganymede and Callisto, which encompass a pure-ice upper layer, a mixed ice/rock lower mantle, and a rock core, incorporate three alternative rock component candidates representing various degrees of silicate hydration and oxidation. The three-layered model facilitates close study of the radius increase required for the internal differentiation of an ice-rock satellite; such expansion is determined to be most significant early in the process, and less so as differentiation approaches completion. The probability of postaccretional melting due to radiogenic heating is calculated.
To enable thermal behavior prediction and design optimization of molten salt reactors, thermal conductivity of molten salts must be characterized in terms of salt composition and temperature. Current theoretical models fail to provide consistent approximations for all halide mixtures, particularly actinide-bearing melts. This study aims to link the short-range order structure of molten salts to the mean free path of energy carriers through a simple structural coherence model informed by the partial pair distribution function. The proposed method is used to predict the thermal conductivity of 33 alkali and alkaline earth halide salts. Predictions approximate experimental measurements with a mean absolute error of 15.7% for dissociating, complexing, and actinide salts, including unary LiCl, NaCl, and MgCl 2 as well as mixtures LiF–NaF–KF (FLiNaK), LiF–BeF 2 (FLiBe), and NaCl–UCl 3 . The work provides evidence for the validity of energy carrier descriptions of molecular-level heat transfer in molten salts, with implications for improved theories of liquid energy transport in general.
The semianalytical method, widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements, is studied in this paper. The paper shows that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. The errors are shown to be associated with the structural model. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. Moreover, a method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.
This report presents the description for the NASTRAN finite element for the AMSU-A1 module. The purpose of this report is to document the NASTRAN bulk data deck, transmitted under separate cover. The structural Math Model is to be used by the spacecraft contractor for dynamic loads analysis.
Analytical and experimental techniques for the prediction and ground verification of the damped structural dynamics of space structures are developed. The options available for similarity-scaled model testing, including replica and multiple scale approaches, are reviewed. For the case when the distortion of potentially dissipative or nonlinear joints, which would be required in multiple-scale modeling, is impractical, a new type of modeling is introduced, which uses a hybrid of joints at replica scale and connecting elements at a modified multiple scale. The model design requirements for replica, multiple-scale, and hybrid models are developed, and the expected scaling of nonlinear dissipation in joints is derived. A damping prediction scheme is developed that relies on a finite element model of the undamped structure and measurements of the individual joint properties to predict the modal damping of the truss attributable to the joints. A hybrid-scaled model of a segment of the Space Station was built and dynamically tested. The predicted and measured truss damping compared favorably.
The drones for aerodynamic and structural testing (DAST) project was designed to control flutter actively at high subsonic speeds. Accurate knowledge of the structural model was critical for the successful design of the control system. A ground vibration test was conducted on the DAST vehicle to determine the structural model characteristics. This report presents and discusses the vibration and test equipment, the test setup and procedures, and the antisymmetric and symmetric mode shape results. The modal characteristics were subsequently used to update the structural model employed in the control law design process.
The design, construction and instrumentation of plastic model shells used in elastic stability studies is discussed. The work demonstrates that specimens made of such materials are of great value in parametric studies because of their ease of manufacture and the simplicity of modification. Additionally, it is shown that they are most useful in the development of nondestructive test techniques. The paper indicates how deformation measurements made at low load levels can be used to accurately predict critical conditions.
Studies have shown that significant vehicle mass and cost savings are possible with the use of ballutes for aero-capture. Through NASA's In-Space Propulsion program, a preliminary examination of ballute sensitivity to geometry and Reynolds number was conducted, and a single-pass coupling between an aero code and a finite element solver was used to assess the static aeroelastic effects. There remain, however, a variety of open questions regarding the dynamic aeroelastic stability of membrane structures for aero-capture, with the primary challenge being the prediction of the membrane flutter onset. The purpose of this paper is to describe and begin addressing these issues. The paper includes a review of the literature associated with the structural analysis of membranes and membrane utter. Flow/structure analysis coupling and hypersonic flow solver options are also discussed. An approach is proposed for tackling this problem that starts with a relatively simple geometry and develops and evaluates analysis methods and procedures. This preliminary study considers a computationally manageable 2-dimensional problem. The membrane structural models used in the paper include a nonlinear finite-difference model for static and dynamic analysis and a NASTRAN finite element membrane model for nonlinear static and linear normal modes analysis. Both structural models are coupled with a structured compressible flow solver for static aeroelastic analysis. For dynamic aeroelastic analyses, the NASTRAN normal modes are used in the structured compressible flow solver and 3rd order piston theories were used with the finite difference membrane model to simulate utter onset. Results from the various static and dynamic aeroelastic analyses are compared.
A new structural modeling and analysis method, the Distributed Transfer Function Method (DTFM), is presented for the applications of space gossamer structures.
A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.
A technique is presented to use a set of identified complex modes together with an analytical mathematical model of a structure under test to compute improved mass, stiffness and damping matrices. A set of identified normal modes, computed from the measured complex modes, is used in the mass orthogonality equation to compute an improved mass matrix. This eliminates possible errors that may result from using approximated complex modes as normal modes. The improved mass matrix, the measured complex modes and the higher analytical modes are then used to compute the improved stiffness and damping matrices. The number of degrees-of-freedom of the improved model is limited to equal the number of elements in the measured modal vectors. A simulated experiment shows considerable improvements, in the system's analytical dynamic model, over the frequency range of the given measured modal information.
Root system architecture (RSA) is an important measure of how plants navigate and interact with the soil environment. However, current methods in studying RSA must make tradeoffs between precision of data and proximity to natural conditions, with root growth in germination papers providing accessibility and high data resolution. Functional–structural plant models (FSPMs) can overcome this tradeoff, though parameterization and evaluation of FSPMs are traditionally based in manual measurements and visual comparison. Here, we applied a germination paper system to study the adventitious RSA and root phenology of Populus trichocarpa stem cuttings using time-series image-based phenotyping augmented by FSPM. We found a significant correlation between timing of root initiation and thermal time at cutting collection (P value = 0.0061, R 2 = 0.875), but little correlation with RSA. We also present a use of RhizoVision [1] for automatically extracting FSPM parameters from time series images and evaluating FSPM simulations. A high accuracy of the parameterization was achieved in predicting 2D growth with a sensitivity rate of 83.5%. This accuracy was lost when predicting 3D growth with sensitivity rates of 38.5% to 48.7%, while overall accuracy varied with phenotyping methods. Despite this loss in accuracy, the new method is amenable to high throughput FSPM parameterization and bridges the gap between advances in time-series phenotyping and FSPMs.
The dynamics of molecules exposed to intense x-ray radiation involve a large number of multiply ionized and highly excited electronic configurations. To model these dynamics a reliable and efficient electronic structure model is imperative. Employing the Hartree-Fock-Slater electronic structure model in combination with the maximum overlap method, we quantify the associated convergence failures when calculating electronic states of carbon monoxide with multiple vacancies in the core and valence levels. Here, we characterize these cases and describe strategies to overcome the convergence problems. The described techniques not only eliminate all convergence issues for CO but also result in a significant reduction of convergence failures for simulations of the x-ray-induced multiple ionization dynamics of the phenol molecule.
Demanding aerodynamic modelling requirements for military and civilian aircraft have motivated researchers to improve computational and experimental techniques and to pursue closer collaboration in these areas. Model identification and validation techniques are key components for this research. This paper presents mathematical model structures and identification techniques that have been used successfully to model more general aerodynamic behaviours in single-degree-of-freedom dynamic testing. Model parameters, characterizing aerodynamic properties, are estimated using linear and nonlinear regression methods in both time and frequency domains. Steps in identification including model structure determination, parameter estimation, and model validation, are addressed in this paper with examples using data from one-degree-of-freedom dynamic wind tunnel and water tunnel experiments. These techniques offer a methodology for expanding the utility of computational methods in application to flight dynamics, stability, and control problems. Since flight test is not always an option for early model validation, time history comparisons are commonly made between computational and experimental results and model adequacy is inferred by corroborating results. An extension is offered to this conventional approach where more general model parameter estimates and their standard errors are compared.
The sensitivity of photosynthesis to environmental changes is essential for understanding carbon cycle responses to global climate change and for the development of modeling approaches that explains its spatial and temporal variability. We collected a large variety of published sensitivity functions of gross primary productivity (GPP) to different forcing variables to assess the response of GPP to environmental factors. These include the responses of GPP to temperature; vapor pressure deficit, some of which include the response to atmospheric CO 2 concentrations; soil water availability (W); light intensity; and cloudiness. These functions were combined in a full factorial light use efficiency (LUE) model structure, leading to a collection of 5600 distinct LUE models. Each model was optimized against daily GPP and evapotranspiration fluxes from 196 FLUXNET sites and ranked across sites based on a bootstrap approach. The GPP sensitivity to each environmental factor, including CO 2 fertilization, was shown to be significant, and that none of the previously published model structures performed as well as the best model selected. From daily and weekly to monthly scales, the best model's median Nash-Sutcliffe model efficiency across sites was 0.73, 0.79 and 0.82, respectively, but poorer at annual scales (0.23), emphasizing the common limitation of current models in describing the interannual variability of GPP. Although the best global model did not match the local best model at each site, the selection was robust across ecosystem types. The contribution of light saturation and cloudiness to GPP was observed across all biomes (from 23% to 43%). Temperature and W dominates GPP and LUE but responses of GPP to temperature and W are lagged in cold and arid ecosystems, respectively. The findings of this study provide a foundation towards more robust LUE-based estimates of global GPP and may provide a benchmark for other empirical GPP products.
Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.