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,189 records · Page 66

Acoustic propagation in a thermally stratified atmosphere

This report describes the activities during the fourth six month period of the investigation of acoustic propagation in the atmosphere with a realistic lapse temperature profile. A significant error was detected since the previous semi-annual report and has been corrected in both the plane wave and point source solutions. This report then describes both of these problems in some detail along with presenting some numerical results from the model. Work will begin this summer on the model of propagation in an inversion.

Vanmoorhem, W. K.↗

Double-diffusive convection with sidewalls

Stommel et al. (1956) have first described an instability, known as thermosolutal convection, thermohaline convection, or double-diffusive convection. This instability may occur in the case of a fluid in a gravitational field with two diffusing components present. The present study is concerned with the effect of sidewalls on flow in the fingering regime in the absence of applied horizontal gradients. The work was motivated by numerical results obtained on the basis of a simulation of thermosolutal convection occurring during the unidirectional solidification of a binary alloy. In this case, the unperturbed solute field in the liquid ahead of the solidifying planar interface has an exponential vertical profile because of the rejection or preferential incorporation of solute by the solid phase.

Mcfadden, G. B.↗

An application of active surface heating for augmenting lift and reducing drag of an airfoil

Application of active control to separated flow on the RC(6)-08 airfoil at high angle of attack by localized surface heating is numerically simulated by integrating the compressible 2-D nonlinear Navier-Stokes equation solver. Active control is simulated by local modification of the temperature boundary condition over a narrow strip of the upper surface of the airfoil. Both mean and perturbed profiles are favorably altered when excited with the same natural frequency of the shear layer by moderate surface heating for both laminar and turbulent separation. The shear layer is found to be very sensitive to localized surface heating in the vicinity of the separation point. The excitation field at the surface sufficiently altered both the local as well as the global circulation to cause a significant increase in lift and reduction in drag.

Maestrello, Lucio↗

Thermal calibration and analysis of a Bridgman flight furnace

A combination of computer modeling and experimental measurements is proposed which permits the calibration of the temperature profile for a spacecraft flight furnace. The Fluent code based on the finite difference method uses the furnace set-point temperatures as boundary conditions to model a Bridgman directional solidification furnace. Experimental measurements are employed as a comparison for the numerical analysis, and the radial gradients and wall temperatures are selected when the calculated values are sufficiently close to those from measurements. The calculations are found to agree with experimental data, and some results of the modeling are given which relate to specific furnace parameters. The CFD approach facilitates the transposition of furnace calibration to microgravitational operation in addition to providing precise initial radial temperature gradients and furnace-wall temperatures.

Knuteson, D. J.↗

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗

COSIM: A Finite-Difference Computer Model to Predict Ternary Concentration Profiles Associated with Oxidation and Interdiffusion of Overlay-Coated Substrates

A finite-difference computer program (COSIM) has been written which models the one-dimensional, diffusional transport associated with high-temperature oxidation and interdiffusion of overlay-coated substrates. The program predicts concentration profiles for up to three elements in the coating and substrate after various oxidation exposures. Surface recession due to solute loss is also predicted. Ternary cross terms and concentration-dependent diffusion coefficients are taken into account. The program also incorporates a previously-developed oxide growth and spalling model to simulate either isothermal or cyclic oxidation exposures. In addition to predicting concentration profiles after various oxidation exposures, the program can also be used to predict coating fife based on a concentration dependent failure criterion (e.g., surface solute content drops to two percent). The computer code, written in an extension of FORTRAN 77, employs numerous subroutines to make the program flexible and easily modifiable to other coating oxidation problems.

Nesbitt, James A.↗

Numerical Issues for Circulation Control Calculations

Steady-state and time-accurate two-dimensional solutions of the compressible Reynolds-averaged Navier- Stokes equations are obtained for flow over the Lockheed circulation control (CC) airfoil and the General Aviation CC (GACC) airfoil. Numerical issues in computing circulation control flows such as the effects of grid resolution, boundary and initial conditions, and unsteadiness are addressed. For the Lockheed CC airfoil computed solutions are compared with detailed experimental data, which include velocity and Reynolds stress profiles. Three turbulence models, having either one or two transport equations, are considered. Solutions are obtained on a sequence of meshes, with mesh refinement primarily concentrated on the airfoil circular trailing edge. Several effects related to mesh refinement are identified. For example, sometimes sufficient mesh resolution can exclude nonphysical solutions, which can occur in CC airfoil calculations. Also, sensitivities of the turbulence models with mesh refinement are discussed. In the case of the GACC airfoil the focus is on the difference between steady-state and time-accurate solutions. A specific objective is to determine if there is self-excited vortex shedding from the jet slot lip.

Swanson, Roy C., Jr.↗

Advanced Intensity-Modulation Continuous-Wave Lidar Techniques for ASCENDS O2 Column Measurements

Global atmospheric carbon dioxide (CO2) measurements for the NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) space mission are critical for improving our understanding of global CO2 sources and sinks. Advanced Intensity- Modulated Continuous-Wave (IM-CW) lidar techniques are investigated as a means of facilitating CO2 measurements from space to meet the ASCENDS measurement requirements. In recent numerical, laboratory and flight experiments we have successfully used the Binary Phase Shift Keying (BPSK) modulation technique to uniquely discriminate surface lidar returns from intermediate aerosol and cloud contamination. We demonstrate the utility of BPSK to eliminate sidelobes in the range profile as a means of making Integrated Path Differential Absorption (IPDA) column CO2 measurements in the presence of optically thin clouds, thereby eliminating the need to correct for sidelobe bias errors caused by the clouds. Furthermore, high accuracy and precision ranging to the surface as well as to the top of intermediate cloud layers, which is a requirement for the inversion of column CO2 number density measurements to column CO2 mixing ratios, has been demonstrated using new hyperfine interpolation techniques that takes advantage of the periodicity of the modulation waveforms. This approach works well for both BPSK and linear swept-frequency modulation techniques. The BPSK technique under investigation has excellent auto-correlation properties while possessing a finite bandwidth. A comparison of BPSK and linear swept-frequency is also discussed in this paper. These results are extended to include Richardson-Lucy deconvolution techniques to extend the resolution of the lidar beyond that implied by limit of the bandwidth of the modulation, where it is shown useful for making tree canopy measurements.

Campbell, Joel F.↗

Assimilation of Microwave Cloudy Observations over the Rainband of Hurricanes Using a Novel Bayesian Monte Carlo Technique

We propose a novel Bayesian Monte Carlo Integration (BMCI) technique to retrieve the profiles of temperature, water vapor, and cloud liquid/ice water content from microwave cloudy measurements in the rainbands of tropical cyclones (TC). These retrievals then can either be directly used by meteorologists to analyze the structure of TCs or be assimilated into numerical models to provide accurate initial conditions for the NWP models. The BMCI technique is applied to the data from the Advanced Technology Microwave Sounder (ATMS) onboard Suomi National Polar-orbiting Partnership (NPP) and Global Precipitation Measurement (GPM) Microwave Imager (GMI).

Moradi, Isaac↗

Global NWP Impacts of Infrared Sounders from Geostationary Orbit

The Geostationary eXtended Observations (GeoXO) program, expected to launch in the 2030s, includes a proposed infrared (IR) sounder that would provide persistent atmospheric profile observations over much of the western hemisphere. In preparation, the National Aeronautics and Space Administration (NASA) Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework was used to assess the impact of such an instrument on global numerical weather prediction (NWP) as part of a global “ring” of such instruments. Building on previous preliminary studies, an evaluation of the impact of geostationary IR sounders will be presented with foci on the analysis, forecasts, and the forecast sensitivity observation impact (FSOI) metric, including an examination of the consequences for tropical cyclone representation. Overall, assimilation of geostationary IR provide a beneficial impact for NWP applications.

Erica L. Mcgrath-spangler↗

Numerical simulation of steady supersonic flow over an ogive-cylinder-boattail body

A recently reported parabolized Navier-Stokes code has been employed to compute the supersonic flow field surrounding an ogive-cylinder-boattail body at incidence. The computations were performed for flow conditions where an extensive series of experimental surface pressure and turbulent boundary-layer profile measurements had been obtained. Comparison between the computational results and experimental measurements for angles of attack up to 6 deg show excellent agreement. At angles greater than 6 deg discrepancies are observed which are tentatively attributed to three-dimensional turbulence modeling errors.

Schiff, L. B.↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GFCCLib: Scalable and efficient coupled-cluster Green's function library for accurately tackling many-body electronic structure problems

Coupled-cluster Green’s function (GFCC) calculation has drawn much attention in the recent years for targeting the molecular and material electronic structure problems from a many-body perspective in a systematically improvable way. However, GFCC calculations on scientific computing clusters usually suffer from expensive higher di- mensional tensor contractions in the complex space, expensive inter-process communi- cation, and severe load imbalance, which limits it’s routine use for tackling electronic structure problems. Here we present a numerical library prototype that is specifically designed for large-scale GFCC calculations. The design of the library is focused on a systematically optimal computing strategy to improve its scalability and efficiency. The performance of the library is demonstrated by the relevant profiling analysis of running GFCC calculations on remote giant computing clusters. The capability of the library is highlighted by computing a wide near valence band of a fullerene C60 molecule for the first time at the GFCCSD level that shows excellent agreement with the experimental spectrum.

Peng, Bo↗

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗

Quantifying the thermal effect and methyl radical production in nanosecond repetitively pulsed glow discharges applied to a methane-air flame

In this work, we investigated non-equilibrium plasma produced by nanosecond repetitively pulsed glow discharges applied across a lean premixed methane-air flame. The flame is stationary, axisymmetric, and laminar. The discharges are applied on the symmetry axis crossing the reactant gases, flame front, and product gases, allowing phase-locked averaged measurements and comparisons with axisymmetric numerical simulations. The thermal effect and methyl radical production are quantified in the discharge in the reactant gas region. One-dimensional, two-beam, hybrid, femtosecond-picosecond, coherent anti-Stokes Raman scattering is used to acquire spatial and temporal profiles of temperature and oxygen-to-nitrogen concentration ratio. Photo-fragmentation laser-induced fluorescence is used to acquire quantitative two-dimensional profiles of methyl radicals in the discharge providing the first quantitative imaging of methyl produced ahead of a flame by plasma-induced methane dissociation. The spatial profiles of temperature and oxygen-to-nitrogen concentration ratio are in steady state, indicating that individual discharges have an insignificant heating effect. Upper and lower bounds of the produced mole fraction of methyl radicals in the plasma are obtained due to uncertainties in the collisional quenching rates of excited state methylidyne radicals in the plasma. The discharges produce a maximum of 600–1100 ppm of methyl radicals upstream of the flame front within 25 ns. This amount is similar to the predicted methyl mole fraction for the flame without plasma and thus represents a significant chemical perturbation to the reactants upstream of the flame front. The produced methyl follows an exponential decay in the first microsecond after the discharge with a decay constant of 8 µs close to the flame, and 0.8 µs further from the flame. The decay then deviates from the exponential curve and the methyl persists for tens of microseconds. The results suggest that for the tested configuration, the thermal effect of individual discharges through fast gas heating is negligible, while active chemical species are produced in large quantities in the reactant gases, upstream of the flame front.

hybrid fs-ps CARS↗

Characterizing extreme laser intensities by ponderomotive acceleration of protons from rarified gas

Abstract A new method to diagnose extreme laser intensities through measurement of angular and spectral distributions of protons directly accelerated by the laser focused into a rarefied gas is proposed. We simulated a laser pulse focused by an off-axis parabolic mirror by Stratton–Chu integrals, that enables description of laser pulse with different spatial-temporal profiles focusing in a focal spot down to the diffraction limit, that makes our theoretical predictions be a basis for experimental realization. The relationship between characteristics of the proton distributions and parameters of the laser pulse have been analyzed. The analytical and numerical results obtained justify the new method of laser diagnostics. The proposed scheme should be valuable for the commissioning of new extreme intensity laser facilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗