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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Asymmetric Base-Bleed Effect on Aerospike Plume-Induced Base-Heating Environment

A computational heat transfer design methodology was developed to study the 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 selected 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 source of its impact comes from the asymmetric and reduced base bleed.

Wang, Ten-See↗
NASA NTRS2011

Space Shuttle Ascent Flight Design Process: Evolution and Lessons Learned

The Space Shuttle Ascent Flight Design team is responsible for defining a launch to orbit trajectory profile that satisfies all programmatic mission objectives and defines the ground and onboard reconfiguration requirements for this high-speed and demanding flight phase. This design, verification and reconfiguration process ensures that all applicable mission scenarios are enveloped within integrated vehicle and spacecraft certification constraints and criteria, and includes the design of the nominal ascent profile and trajectory profiles for both uphill and ground-to-ground aborts. The team also develops a wide array of associated training, avionics flight software verification, onboard crew and operations facility products. These key ground and onboard products provide the ultimate users and operators the necessary insight and situational awareness for trajectory dynamics, performance and event sequences, abort mode boundaries and moding, flight performance and impact predictions for launch vehicle stages for use in range safety, and flight software performance. These products also provide the necessary insight to or reconfiguration of communications and tracking systems, launch collision avoidance requirements, and day of launch crew targeting and onboard guidance, navigation and flight control updates that incorporate the final vehicle configuration and environment conditions for the mission. Over the course of the Space Shuttle Program, ascent trajectory design and mission planning has evolved in order to improve program flexibility and reduce cost, while maintaining outstanding data quality. Along the way, the team has implemented innovative solutions and technologies in order to overcome significant challenges. A number of these solutions may have applicability to future human spaceflight programs.

Picka, Bret A.↗
NASA NTRS2015

Non-Linear Slosh Damping Model Development and Validation

Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can lead to significant savings by reducing the number and size of slosh baffles in liquid propellant tanks.

Yang, H. Q.↗
DOE OSTI2024

Multi-scale modeling of fuel fragmentation and microstructural evolution

An overview of NEAMS program activities to model fuel fragmentation and microstructual evolution to a group of industry attendees at a meeting hosted by EPRI. Topics include atomistic methods to determine grain boundary strength, phase-field fracture modeling to determine critical pressure at grain boundaries, and coupled phase-field-cluster dynamics to model microstructural evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1993

Models of thermal/chemical boundary layer convection: Potential application to Venus

The upper boundary layer of Venus is comprised of at least two distinct chemical components, mantle and crust. Fluid dynamical models of convection within Venus' mantle were primarily of the thermal boundary layer type. Models assessing the ability of convective mantle flows to deform the crust were undertaken, but models exploring the effects of a variable thickness crust on mantle convection were largely lacking. A Venusian crust of variable thickness could couple back into, and alter, the mantle flow patterns that helped create it, leading to deformation mechanisms not predicted by purely thermal boundary layer convection models. This possibility is explored through a finite element model of thermal/chemical boundary layer convection. Model results suggest that a crust of variable thickness can serve as a mantle flow driver by perturbing lateral temperature gradients in the upper mantle. Resulting mantle flow is driven by the combination of free convective and nonuniform crustal distribution. This combination can lead to a flow instability manifest in the occurrence of episodic mantle lithosphere subduction initiated at the periphery of a crustal plateau. The ability of a light, near surface, chemical layer to potentially alter mantle flow patterns suggest that mantle convection and the creation and/or deformation of such a chemical layer may be highly nonseparable problems on time scales of 10(exp 8) years.

Lenardic, A.↗
NASA NTRS2017

A Tree-Ring Chronology and Paleoclimate Record for the Younger Dryas-Early Holocene Transition from Northeastern North America

Spruce and tamarack logs dating from the Younger Dryas and Early Holocene (YDEH; approx. 12.9 - 11.3k cal a BP) were found at Bell Creek in the Lake Ontario lowlands of the Great Lakes region, North America. A 211-year tree-ring chronology dates to approx. 11 755 -11 545 cal a BP, across the YDEH transition. A 23-year period of higher year-to-year ring-width variability dates to around 11 650 cal a BP, infers strong regional climatic perturbations and may represent the end of the YD. Tamarack and spruce were dominant species throughout the YD - EH interval at the site, indicating that boreal conditions persisted into the EH, in contrast to geographical regions immediately south and east of the lowlands, but consistent with the Great Lakes interior lowlands. This infers that Bell Creek was at the eastern boundary of a boreal ecotone, perhaps a result of its lower elevation and the non-analog dynamics of the Laurentide Ice Sheet. This finding suggests that the ecotone boundary extended farther east during the YD - EH transition than previously thought.

dendroclimatology;northeastern North America;spruc↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

The Dynamical Role of the Chesapeake Bay on the Local Ozone Pollution Using Mesoscale Modeling—A Case Study

This study investigated the dynamic influence of the Chesapeake Bay (CB) on local ozone (O3) concentration and distribution using a weather forecasting model. The Weather Research and Forecasting model coupled with Chemistry (WRF–Chem) was employed to simulate O3 production and transportation near the CB. Baseline (water) as well as sensitivity (nowater) model experiments of bay circulation were carried out with and without bay water by changing the water surface from water to land (loam). First, the model performance simulating O3 was evaluated using the baseline experiment results and AirNow surface wind and O3 observations. The results showed that the model overestimates surface O3 by up to 20–30%. Further, the comparisons of the baseline and sensitivity experiments revealed higher O3 mixing ratios, primarily due to the resulting bay breeze circulation. These increases, after considering model overestimation, represent a mean bay dynamics circulation-induced contribution of up to 10% at night and 5% during the day. Furthermore, the boundary layer over northern CB, where it is at its narrowest width, was higher (by 1.2 km on average) during daytime due to higher surface temperatures observed. The boundary layer depth difference between the northern, central, and southern regions of the bay leads to a differential in the role of bay circulation dynamics in the observed O3 increase. The relatively wider swath of water surface over southern CB resulted in a lower boundary layer depth and stronger breeze circulation and this circulation contributed to O3 concentrations. Moreover, since the case selected has a minimal bay breeze circulation, the associated surface ozone enhancements represent what is expected at least at a minimum.

Chesapeak Bay↗
DOE OSTI2022

Characterization of Disconnection Content Associated with Grain Boundary Motion in Tungsten

The objective of this study is to develop a new computational methodology to predict the motion of grain boundaries in tungsten associated with recrystallization. This approach accounts recently discovered high-temperature grain boundary phases in tungsten and other metals and alloys. This paper reports the initial work in the project that developed a novel technique for molecular dynamics to characterize disconnections, dislocation-like defects in grain boundaries that mediate grain boundary motion. Grain boundary properties are needed for a thermomechanical model of recrystallization of tungsten for magnetic fusion applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A simulation study of test fixture attachment points for vibration test optimization (MEEN 5440 Final Project)

The aerospace industry uses vibration shaker tables to perform component durability testing. These dynamic environments tests often have a hard time replicating the service environment of a component due to a mismatch in the applied force or the boundary condition. Figure 1 shows the desired outcome of a dynamic environments test. The measured responses in a lab test, ideally, should match the measured response of the service environment field test. One way to address the mismatch between field environment and lab environment is with test fixture design. The test fixture attaches the component under test to the shaker table. Single-axis shakers are the industry standard type of shaker table, and, therefore, test fixtures are often rigid in order to minimize cross-axis responses.

42 ENGINEERING↗
NASA NTRS2011

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS1999

A Generalized Wall Function

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Shih, Tsan-Hsing↗
NASA NTRS2017

Development of a Flow Field for Testing a Boundary-Layer-Ingesting Propulsor

The test section of the 8- by 6-Foot Supersonic Wind Tunnel at NASA Glenn Research Center was modified to produce the test conditions for a boundary-layer-ingesting propulsor. A test was conducted to measure the flow properties in the modified test section before the propulsor was installed. Measured boundary layer and freestream conditions were compared to results from computational fluid dynamics simulations of the external surface for the reference vehicle. Testing showed that the desired freestream conditions and boundary layer thickness could be achieved; however, some non-uniformity of the freestream conditions, particularly the total temperature, were observed.

embedded propulsion↗
NASA NTRSDate not supplied

PIV Measurements of Shock-Wave/Boundary-Layer Interactions at Mach 2.5 in a Circular Test Section

The particle image velocimetry (PIV) measurement technique has been applied to make detailed measurements in an impinging shock-wave/boundary-layer interaction (SWBLI) flowfield in a test section of circular cross-section. By virtue of this test section geometry, the flowfields studied are free from effects of the complex coupling of sidewall and corner interaction into the primary impinging shock interaction typical of a rectangular test section. The cone-cylinder shock generator hardware was arranged in both coaxial alignment with the test section to generate axisymmetric SWBLIs and in offset positions to generate swept SWBLIs. The angle of the conical section was also varied to alter the strength of the impinging shock wave. The application of the PIV technique in this test section required meticulous design to provide adequate optical access without disturbing the SWBLI flow. Details of the stereoscopic PIV apparatus and the velocity measurements obtained for a Mach 2.5 freestream are presented. Key features of the mean velocity and turbulent stress fields are highlighted. This data set should enable validation of computational fluid dynamics results for these flowfields.

Shockwave Boundary Layer Interaction↗
NASA NTRSDate not supplied

PIV Measurements of Shock-Wave/Boundary-Layer Interactions at Mach 2.5 in a Circular Test Section

The particle image velocimetry (PIV) measurement technique has been applied to make detailed measurements in an impinging shock-wave/boundary-layer interaction (SWBLI) flowfield in a test section of circular cross-section. By virtue of this test section geometry, the flowfields studied are free from effects of the complex coupling of sidewall and corner interaction into the primary impinging shock interaction typical of a rectangular test section. The cone-cylinder shock generator hardware was arranged in both coaxial alignment with the test section to generate axisymmetric SWBLIs and in offset positions to generate swept SWBLIs. The angle of the conical section was also varied to alter the strength of the impinging shock wave. The application of the PIV technique in this test section required meticulous design to provide adequate optical access without disturbing the SWBLI flow. Details of the stereoscopic PIV apparatus and the velocity measurements obtained for a Mach 2.5 freestream are presented. Key features of the mean velocity and turbulent stress fields are highlighted. This data set should enable validation of computational fluid dynamics results for these flowfields.

Shockwave Boundary Layer Interaction↗
NASA NTRS2006

Three-Dimensional Thermal Boundary Layer Corrections for Circular Heat Flux Gauges Mounted in a Flat Plate with a Surface Temperature Discontinuity

Three-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been performed in an effort to determine thermal boundary layer correction factors for circular convective heat flux gauges (such as Schmidt-Boelter and plug type)mounted flush in a flat plate subjected to a stepwise surface temperature discontinuity. Turbulent flow solutions with temperature-dependent properties are obtained for a free stream Reynolds number of 1E6, and freestream Mach numbers of 2 and 4. The effect of gauge diameter and the plate surface temperature have been investigated. The 3-D CFD results for the heat flux correction factors are compared to quasi-21) results deduced from constant property integral solutions and also 2-D CFD analysis with both constant and variable properties. The role of three-dimensionality and of property variations on the heat flux correction factors has been demonstrated.

Kandula, M.↗
DOE OSTI2021

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗