Engineering topics

Mansour, Nagi N.

Publications and source records attributed to Mansour, Nagi N..

At least 73 records · Page 4

NASA NTRS2013

Approach to Integrate Global-Sun Models of Magnetic Flux Emergence and Transport for Space Weather Studies

The Sun lies at the center of space weather and is the source of its variability. The primary input to coronal and solar wind models is the activity of the magnetic field in the solar photosphere. Recent advancements in solar observations and numerical simulations provide a basis for developing physics-based models for the dynamics of the magnetic field from the deep convection zone of the Sun to the corona with the goal of providing robust near real-time boundary conditions at the base of space weather forecast models. The goal is to develop new strategic capabilities that enable characterization and prediction of the magnetic field structure and flow dynamics of the Sun by assimilating data from helioseismology and magnetic field observations into physics-based realistic magnetohydrodynamics (MHD) simulations. The integration of first-principle modeling of solar magnetism and flow dynamics with real-time observational data via advanced data assimilation methods is a new, transformative step in space weather research and prediction. This approach will substantially enhance an existing model of magnetic flux distribution and transport developed by the Air Force Research Lab. The development plan is to use the Space Weather Modeling Framework (SWMF) to develop Coupled Models for Emerging flux Simulations (CMES) that couples three existing models: (1) an MHD formulation with the anelastic approximation to simulate the deep convection zone (FSAM code), (2) an MHD formulation with full compressible Navier-Stokes equations and a detailed description of radiative transfer and thermodynamics to simulate near-surface convection and the photosphere (Stagger code), and (3) an MHD formulation with full, compressible Navier-Stokes equations and an approximate description of radiative transfer and heating to simulate the corona (Module in BATS-R-US). CMES will enable simulations of the emergence of magnetic structures from the deep convection zone to the corona. Finally, a plan will be summarized on the development of a Flux Emergence Prediction Tool (FEPT) in which helioseismology-derived data and vector magnetic maps are assimilated into CMES that couples the dynamics of magnetic flux from the deep interior to the corona.

Mansour, Nagi N.↗

NASA NTRS2012

Performance of Reynolds Averaged Navier-Stokes Models in Predicting Separated Flows: Study of the Hump Flow Model Problem

Separation can be seen in most aerodynamic flows, but accurate prediction of separated flows is still a challenging problem for computational fluid dynamics (CFD) tools. The behavior of several Reynolds Averaged Navier-Stokes (RANS) models in predicting the separated ow over a wall-mounted hump is studied. The strengths and weaknesses of the most popular RANS models (Spalart-Allmaras, k-epsilon, k-omega, k-omega-SST) are evaluated using the open source software OpenFOAM. The hump ow modeled in this work has been documented in the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control. Only the baseline case is treated; the slot flow control cases are not considered in this paper. Particular attention is given to predicting the size of the recirculation bubble, the position of the reattachment point, and the velocity profiles downstream of the hump.

Cappelli, Daniele↗

NASA NTRS2011

A Short Review of Ablative-Material Response Models and Simulation Tools

A review of the governing equations and boundary conditions used to model the response of ablative materials submitted to a high-enthalpy flow is proposed. The heritage of model-development efforts undertaken in the 1960s is extremely clear: the bases of the models used in the community are mathematically equivalent. Most of the material-response codes implement a single model in which the equation parameters may be modified to model different materials or conditions. The level of fidelity of the models implemented in design tools only slightly varies. Research and development codes are generally more advanced but often not as robust. The capabilities of each of these codes are summarized in a color-coded table along with research and development efforts currently in progress.

Lachaud, Jean↗

NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

Ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic matrix is modeled at microscopic scale. Direct numerical simulations show that the carbonized phenolic matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolic is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000 K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean Romain↗

NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

This slide presentation reviews the ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic-formaldehyde matrix is modeled at microscopic scale. Direct numerical simulations show that the matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolics is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean↗

NASA NTRS2008

Multiscale Modeling of Ablation and Pyrolysis in PICA-Like materials

During atmospheric entry of planetary probes, the thermal protection system (TIPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous fibrous TIPS materials may undergo oxidation leading to mass loss and wall recession called ablation. This work aims to improve the understanding of material/environment interactions through a study of the coupling between oxygen transport in the Knudsen regime, heterogeneous oxidation of carbon, and surface recession. A 3D Random Walk Monte Carlo simulation tool is used for this study. The fibrous architecture of a model material, consisting of high porosity random array of carbon fibers, is numerically represented on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession of the fibers is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The model is applied to Stardust mission reentry conditions and explains the unexpected behavior of the TIPS material that underwent mass loss in volume.

Lachaud, Jean↗

NASA NTRS2008

Modeling Ablation of Fibrous Materials from Bulk to Knudsen Regime

Material-environment interactions are analyzed at microscopic scale to explain the lower than expected density observed by post-flight analysis of the char layer on the Stardust shield. Mass transfer, ablation (oxidation), and surface recession of fibrous material is simulated in 3D using a Monte-Carlo simulation tool. Ablation is found to occur either at the surface or in volume depending on Knudsen and Thiele number values. This study supports the idea of volume ablation followed by possible carbon fiber spallation that may explain post-flight analyses.

Lachaud, Jean↗

NASA NTRS2004

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗

NASA NTRS2004

Annual Research Briefs, 2004: Center for Turbulence Research

This report contains the 2004 annual progress reports of the Research Fellows and students of the Center for Turbulence Research in its eighteenth year of operation. Since its inception in 1987, the objective of the CTR has been to advance the physical understanding of turbulent flows and development of physics based predictive tools for engineering analysis and turbulence control. Turbulence is ubiquitous in nature and in engineering devices. The studies at CTR have been motivated by applications where turbulence effects are significant; these include a broad range of technical areas such as planetary boundary layers, formation of planets, solar convection, magnetohydrodynamics, environmental and eco systems, aerodynamic noise, propulsion systems and high speed transportation. Numerical simulation has been the predominant research tool at CTR which has required a critical mass of researchers in numerical analysis and computer science in addition to core disciplines such as applied mathematics, chemical kinetics and fluid mechanics. Maintaining and promoting this interdisciplinary culture has been a hallmark of CTR and has been responsible for the realization of the results of its basic research in applications. The first group of reports in this volume are directed towards development, analysis and application of novel numerical methods for ow simulations. Development of methods for large eddy simulation of complex flows has been a central theme in this group. The second group is concerned with turbulent combustion, scalar transport and multi-phase ows. The nal group is devoted to geophysical turbulence where the problem of solar convection has been a new focus of considerable attention recently at CTR.

Moin, Parviz↗

NASA NTRS2002

Annual Research Briefs, 2002: Center for Turbulence Research

Turbulent combustion remains the largest component of the CTR's core program. This program and several related activities at CTR are supported by NASA's Ultra Efficient Engine Technology Program. It is also intimately connected with the Department of Energy's ASCI program at Stanford which develops the technology for numerical simulation of realistic aircraft engines using state of the art massively parallel computers. In combustion modeling the attention has been directed to the modeling of higher levels of complexity such as spray dynamics, radiation and soot formation. Major aircraft engine manufacturers have shown considerable interest in this program; in particular, a significant active collaboration exists between CTR and the Pratt & Whitney Corporation. CTR's combustion program is essentially based on the large-eddy simulation technique, LES, which is actively being pursued at CTR for this and many other applications. Important accomplishments in LES included simulations with three-dimensional filters which result in grid independent calculations (that is why we call it "true" LES), and the development of the methodology for integration of LES and Reynolds Averaged computations. Optimization techniques are being studied and used for the important problem of wall boundary conditions for LES as well as for optimal shape design for aeroacoustic and aerodynamic performance gains.

ULTRA EFFICIENT ENGINE TECHNOLOGY↗

NASA NTRS2002

The Speed of Axial Propagation of a Cylindrical Bubble Through a Cylindrical Vortex

Inspired by the rapid elongation of air columns injected into vortices by dolphins, we present an exact inviscid solution for the axial speed (assumed steady) of propagation of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The bubble is assumed to be held at constant pressure by being connected to a reservoir, the lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, the propagation speed of the bubble (relative to axial velocity within the vortex) varies between 0.5 and 0.6 of the maximum rotational speed of the vortex. Surprisingly, the bubble tip can propagate (almost as rapidly) even when the pressure minimum in the vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple on the nose of the bubble. A situation important for incipient vortex cavitation, and one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose internal pressure may vary. Under the assumption that the acceleration term is small (checked a posteriori), the steady solution is applied at each instant during the elongation. Three types of behavior are then possible depending on physical parameters and initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and nearly constant volume. Volume begins to decrease in the late stages. (B1) Elongation with decreasing bubble pressure. A limit point of the steady solution is encountered at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure and indefinite creation of volume. This is made possible by the existence of propagating solutions at bubble pressures below the minimum vortex pressure. As the bubble stretches, its radius initially decreases but then becomes constant; this is also observed in experiments on incipient vortex cavitation.

Shariff, Karim↗

NASA NTRS2002

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This talk describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗

NASA NTRS2001

B-spline Method in Fluid Dynamics

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Botella, Olivier↗

NASA NTRS2001

Energy Stable Flux Formulas For The Discontinuous Galerkin Discretization Of First Order Nonlinear Conservation Laws

We consider the discontinuous Galerkin (DG) finite element discretization of first order systems of conservation laws derivable as moments of the kinetic Boltzmann equation. This includes well known conservation law systems such as the Euler For the class of first order nonlinear conservation laws equipped with an entropy extension, an energy analysis of the DG method for the Cauchy initial value problem is developed. Using this DG energy analysis, several new variants of existing numerical flux functions are derived and shown to be energy stable.

Barth, Timothy↗

NASA NTRS2001

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele↗

NASA NTRS2001

Large Eddy Simulation of a Turbulent Jet

Here we present the results of a Large Eddy Simulation of a non-buoyant jet issuing from a circular orifice in a wall, and developing in neutral surroundings. The effects of the subgrid scales on the large eddies have been modeled with the dynamic large eddy simulation model applied to the fully 3D domain in spherical coordinates. The simulation captures the unsteady motions of the large-scales within the jet as well as the laminar motions in the entrainment region surrounding the jet. The computed time-averaged statistics (mean velocity, concentration, and turbulence parameters) compare well with laboratory data without invoking an empirical entrainment coefficient as employed by line integral models. The use of the large eddy simulation technique allows examination of unsteady and inhomogeneous features such as the evolution of eddies and the details of the entrainment process.

Webb, A. T.↗

NASA NTRS2001

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗

NASA NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗