SEARCH · Search NASA
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.
Parallel adaptive mesh refinement (AMR) is an important numerical technique that leads to the efficient solution of many physical and engineering problems. In this paper, we describe how AMR programing can be performed in an object-oreinted way using the modern aspects of Fortran 90 combined with the parallelization features of OpenMP.
Computational methods for global/local analysis of structures which include both uncoupled and coupled methods are described. In addition, global/local analysis methodology for automatic refinement of incompatible global and local finite element models is developed. Representative structural analysis problems are presented to demonstrate the global/local analysis methods.
High performance control design for a flexible space structure is challenging since high fidelity plant models are difficult to obtain a priori. Uncertainty in the control design models typically require a very robust, low performance control design which must be tuned on-orbit to achieve the required performance. A new procedure for refining a multivariable open loop plant model based on closed-loop response data is presented. Using a minimal representation of the state space dynamics, a least squares prediction error method is employed to estimate the plant parameters. This control-relevant system identification procedure stresses the joint nature of the system identification and control design problem by seeking to obtain a model that minimizes the difference between the predicted and actual closed-loop performance. This paper presents an algorithm for iterative closed-loop system identification and controller redesign along with illustrative examples.
Research for this grant focused on the following tasks: (1) the prediction of severe, localized aerodynamic heating for complex, high speed flows; (2) finite element adaptive refinement methodology for multi-disciplinary analyses; (3) the prediction of thermoviscoplastic structural response with rate-dependent effects and large deformations; (4) thermoviscoplastic constitutive models for metals; and (5) coolant flow/structural heat transfer analyses.
Parallel Adaptive Mesh Refinement Library (PARAMESH) is a package of Fortran 90 subroutines designed to provide a computer programmer with an easy route to extension of (1) a previously written serial code that uses a logically Cartesian structured mesh into (2) a parallel code with adaptive mesh refinement (AMR). Alternatively, in its simplest use, and with minimal effort, PARAMESH can operate as a domain-decomposition tool for users who want to parallelize their serial codes but who do not wish to utilize adaptivity. The package builds a hierarchy of sub-grids to cover the computational domain of a given application program, with spatial resolution varying to satisfy the demands of the application. The sub-grid blocks form the nodes of a tree data structure (a quad-tree in two or an oct-tree in three dimensions). Each grid block has a logically Cartesian mesh. The package supports one-, two- and three-dimensional models.
Electronic integration of design and analysis processes was achieved and refined at Langley Research Center (LaRC) during the development of an optical bench for a laser-based aerospace experiment. Mechanical design has been integrated with thermal, structural and optical analyses. Electronic import of the model geometry eliminates the repetitive steps of geometry input to develop each analysis model, leading to faster and more accurate analyses. Guidelines for integrated model development are given. This integrated analysis process has been built around software that was already in use by designers and analysis at LaRC. The process as currently implemented used Pro/Engineer for design, Pro/Manufacturing for fabrication, PATRAN for solid modeling, NASTRAN for structural analysis, SINDA-85 and P/Thermal for thermal analysis, and Code V for optical analysis. Currently, the only analysis model to be built manually is the Code V model; all others can be imported for the Pro/E geometry. The translator from PATRAN results to Code V optical analysis (PATCOD) was developed and tested at LaRC. Directions for use of the translator or other models are given.
We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.
The crystal structure of brimonidine hydrogen tartrate has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional techniques. Brimonidine hydrogen tartrate crystallizes in space groupP2 1 (#4) witha= 7.56032(2),b= 7.35278(2),c= 30.10149(9) Å,β= 90.1992(2)°,V= 1673.312(10) Å 3 , andZ= 4 at 295 K. The crystal structure consists of alternating layers of cations and anions parallel to theab-plane. Each of the hydrogen tartrate anions is linked to itself by very strong charge-assisted O–H⋯O hydrogen bonds into chains along thea-axis. Each hydroxyl group of each tartrate acts as a donor in an O–H⋯O or O–H⋯N hydrogen bond. One of these is intramolecular, but the other three are intermolecular. These hydrogen bonds link the hydrogen tartrate anions into layers parallel to theab-plane and also link the anion–cation layers. The protonated N atoms act as donors in N–H⋯O or N–H⋯N hydrogen bonds to the carboxyl groups of the tartrates and to a ring nitrogen atom. These link the cations and anions, as well as providing cation–cation links. The amino N atoms of the cations form N–H⋯O hydrogen bonds to hydroxyl groups of the anions. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)
The analytical and experimental results of an all-composite wing stub box are presented in this report. The wing stub box, which is representative of an inboard portion of a commercial transport high-aspect-ratio wing, was fabricated from stitched graphite-epoxy material with a Resin Film Infusion manufacturing process. The wing stub box was designed and constructed by the McDonnell Douglas Aerospace Company as part of the NASA Advanced Composites Technology program. The test article contained metallic load-introduction structures on the inboard and outboard ends of the graphite-epoxy wing stub box. The root end of the inboard load introduction structure was attached to a vertical reaction structure, and an upward load was applied to the outermost tip of the outboard load introduction structure to induce bending of the wing stub box. A finite element model was created in which the center portion of the wing-stub-box upper cover panel was modeled with a refined mesh. The refined mesh was required to represent properly the geometrically nonlinear structural behavior of the upper cover panel and to predict accurately the strains in the stringer webs of the stiffened upper cover panel. The analytical and experimental results for deflections and strains are in good agreement.
The potential is explored for using composite materials in urban mass transit systems. The emphasis was to identify specific advantages of composite materials in order to determine their actual and potential usage for carbody and guideway structure applications. The literature was reviewed, contacts were made with major domestic system operators, designers, and builders, and an analysis was made of potential composite application to railcar construction. Composites were found to be in use throughout the transit industry, usually in secondary or auxiliary applications such as car interior and nonstructural exterior panels. More recently, considerable activity has been initiated in the area of using composites in the load bearing elements of civil engineering structures such as highway bridges. It is believed that new and improved manufacturing refinements in pultrusion and filament winding will permit the production of beam sections which can be used in guideway structures. The inherent corrosion resistance and low maintenance characteristics of composites should result in lowered maintenance costs over a prolonged life of the structure.
Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.
This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.
The structures of K3Sb5O14 and K2Sb4O11 have been solved by the single-crystal X-ray direct method and the heavy-atom method, respectively. The structure of K3Sb5O14 is orthorhombic, with space group Pbam and cell parameters a = 24.247 (4), b = 7.157 (2), c = 7.334 (2) A, Z = 4. The structure of K2Sb4O11 is monoclinic, with space group C2/m and cell parameters a = 19.473 (4), b = 7.542 (1), c = 7.198 (1) A, beta = 94.82 (2) deg, Z = 4. A full-matrix least-squares refinement gave R = 0.072 and R = 0.067, respectively. In both structures, oxygen atoms form an octahedron around each Sb atom and an irregular polyhedron around each K atom. By sharing corners and edges, the octahedra form a skeleton network having intersecting b-axis and c-axis tunnels. The K(+) ions, which have more than ten oxygen near neighbors, are located in these tunnels. Evidence for K(+)-ion transport within and between tunnels comes from ion exchange of the alkali ions in molten salts and anisotropic temperature factors that are anomalously large in the direction of the tunnels.
The noise induced inside a cylindrical fuselage model by shaker excitation is investigated theoretically and experimentally. The NASTRAN modal-synthesis program is used in the theoretical analysis, and the predictions are compared with experimental measurements in extensive graphs. Good general agreement is obtained, but the need for further refinements to account for acoustic-cavity damping and structural-acoustic interaction is indicated.
We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.
The Supersonic Cruise Airplane Research (SCAR) technology program is reported. Mission requirements, cycle considerations, engine evaluation, and SCAR test results are discussed. It is concluded that improved aerodynamics and structures, greater use of composites, and more-refined active controls are needed in the airframe area for designing an economically viable, environmental acceptable supersonic aircraft.
The structures and transition pressures of Group III-V and II-VI semiconductors and of a pseudobinary system (Ga/x/In/1-x/Sb) have been investigated. Results indicate that GaP, InSb, GaSb, GaAs and possible AlP assume Metallic structures at high pressures; a tetragonal, beta-Sn-like structure is adopted by only InSb and GaSb. The rocksalt phase is preferred in InP, InAs, AlSb, ZnO and ZnS. The model of Van Vechten (1973) gives transition pressures which are in good agreement with measured values, but must be refined to account for the occurrence of the ionic rocksalt structure in some compounds. In addition, discrepancies between the theoretical scaling values for volume changes at the semiconductor-to-metal transitions are observed.
Severe, low-altitude wind shear is a threat to aviation safety. Airborne sensors under development measure the radial component of wind along a line directly in front of an aircraft. In this paper, optimal estimation theory is used to define a detection algorithm to warn of hazardous wind shear from these sensors. To achieve robustness, a wind shear detection algorithm must distinguish threatening wind shear from less hazardous gustiness, despite variations in wind shear structure. This paper presents statistical analysis methods to refine wind shear detection algorithm robustness. Computational methods predict the ability to warn of severe wind shear and avoid false warning. Comparative capability of the detection algorithm as a function of its design parameters is determined, identifying designs that provide robust detection of severe wind shear.