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.
Postprocessing algorithms devised to facilitate vibrational-mode analyses of dynamics of complicated structures. Yields inertia matrices and elastic/rigid-coupling matrices. Such analyses important in simulation and control in active suppression of vibrations in large building or in precise aiming of large antenna.
This paper describes an aeroelastic analysis program for turbomachines. Unsteady Navier-Stokes equations are solved on dynamically deforming, body fitted, grid to obtain the aeroelastic characteristics. Blade structural response is modeled using a modal representation of the blade and the work-per-cycle method is used to evaluate the stability characteristics. Nonzero interblade phase angle is modeled using phase-lagged boundary conditions. Results obtained showed good correlation with existing experimental, analytical, and numerical results. Numerical analysis also showed that given the computational resources available today, engineering solutions with good accuracy are possible using higher fidelity analyses.
Excessive vibration remains one one of the most difficult problems that faces the helicopter industry today, affecting all production helicopters at some phase of their development. Vibrations in rotating structures may arise from external periodic dynamic airloads whose frequencies are are close to the natural frequencies of the rotating system itself. The goal for the structures engineer would thus be to design a structure as free from resonance effects as possible. In the case of a helicopter rotor blade these dynamic loads are a consequence of asymmetric airload distribution on the rotor blade in forward flight, leading to a rich collection of higher harmonic airloads that force rotor and airframe response. Accurate prediction of the dynamic characteristics of a helicopter rotor blade will provide the opportunity to affect in a positive manner noise intensity, vibration level, durability, reliability and operating costs by reducing objectionable frequencies or moving them to a different frequency range and thus providing us with a lower vibration rotor. In fact, the dynamic characteristics tend to define the operating limits of a rotorcraft. As computing power has increased greatly over the last decade, researchers and engineers have turned to analyzing the vibrational characteristics of aerospace structures at the design and development stage of the production of an aircraft. Modern rotor blade construction methods lead to products with low mass and low inherent damping so careful design and analysis is required to avoid resonance and an undesirable dynamic performance. In addition, accurate modal analysis is necessary for several current approaches in elastic system identification and active control.
Predictive methods/computer codes for the computational simulation of acoustic fatigue resistance of hot composite structures subjected to acoustic excitation emanating from an adjacent vibrating component are discussed. Select codes developed over the past two decades at the NASA Lewis Research Center are used. The codes include computation of acoustic noise generated from a vibrating component, degradation in material properties of a composite laminate at use temperature, dynamic response of acoustically excited hot multilayered composite structure, degradation in the first ply strength of the excited structure due to acoustic loading, and acoustic fatigue resistance of the excited structure, including the propulsion environment. Effects of the laminate lay-up and environment on the acoustic fatigue life are evaluated. The results show that, by keeping the angled plies on the outer surface of the laminate, a substantial increase in the acoustic fatigue life is obtained. The effect of environment (temperature and moisture) is to relieve the residual stresses leading to an increase in the acoustic fatigue life of the excited panel.
Strongly coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute nonequal-time correlation functions and perform entanglement tomography of nonequilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently predicted dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent nonequilibrium dynamics in three ways: (i) overlap echos signaling dynamical transitions, (ii) nonequal-time correlation functions with an underlying topological nature, and (iii) the entanglement structure of nonequilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.
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.
Remdesivir was the first antiviral drug that received emergency use authorization from the United States Food and Drug Administration and is now formally approved to treat COVID-19. Remdesivir is a nucleotide analogue that targets the RNA-dependent RNA polymerase (RdRp) of coronaviruses, including SARS-CoV-2. The solution of multiple RdRp structures has been one of the main axes of research in the race against the SARS-CoV-2 virus. Several hypotheses of the mechanism of inhibition of RdRp by remdesivir have been proposed, although open questions remain. This work uses molecular dynamics simulations to explore the impact of remdesivir and two analogues as incoming nucleotides and of up to four incorporations of remdesivir along the primer strand on RdRp. The simulation results suggest that the overall structure and the dynamical behavior of RdRp are destabilized by remdesivir and the two analogues in the incoming position. The incorporation of remdesivir along the primer strand impacts specific non-bonded interactions between the nascent RNA and the polymerase subunit, as well as the overall dynamical networks on RdRp. The strongest impact on the structure and dynamics are observed after three incorporations, when remdesivir is located at position -A3, in agreement with previously reported experimental and computational results. Our results provide atomic-level details of the role played by remdesivir on the disruption of RNA synthesis by RdRp and the main drivers of these disruptions.
The rotary bacterial flagellar motor is remarkable in biochemistry for its highly synchronized operation and amplification during switching of rotation sense. The motor is part of the flagellar basal body, a complex multi-protein assembly. Sensory and energy transduction depends on a core of six proteins that are adapted in different species to adjust torque and produce diverse switches. Motor response to chemotactic and environmental stimuli is driven by interactions of the core with small signal proteins. The initial protein interactions are propagated across a multi-subunit cytoplasmic ring to switch torque. Torque reversal triggers structural transitions in the flagellar filament to change motile behavior. Subtle variations in the core components invert or block switch operation. The mechanics of the flagellar switch have been studied with multiple approaches, from protein dynamics to single molecule and cell biophysics. The architecture, driven by recent advances in electron cryo-microscopy, is available for several species. Computational methods have correlated structure with genetic and biochemical databases. The design principles underlying the basis of switch ultra-sensitivity and its dependence on motor torque remain elusive, but tantalizing clues have emerged. This review aims to consolidate recent knowledge into a unified platform that can inspire new research strategies.
Topological defects in ferroelectric materials have attracted much attention due to the emergence of conductive, ferroic, and magnetic functionalities. However, many topological configurations dynamically evolve during the switching processes, making them a challenge to characterize via traditional techniques. Here, we implement an automated experimentation approach for the exploration of functional properties in BiFeO 3 thin films. Specifically, we visualize the ferroelectric domain structures via single frequency piezoresponse force microscopy and implement a computer vision-based algorithm to discover features of interest at which spectroscopic measurements are taken. Subsequently, we employ dimensionality reduction techniques to reveal characteristic polarization behaviors at these features. This approach can be extended to other spectroscopies and modalities to probe only specific features of interest, ultimately enabling dynamical processes in ferroelectrics to be studied.
Tracking the multifarious ultrafast electronic and structural changes occurring in a molecule during a photochemical transformation is a challenging endeavor that benefits from recent experimental and computational progress in time-resolved techniques. Measurements of valence electronic states, which provide a global picture of the bonding structure of the molecule, and core electronic states, which provide insight into the local environment, traditionally require different approaches and are often studied separately. Here, we demonstrate that X-ray pulses from a seeded free-electron laser (FEL) enable the measurement of high-resolution, time-resolved X-ray photoelectron spectra (XPS) that capture weak satellite states resulting from shake-down processes in a valence-excited molecule. This approach effectively combines the advantages of both valence- and core-state investigations. We applied this method to investigate photoexcited CS2 molecules, where the role of internal conversion (IC) and intersystem crossing (ISC) in determining the predissociation dynamics is controversial. We present XPS spectra from photoexcited CS2, obtained at the FERMI FEL. High-resolution measurements, compared to the corresponding spectra obtained from accurate multireference quantum chemical calculations, reveal that shake-down satellite channels are highly sensitive to both valence electronic and geometric changes. Previous studies of the predissociation dynamics have led to uncertain assignments of the branching between singlet and triplet excited states. We derive a propensity rule that demonstrates the spin-selectivity of the shake-downs. This selectivity allows us to unequivocally assign contributions from the bright and dark singlet excited states, with populations tracked along the predissociation dynamic pathway.
This paper presents a modal method for the analysis of controlled structural systems that retains the uncoupled nature of the classical transient response analysis of a structure subjected to a prescribed time-varying load. The control force is expanded as a Taylor series that remains on the right side of the equations, and it does not lead to a computational approach that requires coupling between modes on the left side. Retaining a sufficient number of terms in the series produces a solution to the modal equations that is accurate to machine precision. The approach is particularly attractive for large problems in which standard matrix exponential methods become computationally prohibitive. Numerical results are presented to show the accuracy and efficiency of the proposed approach for dynamic feedback compensation of a truss structure with local member modes in the controller bandwidth.
The main objective of the structures group is to provide quality aerospace research with the Center for Aerospace Research - A NASA Center for Excellence at North Carolina Agricultural and Technical State University. The group includes dedicated faculty and students who have a proven record in the area of structures, in particular space structures. The participating faculty developed accurate mathematical models and effective computational algorithms to characterize the flexibility parameters of joint dominated beam-truss structures. Both experimental and theoretical modelling has been applied to the dynamic mode shapes and mode frequencies for a large truss system. During the past few months, the above procedures has been applied to the hypersonic transport plane model. The plane structure has been modeled as a lumped mass system by Doctor Abu-Saba while Doctor Shen applied the transfer matrix method with a piecewise continuous Timoshenko tapered beam model. Results from both procedures compare favorably with those obtained using the finite element method. These two methods are more compact and require less computer time than the finite element method. The group intends to perform experiments on structural systems including the hypersonic plane model to verify the results from the theoretical models.
The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.
Fast 3D X-ray imaging has proved to provide crucial insights into multi-phase dynamic processes in various geomaterials. In particular, many in situ imaging experiments have been made to study gas-hydrate formation in porous sandy samples. Such imaging is challenging for the methane gas-hydrate formation in coal samples because of the coal micro-porosity structure and lower X-ray contrast. Here we present results of the first dynamic in situ micro-computed tomography experiment of methane-hydrate formation in coal samples. Synchrotron phase-contrast tomography techniques allowed to achieve necessary contrast levels to separate all the materials of interest (gas, water, coal, and gas hydrate) in reconstructed images with high spatial and temporal resolution. The imaging results are compared to the ones from a similar tomographic experiment with sand samples. Methane-hydrate formation is accompanied by the water movements caused by cryogenic water suction that happens in sequences of short fast movements with longer equilibrium states in between. Only one type of the hydrate formation was observed in coal (growth as shells on grain boundaries) as opposed to three types in sand (shells, growth into gas pockets, and inside water volumes). In particular, this leads to a slower hydrateformation speed in coal. For the coal sample, we also observed water extraction from grains, and interpreted it as competitive sorption of methane. Lastly, we visualized the dynamic behavior of this water extraction via microchannels inside the coal grains, and performed nano-tomography imaging of these channels for a better understanding of this phenomenon.
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. TDK has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 716142). VRR has been supported by the Swiss National Science Foundation in the form of Ambizione grant No. PZ00P2 174227 and RZK by the Natural Sciences and Engineering Research Council of Canada (NSERC) through Discovery Grants (RGPIN-2016-0505). GKS and CJM are supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. UK based work was funded under the embedded CSE programme of the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), grants eCSE03-011, eCSE06-6, eCSE08-9, eCSE13-17 and the EPSRC (EP/P022235/1) grant “Surface and Interface Toolkit for the Materials Chemistry Community". Computational resources were provided by the Swiss National Supercomputing Centre (CSCS) and Compute Canada. The generous allocation of computing time on the FPGA-based supercomputer “Noctua" at PC2 is kindly acknowledged.
Kolmogorov's theory of turbulence assumes that the small-scale turbulent structures in the energy cascade are universal and are determined by the energy dissipation rate and the kinematic viscosity alone. However, thermal fluctuations, absent from the continuum description, terminate the energy cascade near the Kolmogorov length scale. Here, we propose a simple superposition model to account for the effects of thermal fluctuations on small-scale turbulence statistics. Furthermore, for compressible Taylor–Green vortex flow, we demonstrate that the superposition model in conjunction with data from direct numerical simulation of the Navier–Stokes equations yields spectra and structure functions that agree with the corresponding quantities computed from the direct simulation Monte Carlo method of molecular gas dynamics, verifying the importance of thermal fluctuations in the dissipation range.
The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.
Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.