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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

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↗
NASA NTRS2016

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↗
NASA NTRS2019

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↗
NASA NTRS2019

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↗
NASA NTRS2019

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↗
NASA NTRS2011

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2020

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

Eigenvalue Sensitivity Computations for Linear Stability Theory

To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.

Boundary Layer Transition↗
NASA NTRSDate not supplied

Progress on Quantitative Infrared Thermography at the NASA Langley Aerothermodynamic Laboratory

A quantitative infrared thermography technique is being developed for accurate and high-resolution measurements of surface temperature and heat flux in the hypersonic wind tunnels at NASA Langley Research Center. A hemisphere test campaign was carried out in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel to assess the current performance of the technique and to identify remaining challenges. Measurements were obtained for three different model materials over a range of freestream Reynolds numbers. The raw infrared images were converted to temperature via a radiometric calibration of the camera sensor. Corrections were made for losses due to test article emissivity and window transmissivity. A preliminary uncertainty analysis yielded an estimated accuracy within 4 K at moderate viewing angles, similar to that of a standard thermocouple. Surface heat flux was obtained from the measured surface temperatures via a finite-difference solution of the one-dimensional heat equation that modeled thermal variation of the material properties. Comparisons of the measured heating to both computational and theoretical predictions demonstrated agreement within the expected uncertainty of the technique. The heat transfer measurement near the stagnation point was found to be highly sensitive to the test article emissivity, which establishes the need for improved characterization of the wind tunnel model materials for high temperature applications.

infrared thermography↗
DOE OSTI2023

Enhancing two-photon spontaneous emission in rare earths using graphene and graphene nanoribbons

The enhancement of two-photon spontaneous emission (2PSE) from trivalent and divalent rare earth ions in proximity to graphene and graphene nanoribbons is calculated for achievable experimental conditions using a combination of finite difference time domain simulations and direct computation of transition rates between energy levels in rare earths. For Er 3+ , we find that the 2PSE rate is initially 8 orders lower than the single-photon spontaneous emission rate but that, with enhancement, 2PSE can reach 2.5% of the overall decay. When graphene nanoribbons are used, we also show that the emission of free-space photon pairs from Er 3+ at 3–3.2 μ⁢m via 2PSE can be increased by ~400. Our calculations show significantly less relative graphene-enhanced 2PSE than previous works, and we attribute this variation to differences in emitter size and assumed graphene mobility. Here, we also show that the internal energy structure of the ion can have an impact on the degree of 2PSE enhancement achievable and find that divalent rare earths are more favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗