Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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,585 records · Page 88

Surface cracks in a plate of finite width under tension or bending

The problem of a finite plate containing collinear surface cracks is considered and solved by using the line spring model with plane elasticity and Reissner's plate theory. The main focus is on the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors in an effort to provide extensive numerical results which may be useful in applications. Some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks, and two corner cracks for wide range of relative dimensions. Particularly in corner cracks, the agreement with the finite element solution is surprisingly very good. The results are obtained for semi-elliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.

Erdogan, F.↗

Nonlinear evolution of the auroral electron beam

The nonlinear spatial evolution, from the source to the atmosphere, of the auroral electron beam and the beam-generated electrostatic whistler noise was studied, calculating changes in beam parameters from equations for the conservation of total particle and wave energy and momentum flux density. Wave power fluxes were calculated by numerically integrating the wave kinetic equations, and the levels of beam-generated noise were determined by using thermal levels of Cerenkov radiation as a source. It was found that beam parameters evolve on ionospheric scale lengths, and their positive slope feature in velocity space is maintained over altitudes measured in thousands of kilometers of altitude, even though they can generate wave energy density fluxes sufficient to modify the ionospheric density profile.

Maggs, James E.↗

Modeling of the Earth's gravity field using the New Global Earth Model (NEWGEM)

Traditionally, the global gravity field was described by representations based on the spherical harmonics (SH) expansion of the geopotential. The SH expansion coefficients were determined by fitting the Earth's gravity data as measured by many different methods including the use of artificial satellites. As gravity data have accumulated with increasingly better accuracies, more of the higher order SH expansion coefficients were determined. The SH representation is useful for describing the gravity field exterior to the Earth but is theoretically invalid on the Earth's surface and in the Earth's interior. A new global Earth model (NEWGEM) (KIM, 1987 and 1988a) was recently proposed to provide a unified description of the Earth's gravity field inside, on, and outside the Earth's surface using the Earth's mass density profile as deduced from seismic studies, elevation and bathymetric information, and local and global gravity data. Using NEWGEM, it is possible to determine the constraints on the mass distribution of the Earth imposed by gravity, topography, and seismic data. NEWGEM is useful in investigating a variety of geophysical phenomena. It is currently being utilized to develop a geophysical interpretation of Kaula's rule. The zeroth order NEWGEM is being used to numerically integrate spherical harmonic expansion coefficients and simultaneously determine the contribution of each layer in the model to a given coefficient. The numerically determined SH expansion coefficients are also being used to test the validity of SH expansions at the surface of the Earth by comparing the resulting SH expansion gravity model with exact calculations of the gravity at the Earth's surface.

Kim, Yeong E.↗

Flight, Wind-Tunnel, and Computational Fluid Dynamics Comparison for Cranked Arrow Wing (F-16XL-1) at Subsonic and Transonic Speeds

Geometrical, flight, computational fluid dynamics (CFD), and wind-tunnel studies for the F-16XL-1 airplane are summarized over a wide range of test conditions. Details are as follows: (1) For geometry, the upper surface of the airplane and the numerical surface description compare reasonably well. (2) For flight, CFD, and wind-tunnel surface pressures, the comparisons are generally good at low angles of attack at both subsonic and transonic speeds, however, local differences are present. In addition, the shock location at transonic speeds from wind-tunnel pressure contours is near the aileron hinge line and generally is in correlative agreement with flight results. (3) For boundary layers, flight profiles were predicted reasonably well for attached flow and underneath the primary vortex but not for the secondary vortex. Flight data indicate the presence of an interaction of the secondary vortex system and the boundary layer and the boundary-layer measurements show the secondary vortex located more outboard than predicted. (4) Predicted and measured skin friction distributions showed qualitative agreement for a two vortex system. (5) Web-based data-extraction and computational-graphical tools have proven useful in expediting the preceding comparisons. (6) Data fusion has produced insightful results for a variety of visualization-based data sets.

John E Lamar↗

Comparison of Experimental Data and Computations Fluid Dynamics Analysis for a Three Dimensional Linear Plug Nozzle

A three dimensional linear plug nozzle of area ratio 12.79 was designed by EADS Space Transportation (former Astrium Space Infrastructure). The nozzle was tested within the German National Technology Program 'LION' in a cold air wind tunnel by TU Dresden. The experimental hardware and test conditions are described. Experimental data was obtained for the nozzle without plug side wall fences at a nozzle pressure ratio of 116 and then with plug side wall fences at NPR 110. Schlieren images were recorded and axial profiles of plug wall static pressures were measured at several spanwise locations and on the plug base. Detailed CFD analysis was performed for these nozzle configurations at NPR 116 by NASA MSFC. The CFD exhibits good agreement with the experimental data. A detailed comparison of the CFD results and the experimental plug wall pressure data are given. Comparisons are made for both the without and with plug side wall fence configurations. Numerical results for density gradient are compared to experimental Schlieren images. Experimental nozzle thrust efficiencies are calculated based on the CFD results. The CFD results are used to illustrate the plug nozzle fluid dynamics. The effect of the plug side wall is emphasized.

Ruf, J. H.↗

Comparison of Experimental Data and Computations Fluid Dynamics Analysis for a Three Dimensional Linear Plug Nozzle

A three dimensional linear plug nozzle of area ratio 12.79 was designed by Astrium. The nozzle was tested within the German National Technology Program LION in a cold air wind tunnel by TU Dresden. The experimental hardware and test conditions are described. Experimental data was obtained for the nozzle without plug side wall fences and then with plug side wall fences. Experimental data for two nozzle pressure ratios (NPR), 116 and 45, are presented for the without fence and with fence configurations. Schlieren images of both NPR were recorded. Axial profiles of plug wall static pressures were measured at several spanwise locations and on the plug base. Detailed computational fluid dynamics (CFD) analysis was performed for these nozzle configurations by NASA MSFC. The CFD exhibits good agreement with the experimental data. A detailed comparison of the CFD results and the experimental plug wall pressure data is given for four test conditions; at both NPRs, without and with plug side wall fences. Numerical schlieren images are compared to experimental schlieren images. Nozzle thrust efficiencies are calculated from the CFD results. The CFD results are used to illustrate the plug nozzle fluid dynamics for all four test conditions. The effect of the plug side wall fences at both NPRs is emphasized.

Ruf, J. H.↗

The Atmospheric Infrared Sounder- An Overview

The Atmospheric Infrared Sounder (AIRS) was launched in May 2002. Along with two companion microwave sensors, it forms the AIRS Sounding Suite. This system is the most advanced atmospheric sounding system to date, with measurement accuracies far surpassing those available on current weather satellites. The data products are calibrated radiances from all three sensors and a number of derived geophysical parameters, including vertical temperature and humidity profiles, surface temperature, cloud fraction, cIoud top pressure, and profiles of ozone. These products are generated under cloudy as well as clear conditions. An ongoing calibration validation effort has confirmed that the system is very accurate and stable, and many of the geophysical parameters have been validated. AIRS is in some cases more accurate than any other source and can therefore be difficult to validate, but this offers interesting new research opportunities. The applications for the AIRS products range from numerical weather prediction to atmospheric research - where the AIRS water vapor products near the surface and in the mid to upper troposphere will make it possible to characterize and model phenomena that are key for short-term atmospheric processes, such as weather patterns, to long-term processes, such as interannual cycles (e.g., El Nino) and climate change.

microwave↗

High-Speed Edge-Detecting Line Scan Smart Camera

A high-speed edge-detecting line scan smart camera was developed. The camera is designed to operate as a component in a NASA Glenn Research Center developed inlet shock detection system. The inlet shock is detected by projecting a laser sheet through the airflow. The shock within the airflow is the densest part and refracts the laser sheet the most in its vicinity, leaving a dark spot or shadowgraph. These spots show up as a dip or negative peak within the pixel intensity profile of an image of the projected laser sheet. The smart camera acquires and processes in real-time the linear image containing the shock shadowgraph and outputting the shock location. Previously a high-speed camera and personal computer would perform the image capture and processing to determine the shock location. This innovation consists of a linear image sensor, analog signal processing circuit, and a digital circuit that provides a numerical digital output of the shock or negative edge location. The smart camera is capable of capturing and processing linear images at over 1,000 frames per second. The edges are identified as numeric pixel values within the linear array of pixels, and the edge location information can be sent out from the circuit in a variety of ways, such as by using a microcontroller and onboard or external digital interface to include serial data such as RS-232/485, USB, Ethernet, or CAN BUS; parallel digital data; or an analog signal. The smart camera system can be integrated into a small package with a relatively small number of parts, reducing size and increasing reliability over the previous imaging system..

Prokop, Norman F.↗

A Twenty-Year Analysis of Winds in California for Offshore Wind Energy Production Using WRF v4.1.2

Offshore wind resource characterization in the United States relies heavily on simulated winds from numerical weather prediction (NWP) models, given the lack of hub-height observations offshore. One such NWP data set used extensively by U.S. stakeholders is the Wind Integration National Dataset (WIND) Toolkit, a 7-year time-series data set produced in 2013 by the National Renewable Energy Laboratory. In this study, we present an update to that data set for offshore California that leverages recent advancements in NWP modeling capabilities and extends the period of record to a full 20 years. The data set predicts a significantly larger wind resource (0.25–1.75 m s-1 stronger), including in three Call Areas that the Bureau of Ocean Energy Management is considering for commercial activity. We conduct a set of yearlong simulations to study factors that contribute to this increase in the modeled wind resource. The largest impact arises from a change in the planetary boundary layer parameterization from the Yonsei University scheme to the Mellor-Yamada-Nakanishi-Niino scheme and their diverging wind profiles under stable stratification. Additionally, we conduct a refined wind resource assessment at the three Call Areas, characterizing distributions of wind speed, shear, veer, stability, frequency of wind droughts, and power production. We find that, depending on the attribute, the new data set can show substantial disagreement with the WIND Toolkit, thereby driving important changes in predicted power.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Evaluation of the Rapid Refresh Numerical Weather Prediction Model over Arctic Alaska

Abstract Despite a need for accurate weather forecasts for societal and economic interests in the U.S. Arctic, thorough evaluations of operational numerical weather prediction in the region have been limited. In particular, the Rapid Refresh Model (RAP), which plays a key role in short-term forecasting and decision-making, has seen very limited assessment in northern Alaska, with most evaluation efforts focused on lower latitudes. In the present study, we verify forecasts from version 4 of the RAP against radiosonde, surface meteorological, and radiative flux observations from two Arctic sites on the northern Alaskan coastline, with a focus on boundary layer thermodynamic and dynamic biases, model representation of surface inversions, and cloud characteristics. We find persistent seasonal thermodynamic biases near the surface that vary with wind direction, and may be related to the RAP’s handling of sea ice and ocean interactions. These biases seem to have diminished in the latest version of the RAP (version 5), which includes refined handling of sea ice, among other improvements. In addition, we find that despite capturing boundary layer temperature profiles well overall, the RAP struggles to consistently represent strong, shallow surface inversions. Further, while the RAP seems to forecast the presence of clouds accurately in most cases, there are errors in the simulated characteristics of these clouds, which we hypothesize may be related to the RAP’s treatment of mixed-phase clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dilution jet mixing program, supplementary report

The velocity and temperature distributions predicted by a 3-D numerical model and experimental measurements are compared. Empirical correlations for the jet velocity trajectory developed are presented. The measured velocity distributions for all test cases of phase through phase 3 are presented in the form of contour and oblique plots. quantification of the effects of the following on the jet mixing characteristics with a confined crossflow are: (1) orifice geometry momentum flux ratio and density ratio; (2) nonuniform mainstream temperature and velocity profiles upstream of dilution orifices; (3) cold versus hot jet injection; (4) cross-stream flow are a convergence as encountered in practical dilution zone geometries; (5) 2-D slot versus circular orifices; (6) discrete noncirculcer orifices; (7) single-sided versus opposed jets; (8) single row of jets.

Srinivasan, R.↗

Remote sensing of land-surface temperature from HIRS/MSU data

A relaxation algorithm which permits meteorological parameters to be obtained from satellite data, without a priori assumptions about the properties of the other unknowns in the field of view, was developed. Atmospheric temperature profiles, atmospheric humidity, cloud cover, cloud top height, cloud top temperature, sea-surface temperature, land-surface temperature, snow cover, and ice cover are derived. Simultaneous determination of atmospheric and surface thermal structure and the cloud distribution provides information on heat sources and sinks, storage rates, and transport phenomena in the atmosphere. Such information is critical in determining the driving mechanisms for motions in the atmosphere and oceans and in improving numerical weather prediction.

Chahine, M. T.↗

Quasilinear saturation of forced current sheet tearing modes

Numerical studies of tearing modes in a nearly singular forced current sheet equilibrium (Liewer and Payne, 1990) show that the modes saturate quasilinearly when the width of the magnetic island formed by the reconnection is on the order of several times the linear mode width which scales as approximately (kS) exp -2/5, where S is the Lundquist number and k is the wavenumber. The modes saturate quasilinearly by flattening the current profile, converting magnetic energy into plasma energy. The longer wavelength modes, which saturate at higher levels, release the most energy. These modes may, nonlinearly, play a role in coronal heating when sharp current sheets form as a result of global magnetic stresses.

Liewer, Paulett C.↗

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗

SHarD: A beam dynamics simulation code for dielectric laser accelerators based on spatial harmonic field expansion

In order to demonstrate acceleration of electrons to relativistic scales by an on chip dielectric laser accelerator (DLA), a ponderomotive focusing scheme capable of capturing and transporting electrons through nanometer-scale apertures over extended interaction lengths has been proposed. Here we present a Matlab-based numerical code (SHarD) utilizing a spatial harmonic expansion of the fields within the dielectric structure to simulate the evolution of the beam phase space distribution in this scheme. The code can be used to optimize key-parameters for the accelerator performance such as the final energy, transverse spot size evolution and total number of electrons accelerated through currently fabricated structures. Eventually, the simulation model will be applied to inform the phase mask profile to be added to a pulse front tilt drive laser pulse using a liquid crystal mask in the experimental setup being assembled at UCLA Pegasus Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fingerprinting shock-induced deformations via diffraction

Abstract During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

36 MATERIALS SCIENCE↗

Modelling the limiter ramp-up of WEST for addressing the future challenges of ITER

This paper presents a joint experimental and numerical investigation into the physics of long limited plasma ramp-up in tokamaks with tungsten (W) first walls, a critical phase for ITER operation. The comparison between the average plasma quantities simulated using the SolEdge-HDG code and the measurements taken during three successive WEST discharges after boronisation shows how challenging it is to predict this phase. While simulations reproduce the general trends at the midplane, with reasonable match in density profiles, they consistently underestimate core temperatures possibly due to too large perpendicular heat conductivity. On the contrary, at the high field side (HFS) limiter, simulations overestimate the measured quantities, and highlights the limitation of using Bohm boundary conditions at grazing magnetic angles. Experimental measurements reveal that the boron layer is rapidly eroded, on a timescale comparable to a single ITER discharge. The subsequent transition from a boron-coated to a tungsten wall increases recycling and significantly degrades the core plasma, reducing the electron temperature by nearly half due to W contamination, despite wall parameters remaining stable. Furthermore, comparisons with Langmuir probes, bolometry, reflectometry, and spectroscopy indicate that the experimental far scrape-Off layer (SOL) is significantly wider than simulated. This wide SOL implies that boron erosion extends along the entire HFS limiter rather than being confined to the contact point. This work highlights some characteristics of the plasma during this phase of the discharge and emphasizes the current modelling issues that need to be resolved in order to obtain reliable predictions concerning the ITER ramp-up.

ITER↗