Combining asymptotic analysis and weighted L2 stability estimates, the convergence of lattice Boltzmann methods for the approximation of 1D convection-diffusion-reaction equations is proved. Unlike previous approaches, the proof does not require transformations to equivalent macroscopic equations.

Currently, there exists a lack of confidence in the computational simulation of multiple body high-speed air delivered systems. Of particular interest is the ability to accurately predict the dispersion pattern of these systems under various deployment configurations. Classical engineering-level methods may not be able to predict these patterns with adequate confidence due primarily to accuracy errors attributable to reduced order modeling. In the current work, a new collision modeling capability has been developed to enable multiple-body proximate-flight simulation in the Loci/CHEM framework. This approach maintains high-fidelity aerodynamics and incorporates six degrees of freedom modeling with collision response, and is well-suited for simulation of a large number of projectiles. The proposed simulation system is intended to capture the strong interaction phase early in the projectile deployment, with subsequent transfer of projectile positions and flight states to the more economical engineering-level methods. Collisions between rigid bodies are modeled using an impulse-based approach with either an iterative propagation method or a simultaneous method. The latter is shown to be more accurate and robust for cases involving multiple simultaneous collisions as it eliminates the need to sort and resolve the collisions sequentially. The implementation of both the collision detection methodology and impact mechanics are described in detail with validation studies to demonstrate the efficiency and accuracy of the developed technologies. The studies chronologically detail the findings for simulating simple impacts and collisions between multiple bodies with aerodynamic interference effects.

An accurate cartesian method is devised to simulate incompressible viscous flows past an arbitrary moving body. The Navier-Stokes equations are spatially discretized onto a fixed Cartesian mesh. The body is taken into account via the ghost-cell method and the so-called penalty method, resulting in second-order accuracy in velocity. The accuracy and the efficiency of the solver are tested through two-dimensional reference simulations. To show the versatility of this scheme we simulate a three-dimensional self propelled jellyfish prototype.

Recently, a new differential discontinuous formulation for conservation laws named the Correction Procedure via Reconstruction (CPR) is developed, which is in-spired by several other discontinuous methods such as the discontinuous Galerkin (DG), the spectral volume (SV)/spectral difference (SD) methods. All of them can be unified under the CPR formulation, which is relatively simple to implement due to its finite-difference-like framework. In this paper, a different discontinuous solution space including both polynomial and Fourier basis functions on each element is employed to compute broad-band waves. Free-parameters introduced in the Fourier bases are optimized to minimize both dispersion and dissipation errors through a wave propagation analysis. The optimization procedure is verified with a mesh resolution analysis. Numerical results are presented to demonstrate the performance of the optimized CPR formulation.

In this paper we demonstrate the accuracy and robustness of combining the advection upwind splitting method (AUSM), specifically AUSM+-UP, with high-order upwind-biased interpolation procedures, the weighted essentially non-oscillatory (WENO-JS) scheme and its variations, and the monotonicity preserving (MP) scheme, for solving the Euler equations. MP is found to be more effective than the three WENO variations studied. AUSM+-UP is also shown to be free of the so-called “carbuncle” phenomenon with the high-order interpolation. The characteristic variables are preferred for interpolation after comparing the results using primitive and conservative variables, even though they require additional matrix-vector operations. Results using the Roe flux with an entropy fix and the Lax-Friedrichs approximate Riemann solvers are also included for comparison. In addition, four reflective boundary condition implementations are compared for their effects on residual convergence and solution accuracy. Finally, a measure for quantifying the efficiency of obtaining high order solutions is proposed; the measure reveals that a maximum return is reached after which no improvement in accuracy is possible for a given grid size.