Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moving grid methods”

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 73 records · Page 4

Application of a novel finite difference method to dynamic crack problems

A versatile finite difference method (HEMP and HEMP 3D computer programs) was developed originally for solving dynamic problems in continuum mechanics. It was extended to analyze the stress field around cracks in a solid with finite geometry subjected to dynamic loads and to simulate numerically the dynamic fracture phenomena with success. This method is an explicit finite difference method applied to the Lagrangian formulation of the equations of continuum mechanics in two and three space dimensions and time. The calculational grid moves with the material and in this way it gives a more detailed description of the physics of the problem than the Eulerian formulation.

Chen, Y. M.↗

Applications of the Space-Time Conservation Element and Solution Element (CE/SE) Method to Computational Aeroacoustic Benchmark Problems

The Internal Propagation problems, Fan Noise problem, and Turbomachinery Noise problems are solved using the space-time conservation element and solution element (CE/SE) method. The problems in internal propagation problems address the propagation of sound waves through a nozzle. Both the nonlinear and linear quasi 1D Euler equations are solved. Numerical solutions are presented and compared with the analytical solution. The fan noise problem concerns the effect of the sweep angle on the acoustic field generated by the interaction of a convected gust with a cascade of 3D flat plates. A parallel version of the 3D CE/SE Euler solver is developed and employed to obtain numerical solutions for a family of swept flat plates. Numerical solutions for sweep angles of 0, 5, 10, and 15 deg are presented. The turbomachinery problems describe the interaction of a 2D vortical gust with a cascade of flat-plate airfoils with/without a downstream moving grid. The 2D nonlinear Euler Equations are solved and the converged numerical solutions are presented and compared with the corresponding analytical solution. All the comparisons demonstrate that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple and efficient manner. Furthermore, the simple non-reflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well in 1D, 2D and 3D problems.

Wang, Xiao-Yen↗

ENSAERO - A multidisciplinary program for fluid/structural interaction studies of aerospace vehicles

A computational procedure is developed that uses a moving zonal grid concept to model complex flexible aerospace vehicles. The Euler/Navier-Stokes equations are used to model the flow, and computations are made using efficient methods based on both central and upwind schemes. The structure is represented by a finite element method which can model general aerospace vehicles. Provisions are made to accommodate other disciplines such as controls and thermal loads. The code is capable of computing unsteady flows on flexible wings with vortical flows. Adaptation of this procedure for parallel processing and validation for complete aerospace configurations is in progress.

Guruswamy, G. P.↗

Enhancing Scalability for FUN3D Rotorcraft Simulations with Yoga: an Overset Grid Assembler

FUN3D, an unstructured grid Navier-Stokes CFD code is capable of overset grid simu- lations, but does not have an internal method for assembling an overset grid system from a group of component grids. FUN3D currently relies on the third party codes Suggar++ and DiRTlib to perform domain assembly and provide intergrid connectivity for overset simulations. Rotorcraft simulations with moving, deforming blades require domain assembly and mesh deformation at each time step. For these simulations, the three primary drivers of computational cost for each time step are: deforming the mesh, performing domain assembly, and performing subiterations of the flow solver. FUN3D exhibits strong and weak scalability for the flow solver subiterations and the mesh deformation. However, FUN3D is currently hardwired directly to the serial version of Suggar++, which has a fixed cost for a given mesh system. Therefore, domain assembly begins to dominate the total cost of each time step as grid systems become larger. An integrated method for parallel domain assembly is presented that addresses scalability for large grid systems. Restructuring within FUN3D to accommodate integrated domain assembly is also discussed, which could enable use of the parallel Suggar++ library.

Cameron T Druyor↗

Unsteady airfoil flow solutions on moving zonal grids

Euler and Navier-Stokes solutions for airfoil flows on zonal grids are presented. The governing equations are solved with an implicit, iterative, factorized numerical scheme. The inviscid fluxes are examined with a third-order accurate upwind method. Zonal grid solutions are compared with experimental measurements for flows over airfoils at fixed angles of incidence. The computed unsteady solutions for rapidly pitching and oscillating airfoils are in good agreement with experiments.

Cricelli, Antonio S.↗

A finite-volume numerical method to calculate fluid forces and rotordynamic coefficients in seals

A numerical method to calculate rotordynamic coefficients of seals is presented. The flow in a seal is solved by using a finite-volume formulation of the full Navier-Stokes equations with appropriate turbulence models. The seal rotor is perturbed along a diameter such that the position of the rotor is a sinusoidal function of time. The resulting flow domain changes with time, and the time-dependent flow in the seal is solved using a space conserving moving grid formulation. The time-varying fluid pressure reaction forces are then linked with the rotor center displacement, velocity and acceleration to yield the rotordynamic coefficients. Results for an annular seal are presented, and compared with experimental data and other more simplified numerical methods.

Athavale, M. M.↗

The need for enhanced initial moisture information in simulations of a complex summertime precipitation event

Initial simulations of the June 28, 1986 Cooperative Huntsville Meteorological Experiment case illustrate the need for mesoscale moisture information in a summertime situation in which deep convection is organized by weak large scale forcing. A methodology is presented for enhancing the initial moisture field from a combination of IR satellite imagery, surface-based cloud observations, and manually digitized radar data. The Mesoscale Atmospheric Simulation Model is utilized to simulate the events of June 28-29. This procedure insures that areas known to have precipitation at the time of initialization will be nearly saturated on the grid scale, which should decrease the time needed by the model to produce the observed Bonnie (a relatively weak hurricane that moved on shore two days before) convection. This method will also result in an initial distribution of model cloudiness (transmissivity) that is very similar to that of the IR satellite image.

Waight, Kenneth T., III↗

Identification of the dynamics of a two-dimensional grid structure using least square lattice filters

It is doubtful whether the dynamics of large space structures (LSS) can be predicted well enough for control system design applications. Hence, dynamic modeling from on-orbit measurements followed by a modification of the control system is of interest, taking into account the utilization of adaptive control concepts. The present paper is concerned with the model determination phase of the adaptive control problem. Using spectral decoupling to determine mode shapes, mode frequency and damping data can be obtained with the aid of an equation error parameter identification method. This method employs a second-order auto-regressive moving average (ARMA) model to represent the natural mode amplitudes. The discussed procedure involves an extension of the application of the least square lattice filter in system identification to a nonintegral, two-dimensional grid structure made of overlapping bars.

Montgomery, R. C.↗

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

A verification of unsteady Navier-Stokes solutions around oscillating airfoils

A finite difference solution code for the two dimensional Navier-Stokes equations was combined with a moving-grid system. The thin layer Navier-Stokes equations with a turbulence model are solved in a time-accurate manner in order to study the unsteady aerodynamics around airfoils undergoing small amplitude pitching or heaving motions in the transonic regime. The accuracy of the solutions obtained by the use of the present moving-grid technqiue is investigated. The effects of the minimum grid size and the integrating time-step size on the solutions are also checked. Some of the solutions obtained by the present method are compared with experimental results. It is demonstrated that the unsteady aerodynamics around oscillating airfoils can be predicted fairly well by the present code for cases in which the dynamic angle of attack or displacement is small.

Nakamichi, J.↗

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

Stirling Analysis Comparison of Commercial Versus High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's compact scheme and Dyson's Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model with sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗

Stirling Analysis Comparison of Commercial vs. High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/ proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's Compact scheme and Dyson s Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model although sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗

Application of a streamwise upwind algorithm for unsteady transonic computations over oscillating wings

A new streamwise upwind algorithm has been derived to compute unsteady flows with a moving grid system and applied to compute flows over oscillating wings at transonic Mach numbers. Comparisons have been made between results obtained from this upwind algorithm, using both temporally nonconservative- and conservative-implicit methods, with the results obtained from a central-difference method, and also with experimental data. The results show (1) the efficiency and practicality of the temporally nonconservative implicit solver and (2) the robustness and accuracy of the upwind method for unsteady computations compared to the central-difference method.

Obayashi, Shigeru↗

The Monotonic Lagrangian Grid for Fast Air-Traffic Evaluation

This paper describes the continued development of a dynamic air-traffic model, ATMLG, intended for rapid evaluation of rules and methods to control and optimize transport systems. The underlying data structure is based on the Monotonic Lagrangian Grid (MLG), which is used for sorting and ordering positions and other data needed to describe N moving bodies, and their interactions. In ATMLG, the MLG is combined with algorithms for collision avoidance and updating aircraft trajectories. Aircraft that are close to each other in physical space are always near neighbors in the MLG data arrays, resulting in a fast nearest-neighbor interaction algorithm that scales as N. In this paper, we use ATMLG to examine how the ability to maintain a required separation between aircraft decreases as the number of aircraft in the volume increases. This requires keeping track of the primary and subsequent collision avoidance maneuvers necessary to maintain a five mile separation distance between all aircraft. Simulation results show that the number of collision avoidance moves increases exponentially with the number of aircraft in the volume.

Alexandrov, Natalia↗

Freestream capturing for moving coordinates in three dimensions

A detailed formula is developed that can be used in both finite-volume (FV) and finite-difference (FD) methods for constructing freestream capturing metrics in space and time. It is shown that, considering an FV cell on the FD grid, the freestream capturing metrics in space and time can be constructed from the FD formulation. The approach is costly but guarantees the global conservation for an arbitrary motion of the grid.

Obayashi, Shigeru↗

Development of a time-dependent incompressible Navier-Stokes solver based on a fractional-step method

The main goals are the development, validation, and application of a fractional step solution method of the time-dependent incompressible Navier-Stokes equations in generalized coordinate systems. A solution method that combines a finite volume discretization with a novel choice of the dependent variables and a fractional step splitting to obtain accurate solutions in arbitrary geometries is extended to include more general situations, including cases with moving grids. The numerical techniques are enhanced to gain efficiency and generality.

Rosenfeld, Moshe↗

Helicopter Rotor Blade Computation in Unsteady Flows Using Moving Overset Grids

An overset grid thin-layer Navier-Stokes code has been extended to include dynamic motion of helicopter rotor blades through relative grid motion. The unsteady flowfield and airloads on an AH-IG rotor in forward flight were computed to verify the methodology and to demonstrate the method's potential usefulness towards comprehensive helicopter codes. In addition, the method uses the blade's first harmonics measured in the flight test to prescribe the blade motion. The solution was impulsively started and became periodic in less than three rotor revolutions. Detailed unsteady numerical flow visualization techniques were applied to the entire unsteady data set of five rotor revolutions and exhibited flowfield features such as blade vortex interaction and wake roll-up. The unsteady blade loads and surface pressures compare well against those from flight measurements. Details of the method, a discussion of the resulting predicted flowfield, and requirements for future work are presented. Overall, given the proper blade dynamics, this method can compute the unsteady flowfield of a general helicopter rotor in forward flight.

Ahmad, Jasim↗