Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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.

At least 739 records · Page 41

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal conductivity and its relation to atomic structure for symmetrical tilt grain boundaries in silicon

We perform a systematic study of thermal resistance and conductance of tilt grain boundaries (GBs) in Si using classical molecular dynamics. The GBs studied are naturally divided into three groups according to the structural units forming the GB core. In this work we find that, within each group, the GB thermal conductivity strongly correlates with the excess GB energy. All three groups predict nearly the same GB conductivity extrapolated to the high-energy limit. This limiting value is close to the thermal conductivity of amorphous Si, suggesting similar heat transport mechanisms. While the lattice thermal conductivity decreases with temperature, the GB conductivity slightly increases. However, at high temperatures it turns over and starts decreasing if the GB structure undergoes a premelting transformation. Analysis of vibrational spectra of GBs resolved along different directions sheds light on the mechanisms of their thermal resistance. The existence of alternating tensile and compressive atomic environments in the GB core gives rise to localized vibrational modes, frequency gaps creating acoustic mismatch with lattice phonons, and anharmonic vibrations of loosely bound atoms residing in open atomic environments.

36 MATERIALS SCIENCE↗

Towards exascale for wind energy simulations

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.

17 WIND ENERGY↗

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗

Calculated global conditions for the unsteady Euler's flow equations

This paper considers the numerical solution of a fluid dynamical system with time-dependent oscillating boundaries. The problem detailed, a compressible and inviscid flow in a one-dimensional pipe with vibrating ends, is examined to develop techniques needed to solve the more complex unsteady flows in turbomachine cascades. The problem is shown to be mathematically ill-posed but by complementing the governing equations with energy conservation, a unique solution is obtained. An artificial damping term is introduced into the governing equations to improve the performance of the time-marching integration.

Kuzo, D. M.↗

The disturbance flow field produced by an evolving vortex

The flow field of a vortex in a viscous shear flow is found by constructing a uniformly valid asymptotic expansion consisting of an inner solution field represented, to lowest order, by a two dimensional, nonliner, inviscid Stuart vortex and an outer solution field represented, to lowest order, by either a two dimensional parallel or self similar viscous flow. The technique involves scaling both the transverse and streamwise coordinates in the vicinity of the vortex as well as allowing for a slow variation of the outer viscous flow. Criteria are established for both the size of the vortical structure and proximity to the boundary surfaces. The composite solution is a consistent mathematical picture of the flow field at a fixed streamwise location as the vortical structure evolves past this point. Such a formulation is also useful in the specification of boundary or initial conditions in numerical fluid dynamic calculations, where an inconsistent setting of these conditions leads to spurious results for rather long computation times.

Gatski, T. B.↗

Fine structure and dynamics of the Harang discontinuity - An event study

Two rockets were launched simultaneously from the Andoya Rocket Range in northern Norway into and through a Harang discontinuity region during a small isolated substorm. Each was equipped with an array of instruments for monitoring ionospheric electric fields and auroral particles. Different sized rockets allowed both to traverse a range of 300 km while introducing altitude separations of up to 440 km, and separations between the feet of the magnetic field lines intersecting the payloads of up to 100 km. The data sets, coupled with multi-station ground-based observations, provide information on the structure and dynamics of the Harang discontinuity region. Two boundaries were encountered. The first electric field reversal was observed simultaneously by both payloads in a region of weak precipitation, while the second reversal was associated with an intense auroral band. Since an unambiguous interpretation of these data is not possible, two alternative scenarios are presented. In one the discontinuity becomes realigned during the flight following a decay in the activity. The second involves a triple cell convection system, possibly the result of an eddy in the flow.

Egeland, A.↗

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗

A new method to predict unsteady aeroelastic behavior

A new method for predicting subsonic flutter and static deflections, including divergence, has been developed. The present method accounts for aspect ratio and, in the case of flutter, static deflections. The angle of attack is limited only by the occurrence of stall or vortex bursting near the wing. The innovation in the present method is to integrate simultaneously and interactively the equations of motion of the structure and the flowfield. The present approach employs an iterative scheme based on the predictor-corrector method. The general unsteady vortex-lattice method (UVLM) is used to predict the aerodynamic loads. Because the UVLM predicts the wakes as part of the solution, the history of the motion is taken into account; hysteresis is predicted. The deflection (for both bending and torsion) is expressed as an expansion in terms of the free-vibration modes. The time-dependent coefficients in these expansions serve as the generalized coordinates. Numerical examples illustrating the calculation of static deflections and transient dynamic responses above and below the flutter boundary are included.

Strganac, Thomas W.↗

An assessment of wind tunnel test data on flexible thermal protection materials and results of new fatigue tests of threads

Advanced Flexible Reusable Surface Insulation (AFRSI) was developed as a replacement for the low-temperature (white) tiles on the Space Shuttle. The first use of the AFRSI for an Orbiter flight was on the OMS POD of Orbiter (OV-099) for STS-6. Post flight examination after STS-6 showed that damage had occurred to the AFRSI during flight. The failure anomaly between previous wind-tunnel tests and STS-6 prompted a series of additional wind tunnel tests to gain an insight as to the cause of the failure. An assessment of all the past STS-6 wind tunnel tests pointed out the sensitivity of the test results to scaling of dynamic loads due to the difference of boundary layer thickness, and the material properties as a result of exposure to heating. The thread component of the AFRSI was exposed to fatigue testing using an apparatus that applied pulsating aerodynamic loads on the threads similar to the loads caused by an oscillating shock. Comparison of the mean values of the number-of-cycles to failure showed that the history of the thread was the major factor in its performance. The thread and the wind tunnel data suggests a mechanism of failure for the AFRSI.

Coe, Charles F.↗

The multiple-function multi-input/multi-output digital controller system for the AFW wind-tunnel model

A real time multiple-function digital controller system was developed for the Active Flexible Wing (AFW) Program. The digital controller system (DCS) allowed simultaneous execution of two control laws: flutter suppression and either roll trim or a rolling maneuver load control. The DCS operated within, but independently of, a slower host operating system environment, at regulated speeds up to 200 Hz. It also coordinated the acquisition, storage, and transfer of data for near real time controller performance evaluation and both open- and closed-loop plant estimation. It synchronized the operation of four different processing units, allowing flexibility in the number, form, functionality, and order of control laws, and variability in the selection of the sensors and actuators employed. Most importantly, the DCS allowed for the successful demonstration of active flutter suppression to conditions approximately 26 percent (in dynamic pressure) above the open-loop boundary in cases when the model was fixed in roll and up to 23 percent when it was free to roll. Aggressive roll maneuvers with load control were achieved above the flutter boundary. The purpose here is to present the development, validation, and wind tunnel testing of this multiple-function digital controller system.

Hoadley, Sherwood T.↗

The multiple-function multi-input/multi-output digital controller system for the AFW wind tunnel model

A real-time multiple-function digital controller system was developed for the Active Flexible Wing (AFW) Program. The digital controller system (DCS) allowed simultaneous execution of two control laws: flutter suppression and either roll trim or a rolling maneuver load control. The DCS operated within, but independently of, a slower host operating system environment, at regulated speeds up to 200 Hz. It also coordinated the acquisition, storage, and transfer of data for near real-time controller performance evaluation and both open- and closed-loop plant estimation. It synchronized the operation of four different processing units, allowing flexibility in the number, form, functionality, and order of control laws, and variability in selection of sensors and actuators employed. Most importantly, the DCS allowed for the successful demonstration of active flutter suppression to conditions approximately 26 percent (in dynamic pressure) above the open-loop boundary in cases when the model was fixed in roll and up to 23 percent when it was free to roll. Aggressive roll maneuvers with load control were achieved above the flutter boundary.

Hoadley, Sherwood T.↗

Subharmonic Route to Boundary-Layer Transition - Critical Layer Nonlinearity

The linear and nonlinear dynamics of a triad of initially linear stability waves comprising a single plane wave at fundamental frequency and two symmetric oblique waves with half the frequency and streamwise wave number of the plane wave are presented. Analysis is performed for the initial nonlinear development of the waves where the order of the oblique waves' amplitude is equal to or less than that of the plane wave. Results show that the fundamental basically follows the linear theory, while the subharmonic follows an exponential-of-an-exponential growth.

Mankbadi, Reda R.↗

Characterization of a hot-film probe for hypersonic flow

The critical issues concerning the application of constant-temperature hot-film anemometry to hypersonic flow are reviewed and extended. Mass-flux static calibrations were conducted in a Mach 10 helium flow, while mass-flux and total-temperature static calibrations were made in a Mach 6 air flow. In addition, comparative hot-film/hot-wire turbulence measurements were made in a Mach 11 helium boundary layer to provide insight into the dynamic response of the hot film. The measurements indicate that substrate conduction 'losses' dominate the static response of the hot-film probe, thus resulting in poor sensitivity to mass-flux and total temperature. Furthermore, it has been found that it is not possible to isolate mass-flux fluctuations at high overheat ratios for the current hot-film design. Thus, the sapphire-substrate hot-film anemometer is a robust, high-bandwidth instrument limited to qualitative transition and turbulence measurements. Finally, the extension of this technique to providing quantitative information is dependent upon the development of lower thermal-conductivity substrate materials.

Sheplak, M.↗

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗

Asymmetric Base-Bleed Effect on Aerospike Plume-Induced Base-Heating Environment During Power-Pack Out

A computational heat transfer design methodology was developed to study tbe dual-engine linear aerospike plume-induced base-heating environment during one power-pack out, in ascent flight. It includes a three-dimensional, finite volume, viscous, chemically reacting, and pressure-based computational fluid dynamics formulation, a special base-bleed boundary condition, and a three-dimensional, finite volume, and spectral-line-based weighted-sum-of-gray-gases absorption computational radiation heat transfer formulation. A separate radiation model was used for diagnostic purposes. The computational methodology was systematically benchmarked. in this study, near-base radiative heat fluxes were computed, and they compared well with those measured during static linear aerospike engine tests. The base-heating environment of 18 trajectory points secected from three power-pack out scenarios was computed. The computed asymmetric base-heating physics were analyzed. The power-pack out condition has the most impact on convective base heating when it happens early in flight. The soume of its impact comes from the asymmetric and reduced base bleed.

Wang, Ten-See↗