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 451 records · Page 25

Multi-Fidelity Simulation of a Turbofan Engine With Results Zoomed Into Mini-Maps for a Zero-D Cycle Simulation

A Zero-D cycle simulation of the GE90-94B high bypass turbofan engine has been achieved utilizing mini-maps generated from a high-fidelity simulation. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled 3D computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the 3D component models are integrated into the cycle model via partial performance maps generated from the CFD flow solutions using one-dimensional mean line turbomachinery programs. This paper highlights the generation of the high-pressure compressor, booster, and fan partial performance maps, as well as turbine maps for the high pressure and low pressure turbine. These are actually "mini-maps" in the sense that they are developed only for a narrow operating range of the component. Results are compared between actual cycle data at a take-off condition and the comparable condition utilizing these mini-maps. The mini-maps are also presented with comparison to actual component data where possible.

Turner, Mark G.↗

Computational Fluid Dynamics (CFD) Design of a Blended Wing Body (BWB) with Boundary Layer Ingestion (BLI) Nacelles

A study is being conducted to improve the propulsion/airframe integration for the Blended Wing-Body (BWB) configuration with boundary layer ingestion nacelles. TWO unstructured grid flow solvers, USM3D and FUN3D, have been coupled with different design methods and are being used to redesign the aft wing region and the nacelles to reduce drag and flow separation. An initial study comparing analyses from these two flow solvers against data from a wind tunnel test as well as predictions from the OVERFLOW structured grid code for a BWB without nacelles has been completed. Results indicate that the unstructured grid codes are sufficiently accurate for use in design. Results from the BWB design study will be presented.

Morehouse, Melissa B.↗

New Challenges of Cloud Remote Sensing from Space

Clouds are recognized as a major source of uncertainties in predicting climate change. Lack of observational constraints on cloud processes has hindered the development of more reliable climate and weather models. With the advanced NASA EOS and A-Train sensors, we now have better knowledge about vertical distribution, water content, and occurrence frequency of global cloudiness. New ice water content (IWC) measurements from CloudSat and MLS have led to several improvements in model physics and parameterization schemes. MISR and GPS high-resolution data start to reveal deep insights on cloud processes and dynamicS in the planetary boundary layer (PBL). However, Earth sciences are still facing urgent needs to measure cloud microphysical properties, interactions between clouds and aerosol/precipitation, and processes that are coupled in 3-D space but not adequately sampled by the A-Train curtains or by the 2-D imageries. Submillimeter- wave, multi-angle imaging, and GPS radio occultation have emerged as promising remote sensing techniques to meet the challenges and enable new sciences for mid-tropospheric and PBL clouds.

Wu, Dong L.↗

Multiscale Analysis of Structurally-Graded Microstructures Using Molecular Dynamics, Discrete Dislocation Dynamics and Continuum Crystal Plasticity

A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.

Saether, Erik↗

Switchback Boundary Dissipation and Relative Age

We examine Parker Solar Probe (PSP) magnetic field and plasma observations during its first encounter with the Sun in early 2018 November. During this perihelion time, impulsive reversals in the magnetic field, called "switchbacks," were found in the data set characterized by a quick rotation in B along with a simultaneous increase in solar wind flow. In this work, we examine the structure and morphology of 920 switchback boundaries as PSP enters and exits the structures, specifically looking for evidence of boundary degradation, dissipation, and associated ultralow frequency (ULF) magnetic wave activity. We find that boundaries with the most abrupt, step-function-like change in Br and Vr also show little evidence of dissipation and ULF wave activity. In contrast, there is a set of boundaries that appears highly degraded with ULF magnetic activity in the vicinity of the boundary. We thus infer that the steep, step-like boundaries with little ULF activity are relatively young in comparison to the degraded boundaries. The distribution in relative ages suggests that the switchback boundary formation process is dynamic and evolving, even occurring near the PSP observation point inside of 40 Rs.

W. M. Farrell↗

Review of numerical procedures for computational surface thermochemistry

Models and equations for surface thermochemistry and near-surface thermophysics of aerodynamically heated thermal protection materials are reviewed, with particular emphasis on computational boundary conditions for surface mass and energy transfer. The surface energy and mass balances, coupled with an appropriate ablation or surface catalysis model, provide complete thermochemical boundary conditions for a true multidisciplinary solution of the fully coupled fluid-dynamics/solid mechanics problem. Practical approximate solutions can be obtained by using a detailed model with full thermophysics for either the solid or fluid phase and a semianalytic method for the other half of the problem. A significant increase in the state-of-the-art in aerothermal computational fluid dynamics is possible by uniting computational fluid dynamic (CFD) methodology with surface thermochemistry boundary conditions and the heat-balance-integral method.

SOLID MECHANICS↗

Unsteady laminar boundary-layer calculations on oscillating configurations including backflow. Part 2: Airfoil in high-amplitude pitching motion. Dynamic stall

A previously developed finite-difference procedure for calculating unsteady, incompressible, laminar boundary layers on an oscillating flat plate is applied to a wing section undergoing high-amplitude pitching oscillations about various mean incidences. To start the entire boundary-layer calculation, appropriate initial conditions and outer boundary conditions are specified, using a stagnation-point fixed frame of reference. The breakdown of the numerical calculation procedure in the x,t-domain is interpreted to coincide with unsteady separation. Details of the boundary-layer behavior in the vicinity of separation are investigated, and a close analogy between the present results and those for a three-dimensional steady separation is found.

Geissler, W.↗

Full Vorticity Budget of the Arabian Sea from a 0.1° Ocean Model: Sverdrup Dynamics, Rossby Waves, and Nonlinear Eddy Effects

Abstract The Arabian Sea, influenced by the Indian monsoon, has many unique features, including its basin-scale seasonally reversing surface circulation and the Great Whirl, a seasonal anticyclonic system appearing during the southwest monsoon close to the western boundary. To establish a comprehensive dynamical picture of the Arabian Sea, we utilize numerical model output and design a full vorticity budget that includes a fully decomposed nonlinear term. The ocean general circulation model has 0.1° resolution and is mesoscale eddy-resolving in the region. In the western boundary current system, we highlight the role of nonlinear eddies in the life cycle of the Great Whirl. The nonlinear eddy term is of leading-order importance in this feature’s vorticity balance. Specifically, it contributes to the Great Whirl’s persistence in boreal fall after the weakening of the southwesterly winds. In the open ocean, Sverdrup dynamics and annual Rossby waves are found to dominate the vorticity balance; the latter is considered as a key factor in the formation of the Great Whirl and the seasonal reversal of the western boundary current. In addition, we discuss different forms of vertically integrated vorticity equations in the model and argue that the bottom pressure torque term can be interpreted analogously as friction in the western boundary and vortex stretching in the open ocean.

58 GEOSCIENCES↗

Computation of incompressible viscous flows through artificial heart devices with moving boundaries

The extension of computational fluid dynamics techniques to artificial heart flow simulations is illustrated. Unsteady incompressible Navier-Stokes equations written in 3-D generalized curvilinear coordinates are solved iteratively at each physical time step until the incompressibility condition is satisfied. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. The efficiency and robustness of the time accurate formulation of the algorithm are tested by computing the flow through model geometries. A channel flow with a moving indentation is computed and validated with experimental measurements and other numerical solutions. In order to handle the geometric complexity and the moving boundary problems, a zonal method and an overlapping grid embedding scheme are used, respectively. Steady state solutions for the flow through a tilting disk heart valve was compared against experimental measurements. Good agreement was obtained. The flow computation during the valve opening and closing is carried out to illustrate the moving boundary capability.

Kiris, Cetin↗

Effects of boundary hydraulics, dissolved oxygen, and dissolved organic carbon on growth and death dynamics of aerobic microbes in riverbed dune-induced hyporheic zones

Surface and groundwater interact in the hyporheic zone beneath and adjacent to rivers in the presence of a diverse microbial community. Heterotrophic bacteria mediate a range of environmentally important reactions, yet few studies have quantified bacterial growth and death dynamics in the hyporheic zone, and none have systematically analyzed their response to variations in hydraulic or chemical conditions. Here we used MODFLOW and SEAM3D to simulate hydraulics; dissolved oxygen (DO) and dissolved organic carbon (DOC) transport; and aerobic microbial metabolism, growth, and death in hyporheic zones induced by riverbed dunes. We ran simulations both with and without growth/death processes, and varied hydraulic parameters and DO/DOC boundary concentrations. Microbial biomass reached steady state (t = 3 days) in every simulation, at which time there was greater biomass and DOC biodegradation rates in the hyporheic flowcell (300% and 85% higher for the base case, respectively) when accounting for microbial growth dynamics. This occurred as microbial biomass tailored its spatial distribution to the availability of DO and DOC, demonstrating the importance of simulating growth/death processes. Biomass generally increased with hyporheic flow cell area as upwelling groundwater decreased. When varying surface water DO and DOC source concentrations relative to the base case, the greatest effect on biomass occurred when increasing DOC and decreasing DO. We determined minimum DO and DOC steady-state concentrations required for microbial growth, but the minimums were not absolute or related by stoichiometry. Increasing DOC created a smaller area of microbes with higher concentrations relative to the base case. Increasing DO slightly increased the area occupied by microbes while keeping the total biomass nearly constant. Overall, microbial growth and death dynamics depend on DO and DOC availability in the hyporheic zone, which is dependent on DOC/DO boundary concentrations and hyporheic flow paths, and in turn the hydraulic interaction between surface water and groundwater.

54 ENVIRONMENTAL SCIENCES↗

Methodology for simulation of unsteady flows to determine the time-dependent interference between stationary and moving boundaries

A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, Kamakhya Prasad↗

On the differential representation and color-kinematics duality of AdS boundary correlators

The AdS boundary correlators and their dual correlation functions of boundary operators have been the main dynamic observables of the holographic duality relating a bulk AdS theory and a boundary conformal field theory. We show that tree-level AdS boundary correlators for generic states can be expressed as nonlocal differential operators of a certain structure acting on contact Witten diagrams. We further write the boundary correlators in a form that is very similar to flat space amplitudes, with Mandelstam variables replaced by certain combinations of single-state conformal generators, prove that all tree-level AdS boundary correlators have a differential representation, and detail the conversion of such differential expressions to position space. We illustrate the construction through the computation of the boundary correlators of scalars coupled to gluons and gravitons; when converted to position space, they reproduce known results. Color-kinematics duality and BCJ relations can be defined in analogy with their flat space counterparts, and are respected by the scalar correlators with a gluon exchange. We also discuss potential approaches to the double copy and find that its direct generalization may require nontrivial extensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling the Effect of Hydrogen on Tungsten Grain Boundary Properties

The objective of this study is to assess the effect of hydrogen isotopes on mechanical properties of tungsten that are relevant to first-wall applications in a tokamak. We focus on the effects of hydrogen segregation to grain boundaries and its impact on the grain boundary structure and properties. Through molecular dynamics and other atomistic techniques, we are able to make predictions of grain boundary structure and properties related to it, including grain boundary mobility and impurity diffusivity along the grain boundary. The grain boundary structure can change significantly and abruptly as conditions change, including due to heating and radiation damage. These grain boundary phase transformations can have a pronounced effect on mobility and impurity diffusivity, even when no hydrogen is present. How this behavior changes in the presence of hydrogen is not known. The phase may affect the level of hydrogen segregation and hydrogen storage, and the hydrogen may induce significant grain boundary structure changes. Our objective is to predict whether they do. The results inform thermomechanical modeling of tungsten at tokamak conditions.

36 MATERIALS SCIENCE↗

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young’s modulus, Poisson’s ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young’s modulus and Poisson’s ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a “safe zone” regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

bimodular materials↗

Fully implicit shock tracking

A tracking method for solving hyperbolic free boundary problems arising in fluid dynamics is presented. The principal feature of the new approach is that the implicit form is used to simultaneously treat both interior points and boundary conditions. In particular, the location of the free boundary (shock) surface is treated implicitly and coupled with all other unknowns. The method is presented in the context of unsteady one-dimensional flow in a variable area duct with an internal shock wave. The fully implicit method and other strategies for advancing the shock are compared for computing a steady solution via a time asymptotic approach. Issues regarding extension of the method to multiple dimensions are also discussed.

Bell, J. B.↗

Atmospheric Boundary-Layer and Flutter Computations Using CFD Model of the Transonic Dynamics Tunnel

This paper presents two Computational Fluid Dynamics (CFD) models of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas as the test medium. The first computational model describes the generation of an atmospheric-boundary-layer (ABL) profile inside the tunnel. An ABL, which includes both a wind profile and turbulence content, is one of the aerodynamic characteristics affecting the occurrence of wind-induced oscillations for a launch vehicle sitting on a pad. The challenging part of this analysis was modeling the turbulent flow inside the tunnel. This is due to the special ABL-generating hardware that was installed at the entrance of the TDT test section in order to change the downstream velocity profile and to introduce velocity fluctuations into the flow. The second CFD model builds on the computational aeroelastic results that were generated in support of the second Aeroelastic Prediction Workshop (AePW) for the NASA Benchmark Supercritical Wing (BSCW) configuration. During the AePW, the wing-only configuration (classical free-air model) was analyzed. In the current study, the flutter computations were conducted on the configuration as it was mounted in the TDT during the experiment. This includes the wing attached to the splitter plate that was attached to the wind-tunnel walls. The preliminary results show that the wind-tunnel walls marginally affect flutter prediction.

Pawel Chwalowski↗

Experimental evaluation of the aerodynamic damping of skin panels at low supersonic Mach numbers.

The aerodynamic damping of the first three vibration modes of a rectangular panel exposed to an attached turbulent boundary layer has been measured over the Mach number range of 1.10 to 1.40. Data showing that aerodynamic damping is a strong function of Mach number, dynamic pressure, mode shape, and boundary layer thickness are presented. At low supersonic Mach numbers, aerodynamic damping was found to be 20 to 100 times as large as common values of structural damping. Means of calculating aerodynamic damping, including effects of the boundary layer profile, are discussed. Calculated aerodynamic damping is compared with experimental data.

Muhlstein, L., Jr.↗