Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic boundary”

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 289 records · Page 16

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.↗

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↗

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↗

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.↗

Computational fluid dynamics overview

The hierarchy of codes; boundary layers and turbulence models; uniqueness and artificial viscosity; efficiency of solution algorithms; data display and analysis; and cost of experimental validation are briefly examined.

Bobbitt, Percy J.↗

Two-Dimensional Thermal Boundary Layer Corrections for Convective Heat Flux Gauges

This work presents a CFD (Computational Fluid Dynamics) study of two-dimensional thermal boundary layer correction factors for convective heat flux gauges mounted in flat plate subjected to a surface temperature discontinuity with variable properties taken into account. A two-equation k - omega turbulence model is considered. Results are obtained for a wide range of Mach numbers (1 to 5), gauge radius ratio, and wall temperature discontinuity. Comparisons are made for correction factors with constant properties and variable properties. It is shown that the variable-property effects on the heat flux correction factors become significant

Kandula, Max↗

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

Orbiter entry trajectory corridors: 32000 pound payload, 67.5 percent center of gravity

Thermal and equilibrium glide boundaries were used to analyze and/or design shuttle orbiter entry trajectories. Plots are presented of orbiter thermal and equilibrium glide boundaries in the drag/mass-relative velocity dynamic pressure-relative velocity, and altitude-relative velocity planes for an orbiter having a 32,000 pound payload and a 67.5% center of gravity location. These boundaries were defined for control points 1 through 4 of the shuttle orbiter for 40 deg-30 deg and 38 deg-28 deg ramped angle of attack entry profiles and 40 deg, 38 deg, 35 deg, 30 deg, 28 deg, and 25 deg constant angle of attack entry profiles each at 20 deg, 15 deg, and 10 deg constant body flap settings.

Treybig, J. H.↗

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

Embedded function methods for compressible high speed turbulent flow

Fundamental issues relating to compressible turbulent flow are addressed. The focus has been on developing methods and testing concepts for attached flows rather than trying to force a conventional law of the wall into a zone of backflow. Although the dynamics of the near-wall flow in an attached turbulent boundary layer are relatively well documented, the dynamical features of a zone of reversed turbulent flow are not, nor are they well understood. Incompressibility introduces effects and issues that have been dealt with only marginally in the literature, therefore, the present work has been focussed on attached high-speed flows. The wall function method has been extended up through the supersonic to hypersonic speeds. Algorithms have been successfully introduced into the code that calculates the flow all the way to the wall, and testing is being carried out for progressively more complex flow situations.

Walker, J. D. A.↗

TopMaker: Technique Developed for Automatic Multiblock Topology Generation Using the Medial Axis

The TopMaker technique was developed in an effort to reduce the time required for grid generation in complex numerical studies. Topology generation accounts for much of the man-hours required for structured multiblock grids. With regard to structured multiblock grids, topology refers to how the blocks are arranged and connected. A two-dimensional multiblock topology generation technique has been developed at the NASA Glenn Research Center. Very general configurations can be addressed by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, a multiblock topology can be generated without user intervention. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns, where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in computational fluid dynamics, where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid-generation process.

Rigby, David L.↗