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.
Numerical solution of the MHD wave equations for stability of the cometary sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is obtained by using a two-point boundary value method.
Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.
Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.
A theoretical model for the interaction of a turbulent boundary layer with an oscillating wavy surface over which a fluid is flowing is developed, with an application to wind-driven water waves and to panel flutter in low supersonic flow. A systematic methodology is developed to obtain the surface pressure distribution by considering separately the effects on the perturbed flow of a mean shear velocity profile, viscous stresses, the turbulent Reynolds stresses, compressibility, and three-dimensionality. The inviscid theory is applied to the wind-water wave problem by specializing to traveling-wave disturbances, and the pressure magnitude and phase shift as a function of the wave phase speed are computed for a logarithmic mean velocity profile and compared with inviscid theory and experiment. The results agree with experimental evidence for the stabilization of the panel motion due to the influence of the unsteady boundary layer.
The onset and development of both dynamically and convectively forced boundary layer rolls are studied with linear and nonlinear analyses of a truncated spectral model of shallow Boussinesq flow. Emphasis is given here on the energetics of the dominant roll modes, on the magnitudes of the roll-induced modifications of the initial basic state wind and temperature profiles, and on the sensitivity of the linear stability results to the use of modified profiles as basic states. It is demonstrated that the roll circulations can produce substantial changes to the cross-roll component of the initial wind profile and that significant changes in orientation angle estimates can result from use of a roll-modified profile in the stability analysis. These results demonstrate that roll contributions must be removed from observed background wind profiles before using them to investigate the mechanisms underlying actual secondary flows in the boundary layer. The model is developed quite generally to accept arbitrary basic state wind profiles as dynamic forcing. An Ekman profile is chosen here merely to provide a means for easy comparison with other theoretical boundary layer studies; the ultimate application of the model is to study observed boundary layer profiles. Results of the analytic stability analysis are validated by comparing them with results from a larger linear model. For an appropriate Ekman depth, a complete set of transition curves is given in forcing parameter space for roll modes driven both thermally and dynamically. Preferred orientation angles, horizontal wavelengths and propagation frequencies, as well as energetics and wind profile modifications, are all shown to agree rather well with results from studies on Ekman layers as well as with studies on near-neutral and convective atmospheric boundary layers.
Wendelstein 7-X (W7-X) (Greifswald, Germany) is an advanced stellarator, which uses the modular coil concept to realize a magnetic configuration optimized for fusion-relevant plasma properties. The magnet system of the machine allows a variation of the rotational transform (iota) at the boundary. In the latest W7-X operational phase a dedicated configuration scan has been performed varying the rotational transform between magnetic configurations with iota = 5/4 and iota = 5/5 at the boundary. This paper presents an overview of the experiments and of the main results with respect to confinement and stability. The main observation is an increase of the plasma energy in several intermediate configurations of the scan when the 5/5-islands are close to the plasma boundary but still inside the last-closed-flux-surface. In addition, these configurations showed marked MHD-activity with a crashing behavior related to the 5/5-islands. The corresponding mode amplitude was correlated with the size of the internal 5/5 islands.
This paper reviews the development of generalized unsteady aerodynamic theory and presents a derivation of the generalized Possio integral equation. Numerical calculations resolve questions concerning subsonic indicial lift functions and demonstrate the generation of Kutta waves at high values of reduced frequency, subsonic Mach number, or both. The use of rational function approximations of unsteady aerodynamic loads in aeroelastic stability calculations is reviewed, and a reformulation of the matrix Pade approximation technique is given. Numerical examples of flutter boundary calculations for a wing which is to be flight tested are given. Finally, a simplified aerodynamic model of transonic flow is used to study the stability of an airfoil exposed to supersonic and subsonic flow regions.
The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.
Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.
Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.
The discovery of an extraordinary number of extrasolar planets, characterized by an unexpected variety of sizes, masses and orbits, challenges our understanding of the formation and evolution of the planets in the Solar System and the perception of the Earth as the prototypical habitable world. Many exoplanets appear to be rocky and yet more massive than the Earth, with expected pressures and temperatures of hundreds of gigapascal and thousands of Kelvin in their deep interiors. At these conditions, the properties of bridgmanite and ferropericlase, expected to dominate their mantles, are largely unconstrained, limiting our knowledge of their interior structure. Here we used nano-second X-ray diffraction and dynamic compression to experimentally investigate the atomic structure and density of iron oxide (FeO), one of the end-members of the (Mg,Fe)O ferropericlase solid solution, up to 700 GPa, a pressure exceeding the core–mantle boundary of a 5 Earth masses planet. Our data document the stability of the high-pressure cesium-chloride B2 structure above 300 GPa, well below the pressure required for magnesium oxide (MgO) to adopt the same phase. These observations, complemented by the calculation of the binary MgO–FeO phase diagram, reveal complex stratification and rheology inside large terrestrial exoplanets.
Stability and dissipation in second and third- order fluid approximations of boundary value problems for unsteady simple shear flow
An analytical investigation was made into the stability of the burning configuration of a single coaxial injector surrounded by similar injectors. The stability criteria was based on an average pressure difference along the boundaries of the adjacent stream tubes as calculated using Spaulding's numerical method. The results indicate qualitatively that there is a tendency for the injectors to have different burning configurations. It is believed that the configuration achieved is random, however once the burning configuration is established, it is believed to persist.
Program computes incompressible linear stability characteristics for swept and tapered wings. Amplification rates of boundary-layer disturbances also are calculated. Program is useful in designing tapered, laminar-flow control wings incorporating suction to prevent boundary layer separation. Program is written in FORTRAN IV and Assembler for batch execution on CYBER 70-series.
Self-adaptive solution algorithms are developed for nonlinear heat conduction problems encountered in analyzing materials for use in high temperature or cryogenic conditions. The nonlinear effects are noted to occur due to convection and radiation effects, as well as temperature-dependent properties of the materials. Incremental successive substitution (ISS) and Newton-Raphson (NR) procedures are treated as extrapolation schemes which have solution projections bounded by a hyperline with an externally applied thermal load vector arising from internal heat generation and boundary conditions. Closed constraints are formulated which improve the efficiency and stability of the procedures by employing closed ellipsoidal surfaces to control the size of successive iterations. Governing equations are defined for nonlinear finite element models, and comparisons are made of results using the the new method and the ISS and NR schemes for epoxy, PVC, and CuGe.
A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.
A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.
In the present analytical procedure for the prediction of helicopter rotor noise generation due to the ingestion of atmospheric turbulence, different models for turbulence fluid mechanics and the ingestion process are combined. The mean flow and turbulence statistics associated with the atmospheric boundary layer are modeled with attention to the effects of atmospheric stability length, windspeed, and altitude. The turbulence field can be modeled as isotropic, locally stationary, and homogeneous. For large mean flow contraction ratios, accurate predictions of turbulence vorticity components at the rotor face requires the incorporation of the differential drift of fluid particles on adjacent streamlines.