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.
The behavior of gas dynamic flows which are perturbations of a uniform stream in terms of information transfer across artificial (computational) boundaries remote from the source of disturbance are discussed. A set of boundary conditions is derived involving vorticity, entropy, and pressure-velocity relationships derived from bicharacteristic equations.
Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.
Superplasticity, the capacity of materials to sustain extraordinary tensile elongations at elevated temperatures, underpins a range of advanced metal-forming technologies. Conventionally, it is achieved in fine-grained alloys where deformation is dominated by grain-boundary sliding, accommodated by diffusion and dislocation activity. The advent of medium- and high-entropy alloys (M/HEAs), with their high chemical complexity and unconventional phase stability, offers new pathways to superplastic behavior beyond traditional alloy systems. Although investigated only recently, several M/HEAs already exhibit elongations that rival or exceed those of classical superplastic materials, particularly when ultrafine or metastable microstructures are engineered. Here, we review progress in understanding superplastic deformation in M/HEAs, emphasizing the interplay among composition, initial microstructure, thermomechanical processing, and microstructural evolution during high-temperature deformation. We discuss approaches to generating the fine-grained structures necessary for grain-boundary sliding, including severe plastic deformation and tailored heat treatments. We further highlight dynamic phenomena such as phase transformations, evolving grain-boundary chemistry, and deformation-induced grain refinement that can enhance plasticity in these systems. These mechanisms often shift the balance of deformation processes, enabling large elongations even outside classical criteria. Finally, we outline key challenges for application, including cost, scalability, recyclability, and microstructural stability.
In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.
The flutter characteristics of a generic arrow-wing supersonic transport configuration are studied. The wing configuration has a 3 percent biconvex airfoil and a leading-edge sweep of 73 deg out to a cranked tip with a 60 deg leading-edge sweep. The ground vibration tests and flutter test procedure are described. The effects of flutter on engine nacelles, fuel loading, wing-mounted vertical fin, wing angle-of-attack, and wing tip mass and stiffness distributions are analyzed. The data reveal that engine nacelles reduce the transonic flutter dynamic pressure by 25-30 percent; fuel loadings decrease dynamic pressures by 25 percent; 4-6 deg wing angles-of-attack cause steep transonic boundaries; and 5-10 percent changes in flutter dynamic pressures are the result of the wing-mounted vertical fin and wing-tip mass and stiffness distributions.
Mason and Weaver developed equations to describe small particles settling under the influence of gravity and Brownian motion, including the limiting case for an infinitely deep suspension. We encountered this common convection–diffusion equation and no-flux boundary conditions in a model for dynamics of adsorbed polymers in dead end pores of a depolymerization catalyst. Close examination reveals that the Mason–Weaver solution is not correct for the infinite domain with a non-uniform initial condition. In this paper, we obtain the time dependent Green’s function for the no flux boundary condition and also for a more general reactive boundary condition. We demonstrate how the results provide solutions, via superposition, which provide solutions for several boundary conditions and all initial conditions.
A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.
The dynamic behavior of flexible rotating beams continues to receive considerable research attention as it constitutes a fundamental problem in applied mechanics. Further, beams comprise parts of many rotating structures of engineering significance. A topic of particular interest at the present time involves the development of techniques for obtaining the behavior in both space and time of a rotor acted upon by a simple airload loading. Most current work on problems of this type use solution techniques based on normal modes. It is certainly true that normal modes cannot be disregarded, as knowledge of natural blade frequencies is always important. However, the present work has considered a computational structural mechanics (CSM) approach to rotor blade dynamics problems in which the physical properties of the rotor blade provide input for a direct numerical solution of the relevant boundary-and-initial-value problem. Analysis of the dynamics of a given rotor system may require solution of the governing equations over a long time interval corresponding to many revolutions of the loaded flexible blade. For this reason, most of the common techniques in computational mechanics, which treat the space-time behavior concurrently, cannot be applied to the rotor dynamics problem without a large expenditure of computational resources. By contrast, the integrating matrix technique of computational mechanics has the ability to consistently incorporate boundary conditions and 'remove' dependence on a space variable. For problems involving both space and time, this feature of the integrating matrix approach thus can generate a 'splitting' which forms the basis of an efficient CSM method for numerical solution of rotor dynamics problems.
The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.
OGO-5 observations show not only that the average position of the magnetopause boundary moves in response to changes in dynamic pressure, but also that it moves in response to changes in the north-south component of the interplanetary field. Further, the boundary often oscillates about its average position as waves propagate along the boundary away from the nose region. The variation of the magnetic field through the boundary at times can be a simple rotation with no change in magnitude, and at other times can resemble a simple tangential discontinuity. However, it often displays complex patterns such as field enhancements on the magnetospheric side of the boundary or apparently uncorrelated field strength and direction changes. One particularly simple boundary crossing has been studied in detail with both positive ion data and magnetic field data. In this case, the electron and ion currents could be separately deduced and the structure agreed with that expected for a Chapman-Ferraro boundary with almost complete neutralization.
A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.
Loci-CHEM to be applied to flow velocities where surface radiation due to heating from compression and friction becomes significant. The module adds a radiation equilibrium boundary condition to the computational fluid dynamics (CFD) code to produce accurate results. The module expanded the upper limit for accurate CFD solutions of Loci-CHEM from Mach 4 to Mach 10 based on Space Shuttle Orbiter Re-Entry trajectories. Loci-CHEM already has a very promising architecture and performance, but absence of radiation equilibrium boundary condition limited the application of Loci-CHEM to below Mach 4. The immediate advantage of the add-on module is that it allows Loci-CHEM to work with supersonic flows up to Mach 10. This transformed Loci-CHEM from a rocket engine- heritage CFD code with general subsonic and low-supersonic applications, to an aeroheating code with hypersonic applications. The follow-on advantage of the module is that it is a building block for additional add-on modules that will solve for the heating generated at Mach numbers higher than 10.
In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, but do not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the frequency range of interest of the test article. Dynamic coupling between the test article and test fixture complicates the model updating process because significant effort is required to model the test fixture and boundary conditions in addition to the test article. If there were a way to correct the modal results for fixture coupling, then setups used for other structural testing could be adequate for modal testing. In the case described in this paper, a partial static loads testing setup was used, which allowed significant schedule and cost savings by eliminating a unique setup for a modal test. A fixed base correction technique was investigated during modal testing of a flexible wing cantilevered from part of a static test fixture. The technique was successfully used to measure the wing modes de-coupled from the dynamically active test fixture. The technique is promising for future aircraft applications, but more research is needed.
The life of the very small, whether in something as complicated as a human cell or as simple as a drop of water, is of fundamental scientific interest: By knowing how a tiny amount of material reacts to changes in its environment, scientists maybe able to answer questions about how a bulk of material would react to comparable changes. NASA is in the forefront of computational research into a broad range of basic scientific questions about fluid dynamics and the nature of liquid boundary instability. For example, one important issue for the space program is how drops of water and other materials will behave in the low-gravity environment of space and how the low gravity will affect the transport and containment of these materials. Accurate prediction of this behavior is among the aims of a set of molecular dynamics experiments carried out on the NCCSs Cray supercomputers. In conventional computational studies of materials, matter is treated as continuous - a macroscopic whole without regard to its molecular parts - and the behavior patterns of the matter in various physical environments are studied using well-established differential equations and mathematical parameters based on physical properties such as compressibility density, heat capacity, and vapor pressure of the bulk material.
The effects of dynamic stall on airfoils oscillating in pitch were investigated by experimentally determining the viscous and inviscid characteristics of the airflow on the NACA 0012 airfoil and on several leading-edge modifications. The test parameters included a wide range of frequencies, Reynolds numbers, and amplitudes-of-oscillation. Three distinct types of separation development were observed within the boundary layer, each leading to classical dynamic stall. The NACA 0012 airfoil is shown to stall by the mechanism of abrupt turbulent leading-edge separation. A detailed step-by-step analysis of the events leading to dynamic stall, and of the results of the stall process, is presented for each of these three types of stall. Techniques for flow analysis in the dynamic stall environment are discussed. A method is presented that reduces most of the oscillating airfoil normal force and pitching-moment data to a single curve, independent of frequency or Reynolds number.
A major challenge in solar and heliospheric physics is understanding the origin and nature of the so-called slow solar wind. The Sun's atmosphere is divided into magnetically open regions, known as coronal holes, where the plasma streams out freely and fills the solar system, and closed regions, where the plasma is confined to coronal loops. The boundary between these regions extends outward as the heliospheric current sheet (HCS). Measurements of plasma composition strongly imply that much of the slow wind consists of plasma from the closed corona that escapes onto open field lines, presumably by field-line opening or by interchange reconnection. Both of these processes are expected to release closed-field plasma into the solar wind within and immediately adjacent to the HCS. Mysteriously, however, slow wind with closed-field plasma composition is often observed in situ far from the HCS. We use high-resolution, three-dimensional, magnetohydrodynamic simulations to calculate the dynamics of a coronal hole with a geometry that includes a narrow corridor flanked by closed field and is driven by supergranule-like flows at the coronal-hole boundary. These dynamics produce giant arcs of closed-field plasma that originate at the open-closed boundary in the corona, but extend far from the HCS and span tens of degrees in latitude and longitude at Earth. We conclude that such structures can account for the long-puzzling slow-wind observations.
The importance of unsteady aerodynamic effects due to rotor blade motion is demonstrated for a helicopter in high-speed forward flight. These unsteady effects are modeled as surface boundary conditions to a three-dimensional, unsteady Computational Fluid Dynamics (CFD) code called the Full-Potential Rotor Code (FPR). These boundary conditions cause significant changes in the computed lift and pitching moment at the front and rear of the rotor disk. Airloads from the modified FPR code are then iteratively coupled with the helicopter comprehensive code, CAMRAD/JA. Computed airloads show good agreement with flight-test data when lift values from the FPR code are used in the coupled calculation. However, the computed airloads from CAMRAD/JA along also show good agreement with the experimental data. Thus for this case one cannot demonstrate a significant improvement in computed airloads with the hybrid coupled scheme. The addition of the pitching moment values from the FPR code into the CAMRAD/JA calculation slows down the overall iterative convergence and does not yield any improvement in the final results.
Abstract In polycrystalline materials, grain boundaries are known to be a critical microstructural component controlling material’s mechanical properties, and their characters such as misorientation and crystallographic boundary planes would also influence the dislocation dynamics. Nevertheless, many of generally used mechanistic models for deformation twin nucleation in fcc metal do not take considerable care of the role of grain boundary characters. Here, we experimentally reveal that deformation twin nucleation occurs at an annealing twin (Σ3{111}) boundary in a high-Mn austenitic steel when dislocation pile-up at Σ3{111} boundary produced a local stress exceeding the twining stress, while no obvious local stress concentration was required at relatively high-energy grain boundaries such as Σ21 or Σ31. A periodic contrast reversal associated with a sequential stacking faults emission from Σ3{111} boundary was observed by in-situ transmission electron microscopy (TEM) deformation experiments, proving the successive layer-by-layer stacking fault emission was the deformation twin nucleation mechanism, different from the previously reported observations in the high-Mn steels. Since this is also true for the observed high Σ-value boundaries in this study, our observation demonstrates the practical importance of taking grain boundary characters into account to understand the deformation twin nucleation mechanism besides well-known factors such as stacking fault energy and grain size.