Reduction of computer usage costs in predicting unsteady aerodynamic loadings caused by control surface motions: Analysis and results

Results of theoretical and numerical investigations conducted to develop economical computing procedures were applied to an existing computer program that predicts unsteady aerodynamic loadings caused by leading and trailing edge control surface motions in subsonic compressible flow. Large reductions in computing costs were achieved by removing the spanwise singularity of the downwash integrand and evaluating its effect separately in closed form. Additional reductions were obtained by modifying the incremental pressure term that account for downwash singularities at control surface edges. Accuracy of theoretical predictions of unsteady loading at high reduced frequencies was increased by applying new pressure expressions that exactly satisified the high frequency boundary conditions of an oscillating control surface. Comparative computer result indicated that the revised procedures provide more accurate predictions of unsteady loadings as well as providing reduction of 50 to 80 percent in computer usage costs

Comparison of macro- and microscopic solutions of the Riemann problem I. Supercritical shock tube and expansion into vacuum

The Riemann problem is a fundamental concept in the development of numerical methods for the macroscopic flow equations. It allows the resolution of discontinuities in the solution, such as shock waves, and provides a powerful tool for the construction of numerical flux functions. A natural extension of the Riemann problem involves two phases, a liquid and a vapour phase which undergo phase change at the material boundary. For this problem, we aim at a comparison with the macroscopic solution from molecular dynamics simulations. In this work, as a first step, the macroscopic solution of two important Riemann problem scenarios, the supercritical shock tube and the expansion into vacuum, were compared to microscopic solutions produced by molecular dynamics simulations. High fidelity equations of state were used to accurately approximate the material behaviour of the model fluid. The results of both scenarios compare almost perfect with each other. During the vacuum expansion, the fluid obtained a state of non-equilibrium, where the microscopic and macroscopic solutions start to diverge. A limiting case was shown, where liquid droplets appeared in the expansion fan, which was approximated by the macroscopic solution, assuming an undercooled vapour.DFG, 84292822, TRR 75: Tropfendynamische Prozesse unter extremen Umgebungsbedingunge

(Un)anticipated Technological Change in an Endogenous Growth Model

This paper examines numerically the impact of a negative exogenous shock to marginal productivity (such as ecological government regulation that becomes effective at some point in time) in an endogenous finite-time growth model with sluggish reallocation of human capital. The policy can be anticipated or unanticipated by firms, and it can also be announced but not implemented. It turns out that these frictions have a very strong long-run effect on output, consumption and on the optimal allocation of capital and labor in particular. The qualitative properties relate to homogenous labor models with positive productivity shocks. The problem is thus to maximize a function of a continuous system, where the system is subject to frictions and stepwise changes; for such a problem the application of calculus of variations necessary conditions is problematic. A numerical optimization method, which has had much success on qualitatively similar problems in engineering, has been employed.two-sector endogenous growth model; unanticipated and anticipated technological change; frictions in reallocation of human capital; Runge-Kutta parallel shooting algorithm

Air-breathing hypersonic vehicle guidance and control studies; An integrated trajectory/control analysis methodology: Phase 1

A tool which generates optimal trajectory/control histories in an integrated manner is generically adapted to the treatment of single-stage-to-orbit air-breathing hypersonic vehicles. The methodology is implemented as a two point boundary value problem solution technique. Its use permits an assessment of an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A simpler analysis strategy that partitions the trajectory into several boundary condition matched segments is also included to construct preliminary trajectory and control history representations with less computational burden than is required for the overall flight profile assessment. A demonstration was accomplished using a tabulated example (winged-cone accelerator) vehicle model that is combined with a newly developed multidimensional cubic spline data smoothing routine. A constrained near-fuel-optimal trajectory, imposing a dynamic pressure limit of 1000 psf, was developed from horizontal takeoff to 20,000 ft/sec relative air speed while aiming for a polar orbit. Previously unspecified propulsive discontinuities were located. Flight regimes demanding rapid attitude changes were identified, dictating control effector and closed-loop controller authority was ascertained after evaluating effector use for vehicle trim. Also, inadequacies in vehicle model representations and specific subsystem models with insufficient fidelity were determined based on unusual control characteristics and/or excessive sensitivity to uncertainty

A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: Triangular grids

A novel wetting and drying treatment for second-order Runge-Kutta discontinuous Galerkin (RKDG2) methods solving the non-linear shallow water equations is proposed. It is developed for general conforming two-dimensional triangular meshes and utilizes a slope limiting strategy to accurately model inundation. The method features a non-destructive limiter, which concurrently meets the requirements for linear stability and wetting and drying. It further combines existing approaches for positivity preservation and well-balancing with an innovative velocity-based limiting of the momentum. This limiting controls spurious velocities in the vicinity of the wet/dry interface. It leads to a computationally stable and robust scheme -- even on unstructured grids -- and allows for large time steps in combination with explicit time integrators. The scheme comprises only one free parameter, to which it is not sensitive in terms of stability. A number of numerical test cases, ranging from analytical tests to near-realistic laboratory benchmarks, demonstrate the performance of the method for inundation applications. In particular, super-linear convergence, mass-conservation, well-balancedness, and stability are verified

Modelling Shear Flows with SPH and Grid Based Methods

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.Comment: 16 pages, 17 figure