Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-differences”

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 1,135 records · Page 63

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver, WRLES, was used to compute the flow from a 2-inch round nozzle. Several different flow conditions, consisting of different jet Mach numbers and temperature ratios, were examined. Predictions of mean and fluctuating velocities were compared to previously obtained particle image velocimetry data. Predictions of mean and fluctuating temperature were compared to new data obtained using Raman spectroscopy. Based on the good agreement with experimental data for the individual quantities, the combined quantity turbulent heat flux was examined.

large-eddy simulations↗

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations (LES) were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver was used to compute the flow from a 2-inch round nozzle. Three different flow conditions of varying jet Mach numbers and temperature ratios were examined. The LES results showed that the temperature field behaves similar to the velocity field, but with a more rapidly spreading mixing layer. Predictions of mean, mu-bar(sub i), and fluctuating, mu'(sub i), velocities were compared to particle image velocimetry data. Predictions of mean, T-bar, and fluctuating, T', temperature were compared to data obtained using Rayleigh scattering and Raman spectroscopy. Very good agreement with experimental data was demonstrated for the mean and fluctuating velocities. The LES correctly predicts the behavior of the turbulent temperature field, but over-predicts the levels of the fluctuations. The turbulent heat flux was examined and compared to Reynolds-averaged Navier-Stokes (RANS) results. The LES and RANS simulations produced very similar results for the radial heat flux. However, the axial heat flux obtained from the LES differed significantly from the RANS result in both structure and magnitude, indicating that the gradient diffusion type model in RANS is inadequate. Finally, the LES data was used to compute the turbulent Prandtl number and verify that a constant value of 0.7 used in the RANS models is a reasonable assumption.

hot jets↗

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver, WRLES, was used to compute the flow from a 2-inch round nozzle. Several different flow conditions, consisting of different jet Mach numbers and temperature ratios, were examined. Predictions of mean and fluctuating velocities were compared to previously obtained particle image velocimetry data. Predictions of mean and fluctuating temperature were compared to new data obtained using Raman spectroscopy. Based on the good agreement with experimental data for the individual quantities, the combined quantity turbulent heat flux was examined.

hot jets↗

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

High-Resolution Modeling of the Dust and Water Cycles with the NASA Ames Mars Global Climate Model

NASA’s Mars Climate Modeling Center at Ames Research Center is currently undergoing an exciting period of growth in personnel, modeling capabilities, and science productivity. We are transitioning from our legacy Arakawa C-grid finite-difference dynamical core to the NOAA/GFDL cubed-sphere finite-volume dynamical core for simulating the climate of Mars in a global framework. This highly parallelized core is scalable and flexible, which allows for significant improvements in the horizontal and vertical resolutions of our simulations. We have implemented the Ames water ice cloud microphysics package described in Haberle et al. (2018) into this new dynamical core. We will present high-resolution simulations of the dust and water cycles that show that sub-degree horizontal resolution improves the agreement between the vertical distribution of dust and water ice and observations. In particular, both water ice clouds and dust are transported to higher altitudes due to stronger topographic circulations at high resolution. Preliminary results suggest that high-resolution global modeling is needed to properly capture critical features of the dust and water cycles, and thus the current Mars climate.

Kahre, Melinda A.↗

Jet Noise Prediction Using Hybrid RANS/LES with Structured Overset Grids

Hybrid RANS/LES simulations using the structured overset grid approach and a low dissipation finite-difference method within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for jet noise prediction. The simulations are part of a validation effort to demonstrate jet noise prediction capability, to assess noise characteristics of next generation quiet supersonic commercial jets. Results are compared with experimental data acquired in the Small Hot Jet Acoustic Rig in the Aeroacoustic Propulsion Laboratory at NASA Glenn Research Center. Details of the structured overset grids, numerical discretization, and turbulence model are provided. Near-field comparisons to PIV data and far-field comparisons to microphone data are discussed. Excellent agreement of time-averaged mean quantities along the jet center-line and lip-line are obtained. Good agreement with root-mean-squared (RMS) quantities along the jet center-line are also obtained. An over prediction of lip-line RMS at the nozzle exit is observed leading to a small over-prediction of side-line far-field noise (less than 2 dB). Good agreement with far-field noise in the Mach wave radiation directions is achieved.

Housman, Jeffrey A.↗

Efficient Near-Field to Mid-Field Sonic Boom Propagation Using a High-Order Space Marching Method

An efficient strategy for propagating sonic boom signatures from a near-field Computational Fluid Dynamics (CFD) solution to the mid-field is presented. The method is based on a high-order accurate finite-difference discretization of the 3D Euler equations on a specially designed curvilinear grid and a single sweep space marching solution algorithm. The new approach leads to more than a factor of two reduction in overall computational resources compared to the current method used to propagate near-field sonic booms to the ground. Accuracy and efficiency of the near-field to mid-field process is demonstrated using a selection of test cases from the AIAA Sonic Boom Prediction Workshops. Azimuthal dependence of nonlinear wave propagation from the near-field to mid-field is analyzed along with its effects on the ground level noise.

Housman, Jeffrey A.↗

Does Orbital Angular Momentum Have Effect on Laser’s Scattering by Molecular Atmosphere?

Lasers with orbital angular momentum (OAM) have potential applications in communication technology, manipulation of particles, and remote sensing. Because of its unusual light-scattering properties, the OAM laser’s interaction with a molecular atmosphere must be studied to ensure that it is not lossy for communication or remote-sensing applications that involve its transmission through an atmospheric environment. In this study, the finite-difference time-domain (FDTD) method [21] is applied to calculate the light scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre-Gaussian (LG) beams with OAM by very small dielectric particles. Not like Lorentz-Mie solutions, the FDTD method can calculate for particles off the central axis of the LG beam. It is found that when the particles are very small, and the topological charge number of the OAM of a laser is not extremely large, the laser’s OAM has little effect on the scattering phase function. This suggests that Rayleigh theory can be applied directly to calculate the light scattering by atmospheric molecules. The transmission of a laser beam with OAM in a molecular atmosphere is not different from that of a regular Gaussian beam.

Molecular atmosphere↗

A FDTD Solution of Scattering of Laser Beam with Orbital Angular Momentum by Dielectric Particles: Far-Field Characteristics

Electromagnetic (EM) beams with orbital angular momentum (OAM) may have great potential applications in communication technology and in remote sensing of the Earth-atmosphere system and outer planets. Study of their interaction with optical lenses and dielectric or metallic objects, or scattering of them by particles in the Earth-atmosphere system, is a necessary step to explore the advantage of the OAM EM beams. In this study, the 3-dimensional (3D) scattered-field (SF) finite-difference time domain (FDTD) technique with the convolutional perfectly matched layer (CPML) absorbing boundary conditions (ABC) is applied to calculate the scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre–Gaussian (LG) beams with the OAM by dielectric particles. We found that for OAM beam׳s interaction with dielectric particles, the forward-scattering peak in the conventional phase function (P11) disappears, and light scattering peak occurs at a scattering angle of ~15° to 45°. The disappearance of forward-scattering peak means that, in laser communications most of the particle-scattered noise cannot enter the receiver, thus the received light is optimally the original OAM-encoded signal. This feature of the OAM beam also implies that in lidar remote sensing of the atmospheric particulates, most of the multiple-scattering energy will be off lidar sensors, and this may result in an accurate profiling of particle layers in the atmosphere or in the oceans by lidar, or even in the ground when a ground penetration radar (GPR) with the OAM is applied. This far-field characteristics of the scattered OAM light also imply that the optical theorem, which is derived from plane-parallel wave scattering case and relates the forward scattering amplitude to the total cross section of the scatterer, is invalid for the scattering of OAM beams by dielectric particles.

Electromagnetic beams↗

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

Using Power Diagrams to Build Optimal Unstructured Meshes for C-Grid Models

The Model for Prediction Across Scales (MPAS) for Ocean (-O), Sea-Ice (-SI) and Land-Ice (-LI), in addition to the Coastal Ocean Marine Prediction Across Scales (COMPAS) are two novel general circulation models designed to resolve coupled ocean-ice dynamics over variable spatial scales using non-uniform unstructured grids. Both models are based on a conservative mimetic finite-difference/volume formulation (TRiSK), in which staggered momentum, vorticity and mass-based degrees- of-freedom are distributed over an orthogonal 'primal-dual' mesh.

Engwirda, Darren↗

Adjoint-Based Anisotropic Mesh Adaptation for a Stabilized Finite-Element Flow Solver

An adjoint solver is implemented in the FUN3D stabilized finite-element flow solver. The adjoint solution is used to generate anisotropic, adapted meshes to control error in scalar out- put functionals, such as lift or drag coefficient. The inviscid and turbulent adjoints are verified with a finite-difference derivative approximation and can be used in design optimization in addition to mesh adaptation. The adjoint capability represents an extension of previous research using the FUN3D stabilized finite-element capability with metric-based mesh adaptation and interpolation-based error estimates to generate highly anisotropic adapted meshes for turbulent flows. In the present work, a metric-based approach is again utilized, where the adjoint and the primal solutions both contribute to the generation of a metric tensor field that is subsequently used to produce the required anisotropic mesh for each adaptation cycle. Adaptive results are then shown for an inviscid supersonic flow over a diamond airfoil, inviscid transonic flow over the ONERA M6 wing, and viscous laminar flow over the NACA 0012 airfoil, all using drag force as the output functional. Adjoint-based adaptation is compared with a multiscale solution-based approach that controls the L norm of Mach number interpolation error to demonstrate the effectiveness and effciency of the adjoint-based adaptive mesh technology.

Balan, Aravind↗

Modeling the effects of hydrogen and dose rate sensitivity in CMOS and bipolar technologies

Low dose rate (LDR) and high dose rate (HDR) experiments on field-oxide-field-effect-transistors (FOXFETs) and gated lateral pnp (GLPNP) bipolar transistors indicate that there is a dose rate enhancement factor (EF) associated with radiation-induced degradation. This EF is also affect by the presence of hydrogen in the oxide. In this work, we developed a one dimensional (1-D) numerical calculations code to investigate the key mechanisms that describe the dose rate sensitivity and the effect of hydrogen on dose rate effects. We used a finite-difference methodology for the numerical calculations that allow for computing solutions for the densities of the mobile species as well as for the electrostatic potential at nodes contained within a mesh superimposed on the solution domain. Results from calculations of damage EF indicate that oxide thickness, distribution of hole traps and hole capture cross-section affect dose rate sensitivity. In addition, calculations show that molecular hydrogen cracking at positively charged defects may be a key reaction relating hydrogen and dose rate response. Comparison to experimental data on bipolar and CMOS devices results in good agreement with the dose rate calculations of interface trap buildup

Barnaby, Hugh↗

A Large-Eddy Simulation Of Turbulent Flow Over A Backward Facing Step

Large-eddy simulation (LES) was used to study the turbulent flow over a backward facing step. The simulations examine the details of the turbulent flow’s separated and reattaching regions created by the step geometry. An explicit high-order of accuracy/high resolution finite-difference computational fluid dynamics solver was used for both LES and Reynolds-averaged Navier-Stokes (RANS) simulations. The experimental data from the backward facing step experiment of Driver and Seegmiller was used for comparison. Data presented are derived from a RANS simulation using the Spalart-Allmaras turbulence model and LES simulations using the implicit LES assumption on three different grids of varying resolution. The RANS solution agrees with past work. It does an adequate job of predicting the flowfield and reattachment point, but misses details such as the strength of the recirculating region and magnitudes of the Reynolds stresses. The synthetic eddy method was used to provide a turbulent inflow to the quasi-2D LES domains. For the LES, the coarse grid was of insufficient resolution and under-predicted the skin friction coefficient and over-predicted the normal stresses. The finer grid LES results accurately predict the skin friction coefficient and provide more accurate Reynolds stresses. Triple product velocity correlations are also computed from the LES and the finer grid results did a good job predicting the magnitudes of the curves and the shapes were reasonably represented. It is recommended that future work include investigating the turbulent inflow domain length, width of the quasi-2D domain, and additional grid refinement

turbulence↗

A Large-Eddy Simulation Of Turbulent Flow Over A Backward Facing Step

Large-eddy simulation (LES) was used to study the turbulent flow over a backward facing step. The simulations examine the details of the turbulent flow's separated and reattaching regions created by the step geometry. An explicit high-order of accuracy/high resolution finite-difference computational fluid dynamics solver was used for both LES and Reynolds-averaged Navier-Stokes (RANS) simulations. The experimental data from the backward facing step experiment of Driver and Seegmiller was used for comparison. Data presented are derived from a RANS simulation using the Spalart-Allmaras turbulence model and LES simulations using the implict LES assumption on three different grids of varying resolution. The RANS solution agrees with past work. It does an adequate job of predicting the flowfield and reattachment point, but misses details such as the strength of the recirculating region and magnitudes of the Reynolds stresses. The synthetic eddy method was used to provide a turbulent inflow to the quasi-2D LES domains. For the LES, the coarse grid was of insufficient resolution and under-predicted the skin friction coefficient and over-predicted the normal stresses. The finer grid LES results accurately predict the skin friction coefficient and provide more accurate Reynolds stresses. Triple product velocity correlations are also computed from the LES and the finer grid results did a good job predicting the magnitudes of the curves and the shapes were reasonably represented. It is recommended that future work include investigating the turbulent inflow domain length, width of the quasi-2D domain, and additional grid refinement.

turbulence↗

Broadband Characterization of Anti-reflection Coated Starshade Optical Edges for Solar Glint Control

A starshade is a large flower-shaped screen designed to enable the direct imaging of exoplanets with a space telescope. The starshade perimeter is composed of sharp, precisely shaped edges to minimize the glint of sunlight into the telescope. Past work has focused on bare edges to minimize the terminal radius. This paper describes the broadband, wide-angle performance of edges coated with a thin multi-layer anti-reflection coating. This coating uses a combination of interference and absorption to reduce the surface reflectivity and to avoid the negative effects associated with a large cross-sectional area. A custom scattered light testbed has been developed to quantify the amount of light scattered from sample edges and to validate Finite-Difference Time-Domain (FDTD) models of the optical scatter. We show that optical edge samples with this coating significantly reduce the solar glint pattern compared to similar uncoated optical edges.

Sheikh, David↗

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS↗