Search NASASearch

SEARCH · Search NASA

Results for “Shifted Boundary Method”

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 37 records · Page 2

Calculated blade response at high tip-speed ratios

This paper presents the initial results of a digital computer study of helicopter rotor-blade-motion stability for a broad range of forward speeds encompassing proposed conventional and compound helicopter designs. The analysis treats the general case of nonlinear, coupled flapping and lagging motions of hinged rigid blades. The results of the study to date indicate that blade-motion stability boundaries cannot be specified for a given blade design in terms of a fixed value of rotor tip-speed ratio. The stability boundaries can shift significantly within the desired operational speed range, with the stability limits depending upon the rotor loading, the mag­nitude of the disturbance encountered by the blade, and the blade position at the instant the external disturbance is encountered. The sample methods suggested for improving blade-motion stability include increasing the effec­tive hinge spring restraint by incorporating pitch-flap coupling or by using hingeless or teetering rotor systems. Reducing rotor loading also has a beneficial effect on the blade-motion stability. The results indi­cate that methods used to deal with the blade-motion stability problem can be expected to diminish the overall vibration levels of the helicopter rotor system across the entire design speed range.

TIP SPEED

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model

Parabolic harmonics in magnetospheric modeling - The main dipole and the ring current

A new representation is developed for that portion of the magnetospheric magnetic field which is produced by currents on the magnetopause or outside it. For a magnetopause which is a paraboloid of rotation and which confines all field lines of magnetospheric sources, the fields due to magnetopause currents are derived for two sources, a dipole and a two-parameter model of the ring current, each arbitrarily inclined to the earth-sun line. Methods are also developed for accomodating changes in the size of the magnetospheric cavity and for using nonparaboloidal boundaries, open at the nightside. Neutral points and polar cusps are derived, and in agreement with observations it is found that a growth of the ring current shifts them equatorward. By extrapolation from these results, it is argued that the observation of such shifts during magnetic storms may help resolve the relative importance of two processes proposed for magnetic storms.

Stern, D. P.

Doubling fusion power with volumetric optimization in magnetic confinement fusion devices

A technique, volumetric power optimization, is presented for enhancing the power output of magnetic confinement fusion devices. Applied to a tokamak, this approach involves shifting the burning plasma region to a larger plasma volume while introducing minimal perturbations to the plasma boundary shape. This edge perturbation—squareness—is analogous to pinching and stretching the edge boundary. Stability calculations confirm that this edge alteration is compatible with maintaining plasma stability. This optimization method for optimizing fusion power output could improve the performance of magnetic confinement fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

The Near Noise Field of Static Jets and Some Model Studies of Devices for Noise Reduction

An experimental study is presented of the pressure fluctuations near jet exhaust streams made during unchoked operation of a turbojet engine and a 1-inch-diameter high-temperature model jet and during choked operation of various sizes of model jets with unheated air. The tests for unchoked operation indicate a random spectrum of rather narrow band width which varies in frequency content with axial position along the jet. Pressure surveys from the model tests along lines parallel to the 15 degree jet boundary indicate that the station of greatest pressure fluctuations is determined by the jet velocity and the radial distance, with a tendency of the maximum to shift downstream as either parameter is increased. From model tests the magnitude of the fluctuations appears to increase as about the second power of jet velocity at points just outside the jet boundary and as increasingly higher powers of jet velocity as distance from the boundary is increased. A laboratory method of noise reduction with model jets was found to produce large decreases in the magnitude of the lower-frequency components of the spectra and thereby also to reduce the total radiated energy. Choked operation of model jets with unheated air indicates the appearance of a discrete-frequency component of very large magnitude. Shadowgraph records of the flow show that this condition is associated with the appearance of flow formations suggestive of partly formed toroidal vortices in the vicinity of the shocks. Elimination of these formations is found to eliminate the discrete component and thereby to reduce the overall noise level.

Lassiter, Leslie W

An atomic cluster expansion potential for twisted multilayer graphene

Twisted multilayer graphene, characterized by its moiré patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. Machine learning interatomic potentials (MLIPs) are a promising approach to model such systems. Our work develops a method to generate training and test datasets for fitting MLIPs that capture all possible misalignments but remain small-scale to facilitate efficient data generation and parameter estimation. To achieve this, we generate configurations with periodic boundary conditions suitable for density functional theory calculations, and then introduce an internal twist and shift within those supercell structures. Using this technique, supplemented with an active learning workflow, we fit an Atomic Cluster Expansion potential for simulating twisted multilayer graphene and test it for accuracy and robustness on a range of simulation tasks.

2D materials

Acceleration of convergence and spectrum transformation of implicit finite difference operators associated with Navier-Stokes equations

Implicit difference operator spectra are presently computed by applying eigensystem analysis techniques to finite-difference formulations of 2D Euler and Navier-Stokes equations, and attention is given to these iterative methods' convergence and stability characteristics by taking into account the effects of grid geometry, time-step, numerical viscosity, and boundary conditions. On the basis of the eigenvalue distributions for various flow configurations, the feasibility of applying such convergence-acceleration techniques as eigenvalue annihilation and relaxation is discussed. Spectrum-shifting is applied to NASA-Ames' ARC2D flow code, achieving a 20-33 percent efficiency.

Saleem, M.

Scale-Resolving Simulations of Low-Pressure Turbine Cascades with Wall Roughness Using A Spectral-Element Method

The accurate prediction of wall-roughness effects in turbomachinery is becoming critical as turbine designers address airfoil surface quality and degradation concerns arising from the shift to advanced ceramic matrix composite (CMC) or additively-manufactured airfoils operating in higher temperature environments. In this paper, a recently developed computational capability for accurate and efficient scale-resolving simulations of turbomachinery is extended to analyze the boundary- layer separation and transition characteristics in a rough-wall low-pressure turbine (LPT) cascade. The computational capability is based on an entropy-stable discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in an efficient manner for a modern high performance computer architecture. Results from the scale-resolving simulations of both smooth and rough airfoil cascades are presented and compared to previous experiments and numerical simulations. The results show that the suction surface boundary layer undergoes laminar separation, transition, and turbulent reattachment for the smooth airfoil cascade, while in the presence of roughness the separation and transition behavior of the suction surface boundary layer is substantially modified. The differences between the smooth and rough airfoil cascades are then highlighted by a detailed analysis of their respective turbulent flow fields.

Spectral-Element

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics

Creation of a Weather Drivers Test Suite for Inclusion in ASHRAE Standard 140

Weather conditions are an important boundary condition for building performance simulation (BPS) calculations. For existing test cases in ASHRAE Standard 140 "Method of Test for Evaluating Building Performance Simulation Software" (ANSI/ASHRAE 2020), it was assumed that the software being tested could adequately read and interpret the weather data in the provided standard weather files. As differences between the programs have been reduced and as more programs have shifted to sub-hourly time steps this assumption has become more stretched. To address these concerns a new test suite testing a program's ability to read and interpret the data from a standard weather file was developed. The purpose of the test suite is to test the use of the typical data used from standard weather files.

54 ENVIRONMENTAL SCIENCES

A Numerical and Experimental Investigation of Flow Separation Control over a Wall-Mounted Hump Model

Numerical simulations and wind tunnel experiments were performed to investigate flow separation and its control over a wall-mounted hump model. Three-dimensional, unsteady fluid dynamic simulations were used to capture the relevant flow physics of unsteady flow separation and its control via active and passive flow control methods. Surface pressure measurements and surface oilflow visualization were obtained in the experiments. Current measurements documented a considerably thinner incoming boundary layer, but similar pressure distribution, compared to the reference data reported in the CFD (Computational Fluid Dynamics) Validation (CFDVAL 2004) Experiments. Consistent with the data available in the literature, numerical simulations predicted a longer separation bubble, which manifested itself as a shift in the pressure distributions. Streamwise vortices generated by passive vortex generators increased the suction pressure and provided substantial pressure recovery. Numerical simulations of vortex generators qualitatively predicted the effects of control, and revealed patterns of attached/separated flow regions. On the other hand, active flow control using steady discrete jets reduced, but was not able to eliminate, flow separation for the tested configuration. Although the major disagreement between the numerical and experimental results was in the pressure recovery region due to the inaccurate prediction of separation bubbles, the relative pressure recoveries from the baselines were comparable. Surface oilflow visualization and simulated surface streamlines were qualitatively in good agreement revealing the key flow features.

Koklu, Mehti

Pair Identity and Smooth Variation Rules Applicable for the Spectroscopic Parameters of H2O Transitions Involving High-J States

Two basic rules (i.e. the pair identity and the smooth variation) applicable for H2O transitions involving high-J states have been discovered. The origins of these rules are the properties of the energy levels and wavefunctions of H2O states with the quantum number J above certain boundaries. As a result, for lines involving high-J states in individually defined groups, all their spectroscopic parameters (i.e. the transition wavenumber, intensity, pressure-broadened half-width, pressure-induced shift, and temperature exponent) must follow these rules. One can use these rules to screen spectroscopic data provided by databases and to identify possible errors. In addition, by using extrapolation methods within the individual groups, one is able to predict the spectroscopic parameters for lines in this group involving very high-J states. The latter are required in developing high-temperature molecular spectroscopic databases such as HITEMP.

Ma, Q.

Improved Discrimination of Volcanic Complexes, Tectonic Features, and Regolith Properties in Mare Serenitatis from Earth-Based Radar Mapping

Radar images at 70 cm wavelength show 4-5 dB variations in backscatter strength within regions of relatively uniform spectral reflectance properties in central and northern Mare Serenitatis, delineating features suggesting lava flow margins, channels, and superposition relationships. These backscatter differences are much less pronounced at 12.6 cm wavelength, consistent with a large component of the 70 cm echo arising from the rough or blocky transition zone between the mare regolith and the intact bedrock. Such deep probing is possible because the ilmenite content, which modulates microwave losses, of central Mare Serenitatis is generally low (2-3% by weight). Modeling of the radar returns from a buried interface shows that an average regolith thickness of 10m could lead to the observed shifts in 70 cm echo power with a change in TiO2 content from 2% to 3%. This thickness is consistent with estimates of regolith depth (10-15m) based on the smallest diameter for which fresh craters have obvious blocky ejecta. The 70 cm backscatter differences provide a view of mare flow-unit boundaries, channels, and lobes unseen by other remote sensing methods. A localized pyroclastic deposit associated with Rima Calippus is identified based on its low radar echo strength. Radar mapping also improves delineation of units for crater age dating and highlights a 250 km long, east-west trending feature in northern Mare Serenitatis that we suggest is a large graben flooded by late-stage mare flows.

Radar echoes

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5 angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for proper validation of the CFD.

Maldonado, Daniel

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5o angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for more thorough validation of the CFD.

Maldonado, Daniel