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.
Relaxations from anelastic and dielectric measurements of fused silica and soda-silica glasses at low temperatures ascribed to vibration modes exceeding Debye spectrum
A practical procedure for optimum design of aerodynamic shapes is demonstrated. The proposed procedure uses an optimization program based on the method of feasible directions coupled with an analysis program that uses a relaxation solution of the inviscid, transonic, small-disturbance equations. Results are presented for low-drag, nonlifting transonic airfoils. Extension of the method to lifting airfoils, other speed regimes, and to three dimensions if feasible.
Highly flare-productive newly emerging active regions are characterized by numerous small low-lying loops which frequently show a chaotic pattern. Flare activity in such a region subsides as the chaotic loop structures relax and expand into a bipolar configuration. The transition zone in such an active region is highly unstable, as shown by broadened and shifted nonthermal line profiles of medium-ionized elements like Si III, Si IV, and C IV. These transition-zone instabilities, which occur as isolated events in active regions of low flare productivity, are often observed prior to flares. Transition-zone instabilities can be traced to the footpoints of active loops and seem to be accompanied by heating of the loop. The loops vary in size and show differing degrees of activity, with the brightest and most compact ones seemingly being in a pre-flare state which results in the catastrophic energy release along the loop during a flare.
With CMC components poised to complete flight certification in turbine engines on commercial aircraft within the near future, there are many efforts within the aerospace community to model the mechanical and environmental degradation of CMCs. Direct observations of damage evolution are needed to support these modeling efforts and provide quantitative measures of damage parameters used in the various models. This study was performed to characterize the damage evolution during tensile loading of a melt infiltrated (MI) silicon carbide reinforced silicon carbide (SiC/SiC) composite. A SiC/SiC tensile coupon was loaded to a maximum global stress of 30 ksi in a tensile fixture within an SEM while observations were made at 5 ksi increments. Both traditional image analysis and DIC (digital image correlation) were used to quantify damage evolution. With the DIC analysis, microscale damage was observed at the fiber-matrix interfaces at stresses as low as 5 ksi. First matrix cracking took place between 20 and 25 ksi, accompanied by an observable relaxation in strain near matrix cracks. Matrix crack opening measurements at the maximum load ranged from 200 nm to 1.5 m. Crack opening along the fiber-matrix interface was also characterized as a function of load and angular position relative to the loading axis. This characterization was funded by NASA GRC and was performed to support NASA GRC modeling of SiC/SiC environmental degradation
Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.
State variable techniques are utilized to generate the RMS vertical load factors for the B-52H and B-1 bombers at low level, mission critical, cruise conditions. A ride quality index is proposed to provide meaningful comparisons between different controls or conditions. Ride quality is shown to be relatively invariant under various popular control laws. Handling quality variations are shown to be major contributors to ride quality variations on both vehicles. Relaxed static stability is artificially implemented on the study vehicles to investigate its effects on ride quality. The B-52H ride quality is generally degraded when handling characteristics are automatically restored by a feedback control to the original values from relaxed stability conditions. The B-1 airplane shows little ride quality sensitivity to the same analysis due to the small rigid body contribution to load factors at the flight condition investigated.
Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.
A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments relies upon the global momentum conservation of the fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. Numerical examples illustrate the method's application to predicting bulk fluid motion including lateral propellant slosh in low-g conditions.
In a recent paper, Kuchner, Crepp, and Ge describe new image-plane coronagraph mask designs that reject to eighth order the leakage of starlight caused by image motion at the mask, resulting in a substantial relaxation of image centroiding requirements compared to previous fourth-order and second-order masks. They also suggest that the new masks are effective at rejecting leakage caused by low-order aberrations (e.g., focus, coma, and astigmatism). In this paper, we derive the sensitivity of eighth-order masks to aberrations of any order and provide simulations of coronagraph behavior in the presence of optical aberrations.We find that the masks leak light as the fourth power of focus, astigmatism, coma, and trefoil. This has tremendous performance advantages for the Terrestrial Planet Finder Coronagraph.
The low speed, high angle of attack flight mechanics of an advanced, canard-configured, supersonic tactical aircraft designed with moderate longitudinal relaxed static stability (Static Margin, SM = 16% C sub W at M = 0.4) was investigated. Control laws were developed for the longitudinal axis (""G'' or maneuver and angle of attack command systems) and for the lateral/directional axes. The performance of these control laws was examined in engineering simulation. A canard deflection/rate requirement study was performed as part of the ""G'' command law evaluation at low angles of attack. Simulated coupled maneuvers revealed the need for command limiters in all three aircraft axes to prevent departure from controlled flight. When modified with command/maneuver limiters, the control laws were shown to be adequate to prevent aircraft departure during aggressive air combat maneuvering.
Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.
Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions
Structural distortions in low-dimensional π-conjugated systems profoundly influence their electronic properties, but the control of such behavior in laterally extended systems remains challenging. Here, in this study, we demonstrate that a one-dimensional conjugated polymer─ poly -(difluorenoheptalene-ethynylene) (PDFHE)─undergoes a pronounced out-of-plane backbone distortion, equivalent to a spontaneous symmetry breaking (SSB) of its mirror symmetry. We synthesized PDFHE on noble metal surfaces and characterized its structure and electronic states using low-temperature scanning tunneling microscopy (STM). Rather than adopting a planar, high-symmetry conformation, PDFHE relaxes into nonplanar isomers stabilized by a pseudo–Jahn–Teller (PJT) distortion having mirror-odd out-of-plane character. The distortion lowers the total energy and increases the band gap, providing a concise rationale for the observed symmetry breaking. Density functional theory calculations corroborate these findings, providing a microscopic explanation for the SSB. Our results show that even in mechanically robust extended π-systems, subtle electron–lattice coupling can spontaneously drive significant structural rearrangements, even in mechanically robust extended π-systems.
Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.
Explore the source record for details and available documents.
In this article we study microwave absorption in superconductors in the presence of a vortex lattice. Here, we show that in addition to the conventional absorption mechanism associated with the vortex core motion, there is another mechanism of microwave absorption, which is caused by the time dependence of the quasiparticle density of states outside the vortex cores. This mechanism exists even in the absence of vortex motion. It provides the dominant contribution to microwave absorption at sufficiently small magnetic fields. At low frequencies, the dissipative part of the microwave conductivity σ(ω) is proportional to the inelastic relaxation time τ in , which is typically much larger than the elastic relaxation time τ el . At high frequencies σ(ω) is proportional to the quasiparticle diffusion time across the intervortex distance τ D , which is still larger than τ el .
Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.
The three-dimensional vortical viscous flow past a sideslipping delta wing at a high angle of attack is numerically investigated. A computational method for calculating low-speed viscous flowfields is developed. The time-iterative method uses an implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. An algebraic turbulence model is implemented to account for the eddy viscosity. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. Computed results are compared with experimental data.