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 541 records · Page 30

A parametric study of supersonic laminar flow for swept wings using linear stability analysis

A parametric study to predict the extent of laminar flow on the upper surface of a generic swept-back wing (NACA 64A010 airfoil section) at supersonic speeds was conducted. The results were obtained by using surface pressure predictions from an Euler/Navier-Stokes computational fluid dynamics code coupled with a boundary layer code, which predicts detailed boundary layer profiles, and finally with a linear stability code to determine the extent of laminar flow. The parameters addressed are Reynolds number, angle of attack, and leading-edge wing sweep. The results of this study show that an increase in angle of attack, for specific Reynolds numbers, can actually delay transition. Therefore, higher lift capability, caused by the increased angle of attack, as well as a reduction in viscous drag due to the delay in transition is possible for certain flight conditions.

Cummings, Russell M.↗

Hot Gas Testing Results of Stagnant Volume Filling Through a Back-Filled RSRM Nozzle Joint

Silicone rubber is back-filled into RSW nozzle joints after assembly. A possible artifact of this process is the formation of small pathways through this filter leading to sensitive sealing materials within the joint. At motor ignition, hot gases fill the stagnant volumes within the joint through this path. A series of 14 test has been completed studying this scenario and providing anchoring data for thermal/flow models. Parameters such as gas path cross section, gas path length, gas path materials, fill volume size, and post path gas spreading prior to unpingement on seal material, have been investigated. Tests were accomplished using geometry similar to RSRM nozzle joint 4 with attached volumes replicating the free volume and flow friction in the actual hardware. The test hardware simulated 8 inches of the full-scale circumference. Testing has pointed to changes required in model boundary condition assumptions and gas dynamics corrections for gas paths of this size and geometry. Areas ,where this date has provided improvement in analysis models will be covered as well as model inadequacies that require separate specialized efforts. Questions remaining after this testing and a possible direction for future testing will be suggested.

Prince, Andrew S.↗

Atmospheric Science Measurements by the EOS Geoscience Laser Altimeter System

Scheduled for Launch in July 2001, the Geoscience Laser Altimeter System (GLAS) is to be the first satellite instrument to provide full global lidar profiling of clouds and aerosol in the earth's atmosphere. GLAS is an EOS program instrument that is on its own satellite, now called the Ice, Cloud and land Elevation Satellite. The instrument is both a surface laser ranging system and an atmospheric profiling lidar. A most important surface measurement for the instrument is to study the change in the mass balance of the polar ice sheets by measuring the change in regional altitudes to an accuracy of 1.5 cm per year. The strategy to combine the surface measurement with a Cloud and aerosol lidar profiling mission is based on the compatibility of the altimetry instrument requirements with those for the required lidar measurements. The primary atmospheric science goal of the GLAS cloud and aerosol measurement is to determine the radiative forcing and vertically resolved atmospheric heating rate due to cloud and aerosol by directly observing the vertical structure and magnitude of cloud and aerosol parameters that are important for the radiative balance of the earth-atmosphere system, but which are ambiguous or impossible to obtain from existing or planned passive remote sensors. A further goal is to directly measure the height of atmospheric transition layers (inversions) which are important for dynamics and mixing, the planetary boundary layer and lifting condensation level.

Spinhirne, James↗

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗

Stochastic Optimal Control via Bellman's Principle

This paper presents a method for finding optimal controls of nonlinear systems subject to random excitations. The method is capable to generate global control solutions when state and control constraints are present. The solution is global in the sense that controls for all initial conditions in a region of the state space are obtained. The approach is based on Bellman's Principle of optimality, the Gaussian closure and the Short-time Gaussian approximation. Examples include a system with a state-dependent diffusion term, a system in which the infinite hierarchy of moment equations cannot be analytically closed, and an impact system with a elastic boundary. The uncontrolled and controlled dynamics are studied by creating a Markov chain with a control dependent transition probability matrix via the Generalized Cell Mapping method. In this fashion, both the transient and stationary controlled responses are evaluated. The results show excellent control performances.

Crespo, Luis G.↗

TURBO-AE: An Aeroelastic Code for Propulsion Applications

NASA's Advanced Subsonic Technology (AST) program is developing new technologies to increase the fuel efficiency of commercial aircraft engines, improve the safety of engine operation, and reduce engine emissions and noise. With the development of new designs for ducted fans, compressors, and turbines to achieve these goals, a basic aeroelastic requirement is that there should be no flutter or high resonant blade stresses in the operating regime. To verify the aeroelastic soundness of these designs, we need an accurate prediction and analysis code. Such a two-dimensional viscous propulsion aeroelastic code, named TURBO-AE, is being developed at the NASA Lewis Research Center. The TURBO-AE aeroelastic code is based on a three-dimensional unsteady aerodynamic Euler/Navier-Stokes turbomachinery code TURBO, developed under a grant from NASA Lewis. TURBO-AE can model viscous flow effects that play an important role in certain aeroelastic problems, such as flutter with flow separation (or stall flutter) and flutter in the presence of shock and boundary-layer interaction. The structural dynamics representation of the blade in the TURBO-AE code is based on a normal mode representation. A finite element analysis code, such as NASTRAN, is used to calculate in-vacuum vibration modes and the associated natural frequency. A work-per-cycle approach is used to determine aeroelastic (flutter) stability. With this approach, the motion of the blade is prescribed to be a harmonic vibration in a specified in vacuum normal mode. The aerodynamic forces acting on the vibrating blade and the work done by these forces on the vibrating blade during a cycle of vibration are calculated. If positive work is being done on the blade by the aerodynamic forces, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the amplitude of the blade's oscillation. Initial calculations have been done for a configuration representative of the Energy Efficient Engine fan rotor. The accompanying figure shows the work-per-cycle after each cycle of vibration. It can be seen that the work-per-cycle does not vary much after the fourth cycle. The negative sign of the converged work-per-cycle shows that the fan blade is dynamically stable and will not flutter. TURBO-AE will provide a useful aeroelastic prediction/analysis capability for engine manufacturers. It will reduce design cycle times by allowing new blade designs to be verified for aeroelastic soundness before blades are built and tested. With this prediction capability, it will be possible to build thinner, lighter, and faster rotating blades without encountering aeroelastic problems like stall flutter and high-cycle fatigue due to forced vibrations.

Bakhle, Milind A.↗

Transient Side Load Analysis of Out-of-Round Film-Cooled Nozzle Extensions

There was interest in understanding the impact of out-of-round nozzle extension on the nozzle side load during transient startup operations. The out-of-round nozzle extension could be the result of asymmetric internal stresses, deformation induced by previous tests, and asymmetric loads induced by hardware attached to the nozzle. The objective of this study was therefore to computationally investigate the effect of out-of-round nozzle extension on the nozzle side loads during an engine startup transient. The rocket engine studied encompasses a regeneratively cooled chamber and nozzle, along with a film cooled nozzle extension. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and transient inlet boundary flow properties derived from an engine system simulation. Six three-dimensional cases were performed with the out-of-roundness achieved by three different degrees of ovalization, elongated on lateral y and z axes: one slightly out-of-round, one more out-of-round, and one significantly out-of-round. The results show that the separation line jump was the primary source of the peak side loads. Comparing to the peak side load of the perfectly round nozzle, the peak side loads increased for the slightly and more ovalized nozzle extensions, and either increased or decreased for the two significantly ovalized nozzle extensions. A theory based on the counteraction of the flow destabilizing effect of an exacerbated asymmetrical flow caused by a lower degree of ovalization, and the flow stabilizing effect of a more symmetrical flow, created also by ovalization, is presented to explain the observations obtained in this effort.

Wang, Ten-See↗

How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?

Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.

Xu, Kuan-Man↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

RELAP5-3D in the OECD-NEA HTGR Thermal Hydraulics Benchmark

The OECD-NEA HTGR Thermal Hydraulics benchmark uses data from the High Temperature Test Facility for a series of code-to-code and code-to-data exercises aimed at improving the state of knowledge on existing thermal hydraulics modeling and simulation tools for prismatic HTGR applications. This benchmark includes problems representing hot gas mixing in the lower plenum, the depressurized conduction cooldown accident and the pressurized conduction cooldown accident. This presentation discusses the benchmark and RELAP5-3D's role in that benchmark, including as a tool for predicting behavior in the core and as a tool for providing boundary conditions to computational fluid dynamics analysis in the lower plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Increased Mach Number Capability for the NASA Glenn 10x10 Supersonic Wind Tunnel

Computational simulations and wind tunnel testing were conducted to explore the operation of the Abe Silverstein Supersonic Wind Tunnel at the NASA Glenn Research Center at test section Mach numbers above the current limit of Mach 3.5. An increased Mach number would enhance the capability for testing of supersonic and hypersonic propulsion systems. The focus of the explorations was on understanding the flow within the second throat of the tunnel, which is downstream of the test section and is where the supersonic flow decelerates to subsonic flow. Methods of computational fluid dynamics (CFD) were applied to provide details of the shock boundary layer structure and to estimate losses in total pressure. The CFD simulations indicated that the tunnel could be operated up to Mach 4.0 if the minimum width of the second throat was made smaller than that used for previous operation of the tunnel. Wind tunnel testing was able to confirm such operation of the tunnel at Mach 3.6 and 3.7 before a hydraulic failure caused a stop to the testing. CFD simulations performed after the wind tunnel testing showed good agreement with test data consisting of static pressures along the ceiling of the second throat. The CFD analyses showed increased shockwave boundary layer interactions, which was also observed as increased unsteadiness of dynamic pressures collected in the wind tunnel testing.

Diffuser↗

Probing chiral Kitaev spin liquids via dangling boundary fermions

Identifying experimental probes capable of diagnosing extreme quantum behavior is widely regarded as one of the foremost challenges in modern condensed matter physics. Here, we propose a novel approach for detecting chiral Kitaev spin-liquid states through measurements of the local dynamical spin structure factor on the boundary using scanning tunneling microscopy (STM). We specifically focus on unpaired (“dangling”) Majorana fermions, which naturally emerge along boundaries of Kitaev spin liquids, and can serve as indicators of chiral boundary modes under broad conditions, thereby offering a clear signature of these exotic quantum states.

Zhang, Shang-Shun [University of Tennessee, Knoxvi↗

Free-Spinning-Tunnel Investigation of a 1/20-Scale Model of an Unswept-Wing Jet-Propelled Trainer Airplane

A flutter analysis employing the kernel function for three- dimensional, subsonic, compressible flow is applied to a flutter-tested tail surface which has an aspect ratio of 3.5, a taper ratio of 0.15, and a leading-edge sweep of 30 deg. Theoretical and experimental results are compared at Mach numbers from 0.75 to 0.98. Good agreement between theoretical and experimental flutter dynamic pressures and frequencies is achieved at Mach numbers to 0.92. At Mach numbers from 0.92 to 0.98, however, a second solution to the flutter determinant results in a spurious theoretical flutter boundary which is at a much lower dynamic pressure and at a much higher frequency than the experimental boundary.

Bowman, James S., Jr.↗

Subsonic Kernel-Function Flutter Analysis of a Highly Tapered Tail Surface and Comparison with Experimental Results

A flutter analysis employing the kernel function for three-dimensional, subsonic, compressible flow is applied to a flutter-tested tail surface which has an aspect ratio of 3.5, a taper ratio of 0.15, and a leading-edge sweep of 30 deg. Theoretical and experimental results are compared at Mach numbers from 0.75 to 0.98. Good agreement between theoretical and experimental flutter dynamic pressures and frequencies is achieved at Mach numbers to 0.92. At Mach numbers from 0.92 to 0.98, however, a second solution to the flutter determinant results in a spurious theoretical flutter boundary which is at a much lower dynamic pressure and at a much higher frequency than the experimental boundary.

Walberg, Gerald D.↗

Implementation of the block-Krylov boundary flexibility method of component synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based on the boundary flexibility vectors of the component. This algorithm is not load-dependent, is applicable to both fixed and free-interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. This methodology was implemented in the MSC/NASTRAN normal modes solution sequence using DMAP. The accuracy is found to be comparable to that of component synthesis based upon normal modes. The block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem.

Carney, Kelly S.↗

Hydrodynamics of two-dimensional compressible fluid with broken parity: Variational principle and free surface dynamics in the absence of dissipation

In this paper, we consider an isotropic compressible nondissipative fluid with broken parity subject to free surface boundary conditions in two spatial dimensions. The hydrodynamic equations describing the bulk dynamics of the fluid and the free surface boundary conditions depend explicitly on the parity-breaking nondissipative odd viscosity term. We construct an effective action which gives both bulk hydrodynamic equations and free surface boundary conditions. The free surface boundary conditions require an additional boundary term in the action which resembles a 1+1D chiral boson field coupled to the background geometry. We solve the linearized hydrodynamic equations for the deep water case and derive the dispersion of chiral surface waves. We show that in the long-wavelength limit the flow profile exhibits an oscillating vortical boundary layer near the free surface. The layer thickness is controlled by the ratio between the odd viscosity (ν o ) and the sound velocity (c s ), δ ~ ν o /c s . In the incompressible limit, c s → ∞, the vortical boundary layer becomes singular with the vorticity within the layer diverging as ω ~ c s . The boundary layer is formed by odd viscosity coupling the divergence of velocity ∇ • v to vorticity ∇ × v. It results in nontrivial chiral free surface dynamics even in the absence of external forces. The structure of the odd-viscosity-induced boundary layer is very different from the conventional free surface boundary layer associated with dissipative shear viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗