Engineering PapersSearch

Engineering topics

Mazaheri, Alireza

Publications and source records attributed to Mazaheri, Alireza.

Hyperbolic Method for Dispersive PDEs: Same High-Order of Accuracy for Solution, Gradient, and Hessian

In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.

Mazaheri, Alireza

High-Order Residual-Distribution Schemes for Discontinuous Problems on Irregular Triangular Grids

In this paper, we develop second- and third-order non-oscillatory shock-capturing hyperbolic residual distribution schemes for irregular triangular grids, extending our second- and third-order schemes to discontinuous problems. We present extended first-order N- and Rusanov-scheme formulations for hyperbolic advection-diffusion system, and demonstrate that the hyperbolic diffusion term does not affect the solution of inviscid problems for vanishingly small viscous coefficient. We then propose second- and third-order blended hyperbolic residual-distribution schemes with the extended first-order Rusanov-scheme. We show that these proposed schemes are extremely accurate in predicting non-oscillatory solutions for discontinuous problems. We also propose a characteristics-based nonlinear wave sensor for accurately detecting shocks, compression, and expansion regions. Using this proposed sensor, we demonstrate that the developed hyperbolic blended schemes do not produce entropy-violating solutions (unphysical stocks). We then verify the design order of accuracy of these blended schemes on irregular triangular grids.

Mazaheri, Alireza

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza

First-Order Hyperbolic System Method for Time-Dependent Advection-Diffusion Problems

A time-dependent extension of the first-order hyperbolic system method for advection-diffusion problems is introduced. Diffusive/viscous terms are written and discretized as a hyperbolic system, which recovers the original equation in the steady state. The resulting scheme offers advantages over traditional schemes: a dramatic simplification in the discretization, high-order accuracy in the solution gradients, and orders-of-magnitude convergence acceleration. The hyperbolic advection-diffusion system is discretized by the second-order upwind residual-distribution scheme in a unified manner, and the system of implicit-residual-equations is solved by Newton's method over every physical time step. The numerical results are presented for linear and nonlinear advection-diffusion problems, demonstrating solutions and gradients produced to the same order of accuracy, with rapid convergence over each physical time step, typically less than five Newton iterations.

Mazaheri, Alireza

Methodology for Flight Relevant Arc-Jet Testing of Flexible Thermal Protection Systems

A methodology to correlate flight aeroheating environments to the arc-jet environment is presented. For a desired hot-wall flight heating rate, the methodology provides the arcjet bulk enthalpy for the corresponding cold-wall heating rate. A series of analyses were conducted to examine the effects of the test sample model holder geometry to the overall performance of the test sample. The analyses were compared with arc-jet test samples and challenges and issues are presented. The transient flight environment was calculated for the Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Earth Atmospheric Reentry Test (HEART) vehicle, which is a planned demonstration vehicle using a large inflatable, flexible thermal protection system to reenter the Earth's atmosphere from the International Space Station. A series of correlations were developed to define the relevant arc-jet test environment to properly approximate the HEART flight environment. The computed arcjet environments were compared with the measured arc-jet values to define the uncertainty of the correlated environment. The results show that for a given flight surface heat flux and a fully-catalytic TPS, the flight relevant arc-jet heat flux increases with the arc-jet bulk enthalpy while for a non-catalytic TPS the arc-jet heat flux decreases with the bulk enthalpy.

Mazaheri, Alireza

Influence of Coupled Radiation and Ablation on the Aerothermodynamic Environment of Planetary Entry Vehicles

A review of recently published coupled radiation and ablation capabilities involving the simulation of hypersonic flowfields relevant to Earth, Mars, or Venus entry is presented. The three fundamental mechanisms of radiation coupling are identified as radiative cooling, precursor photochemistry, and ablation-radiation interaction. The impact of these mechanisms are shown to be significant for a 3 m radius sphere entering Earth at hypothetical Mars return conditions (approximately 15 km/s). To estimate the influence precursor absorption on the radiative flux for a wide range of conditions, a simplified approach is developed that requires only the non-precursor solution. Details of a developed coupled ablation approach, which is capable of treating both massively ablating flowfields in the sublimation regime and weakly ablating diffusion Climited oxidation cases, are presented. A review of the two primary uncoupled ablation approximations, identified as the blowing correction and film coefficient approximations, is made and their impact for Earth and Mars entries is shown to be significant for recession and convective heating predictions. Fully coupled ablation and radiation simulations are presented for the Mars return sphere throughout its entire trajectory. Applying to the Mars return sphere the Pioneer- Venus heritage carbon phenolic heatshield, which has properties available in the open literature, the differences between steady state ablation and coupling to a material response code are shown to be significant.

Johnston, Christopher O.

LAURA Users Manual: 5.5-64987

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintain ability by eliminating the requirement for problem dependent recompilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multi-physics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the Fun3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza

High-Energy Atmospheric Reentry Test Aerothermodynamic Analysis

This paper presents an assessment of the aerothermodynamic environment around an 8.3 meter High Energy Atmospheric Reentry Test (HEART) vehicle. This study generated twelve nose shape configurations and compared their responses at the peak heating trajectory point against the baseline nose shape. The heat flux sensitivity to the angle of attack variations are also discussed. The possibility of a two-piece Thermal Protection System (TPS) design at the nose is also considered, as are the surface catalytic affects of the aeroheating environment of such configuration. Based on these analyses, an optimum nose shape is proposed to minimize the surface heating. A recommendation is also made for a two-piece TPS design, for which the surface catalytic uncertainty associated with the jump in heating at the nose-IAD juncture is reduced by a minimum of 93%. In this paper, the aeroshell is assumed to be rigid and the inflatable fluid interaction effect is left for future investigations

Mazaheri, Alireza

HEART Aerothermodynamic Analysis

This paper presents an assessment of the aerothermodynamic environment around an 8.3 meter High Energy Atmospheric Reentry Test (HEART) vehicle. This study generated twelve nose shape configurations and compared their responses at the peak heating trajectory point against the baseline nose shape. The heat flux sensitivity to the angle of attack variations are also discussed. The possibility of a two-piece Thermal Protection System (TPS) design at the nose is also considered, as are the surface catalytic affects of the aeroheating environment of such configuration. Based on these analyses, an optimum nose shape is proposed to minimize the surface heating. A recommendation is also made for a two-piece TPS design, for which the surface catalytic uncertainty associated with the jump in heating at the nose-IAD juncture is reduced by a minimum of 93%. In this paper, the aeroshell is assumed to be rigid and the inflatable fluid interaction effect is left for future investigations.

Mazaheri, Alireza

Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Johnston, Christopher O.

LAURA Users Manual: 5.4-54166

This users manual provides in-depth information concerning installation and execution of Laura, version 5. Laura is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 Laura code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multi-physics coupling. As a result, Laura now shares gas-physics modules, MPI modules, and other low-level modules with the Fun3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza

LAURA Users Manual: 5.3-48528

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multi-physics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza

LAURA Users Manual: 5.2-43231

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multiphysics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza

Laura Users Manual: 5.1-41601

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multiphysics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza

A Study of Ablation-Flowfield Coupling Relevant to the Orion Heatshield

The coupled interaction between an ablating surface and the surrounding aerothermal environment is studied. An equilibrium ablation model is coupled to the LAURA flowfield solver, which allows the char ablation rate (m-dot(sub c)) to be computed as part of the flowfield solution. The wall temperature (T(sub w)) and pyrolysis ablation rate (m-dot(sub g)) may be specified by the user, obtained from the steady-state ablation approximation, or computed from a a material response code. A 32 species thermochemical nonequilibrium flowfield model is applied, which permits the treatment of C, H, O, N, and Si containing species. Coupled ablation cases relevant to the Orion heatshield are studied. These consist of diffusion-limited oxidation cases with Avcoat as the ablation material. The m-dot(sub c) values predicted from the developed coupled ablation analysis were compared with those obtained from a typical uncoupled ablation analysis. The coupled results were found to be as much as 50% greater than the uncoupled values. This is shown to be a result of the cumulative effect of the two fundamental approximations inherent in the uncoupled analysis.

Johnston, Christopher O.