Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

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 991 records · Page 55

Data Democratization: Challenges and Opportunities

Democratizing Earth data is one of the challenges many organizations around the world face in order to maximize the use of their Earth data for research, applications, education, and societal benefits. For example, at the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC), over 1600 global and regional datasets in several NASA Earth science focus areas, including atmospheric composition, water and energy cycles, and climate variability, are archived and distributed to the public. Giovanni, the Geospatial Interactive Online Visualization and Analysis Infrastructure, was developed by GES DISC to facilitate data access and exploration, especially for novice users of Earth science. With Giovanni, users can analyze and visualize over 2000 Earth science variables (e.g., precipitation, aerosol, surface wind) without downloading data, software, the expert understanding of data formats and structures, and coding skills, lowering the barrier to data analysis/comparison by preprocessing and accessing to the data. Results of data analysis and visualization can be accessed in several popular formats (e.g., NetCDF, CSV). As a result of Giovanni's efforts, more than 3000 referral papers have been published in various fields. In spite of this, Giovanni is still difficult to use for some users. For instance, if one searches for "precipitation," it will return over 150 related variables. The question is, which one to use? Furthermore, variables from different data providers (e.g., satellites and models) are named differently with different units, further confusing users, especially those outside the communities. Data democratization is complex and multifaceted. Challenges include service and data discovery, user experiences, visualization, data quality, trustworthiness, and more. In this presentation, we will examine Giovanni as an example of challenges and opportunities in developing data democratization services.

data democratization↗

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗

Influence of Ice Particle Surface Roughening on the Global Cloud Radiative Effect

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.

clouds↗

Modern Experimental Techniques in Turbine Engine Testing

The paper describes application of two modern experimental techniques, thin-film thermocouples and pressure sensitive paint, to measurement in turbine engine components. A growing trend of using computational codes in turbomachinery design and development requires experimental techniques to refocus from overall performance testing to acquisition of detailed data on flow and heat transfer physics to validate these codes for design applications. The discussed experimental techniques satisfy this shift in focus. Both techniques are nonintrusive in practical terms. The thin-film thermocouple technique improves accuracy of surface temperature and heat transfer measurements. The pressure sensitive paint technique supplies areal surface pressure data rather than discrete point values only. The paper summarizes our experience with these techniques and suggests improvements to ease the application of these techniques for future turbomachinery research and code verifications.

Lepicovsky, J.↗

Investigations of radiation-induced and carrier-enhanced conductivity

A steady-state carrier computer code, PECK (Parker Enhanced Carrier Kinetics), that predicts the radiation-induced conductivity (RIC) produced in a dielectric by an electron beam was developed. The model, which assumes instantly-trapped holes, was then applied to experimental measurements on thin Kapton samples penetrated by an electron beam. Measurements at high bias were matched in the model by an appropriate choice for the trap-modulated electron mobility. A fractional split between front and rear currents measured at zone bias is explained on the basis of beam-scattering. The effects of carrier-enhanced conductivity (CEC) on data obtained for thick, free-surface Kapton samples is described by using an analytical model that incorporates field injection of carriers from the RIC region. The computer code, LWPCHARGE, modified for carrier transport, is also used to predict partial penetration effects associated with CEC in the unirradiated region. Experimental currents and surface voltages, when incorporated in the appropriate models, provide a value for the trap modulated mobility that is in essential agreement with the RIC results.

Meulenberg, A., Jr.↗

An improved semi-global intrinsic kinetics model for high temperature carbon oxidation

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

33 ADVANCED PROPULSION SYSTEMS↗

Derivation of scaled surface reflectances from AVIRIS data

A method for retrieving 'scaled surface reflectances' assuming horizontal surfaces having Lambertian reflectances from spectral data collected by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is presented here. In this method, the integrated water vapor amount on a pixel by pixel basis is derived from the 0.94 micron and 1.14 micron water vapor absorption features. The transmission spectra of H2O, CO2, O3, N2O, CO, CH4, and O2 in the 0.4-2.5 micron region are simulated. The scattering effect due to atmospheric molecules and aerosols is modeled with the 5S computer code. The AVIRIS radiances are divided by solar irradiances above the atmosphere to obtain the apparent reflectances. The scaled surface reflectances are derived from the apparent reflectances using the simulated atmospheric gaseous transmittances and the simulated molecular and aerosol scattering data. The scaled surface reflectances differ from the real surface reflectances by a multiplicative factor. In order to convert the scaled surface reflectances into real surface reflectances, the slopes and aspects of the surfaces must be known.

Gao, Bo-Cai↗

Transonic flow about a thick circular-arc airfoil

An experimental and theoretical study of transonic flow over a thick airfoil, prompted by a need for adequately documented experiments that could provide rigorous verification of viscous flow simulation computer codes, is reported. Special attention is given to the shock-induced separation phenomenon in the turbulent regime. Measurements presented include surface pressures, streamline and flow separation patterns, and shadowgraphs. For a limited range of free-stream Mach numbers the airfoil flow field is found to be unsteady. Dynamic pressure measurements and high-speed shadowgraph movies were taken to investigate this phenomenon. Comparisons of experimentally determined and numerically simulated steady flows using a new viscous-turbulent code are also included. The comparisons show the importance of including an accurate turbulence model. When the shock-boundary layer interaction is weak the turbulence model employed appears adequate, but when the interaction is strong, and extensive regions of separation are present, the model is inadequate and needs further development.

Mcdevitt, J. B.↗

Flow over an all-body hypersonic aircraft - Experiment and computation

The objective of the present investigation is to establish a benchmark experimental data base for a generic hypersonic vehicle shape for validation and/or calibration of advanced computational fluid dynamics computer codes. This paper includes results from the comprehensive test program conducted in the NASA Ames 3.5-ft Hypersonic Wind Tunnel for a generic all-body hypersonic aircraft model. Experimental and computational results on flow visualization, surface pressures, surface convective heat transfer, and pilot-pressure flowfield surveys are presented. Comparisons of the experimental results with computational results from an upwind parabolized Navier-Stokes code developed at NASA Ames demonstrate the capabilities of this code.

William K. Lockman↗

Cosmic topology. Part I. Limits on orientable Euclidean manifolds from circle searches

The Einstein field equations of general relativity constrain the local curvature at every point in spacetime, but say nothing about the global topology of the Universe. Cosmic microwave background anisotropies have proven to be the most powerful probe of non-trivial topology since, within ΛCDM, these anisotropies have well-characterized statistical properties, the signal is principally from a thin spherical shell centered on the observer (the last scattering surface), and space-based observations nearly cover the full sky. The most generic signature of cosmic topology in the microwave background is pairs of circles with matching temperature and polarization patterns. No such circle pairs have been seen above noise in the WMAP or Planck temperature data, implying that the shortest non-contractible loop around the Universe through our location is longer than 98.5% of the comoving diameter of the last scattering surface. We translate this generic constraint into limits on the parameters that characterize manifolds with each of the nine possible non-trivial orientable Euclidean topologies, and provide a code which computes these constraints. In all but the simplest cases, the shortest non-contractible loop in the space can avoid us, and be shorter than the diameter of the last scattering surface by a factor ranging from 2 to at least 6. Further, this result implies that a broader range of manifolds is observationally allowed than widely appreciated. Probing these manifolds will require more subtle statistical signatures than matched circles, such as off-diagonal correlations of harmonic coefficients.

79 ASTRONOMY AND ASTROPHYSICS↗

Vector Radiative Transfer Code SORD: Performance Analysis and Quick Start Guide

We present a new open source polarized radiative transfer code SORD written in Fortran 9095. SORD numerically simulates propagation of monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface using the method of successive orders of scattering (hence the name). Thermal emission is ignored. We did not improve the method in any way, but report the accuracy and runtime in 52 benchmark scenarios. This paper also serves as a quick start users guide for the code available from ftp:maiac.gsfc.nasa.govpubskorkin, from the JQSRT website, or from the corresponding (first) author.

polarized radiative transfer↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Attachment-Line Transition Due to Roughness on a 76 Degree Swept Cylinder at Mach 1.6

Experiments were conducted on a 76 deg. swept cylinder to establish the conditions for the attachment-line transition process with, and without, surface roughness in a low-disturbance ("quiet"), Mach number 1.6 flow. Local flow parameters were estimated from pressure measurements. These were in good agreement with predictions from an Euler code (CFL3D) and a boundary layer code (BL3D). Hot-wires and Schlieren photography were used to determine the state of the boundary layer. It was found that, for a near-adiabatic wall condition and a smooth surface, the attachment-line, boundary-layer remained laminar up to the highest attainable Reynolds number (R-bar of 790). Transition under the influence of trip wires was found to depend on wind-tunnel disturbance levels and the onset conditions have been established. Results suggest that current design practice, which is based upon data from conventional ("noisy") tunnels, may be highly conservative.

Coleman, Colin P.↗

Convection from Hemispherical and Conical Model Ice Roughness Elements in Stagnation Region Flows

To improve ice accretion prediction codes, more data regarding ice roughness and its effects on convective heat transfer are required. The Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research was used to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. In the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with multiple surfaces or sets of roughness panels, each with a different representation of ice roughness. The sets of roughness panels were constructed using two element distribution patterns that were created based on a laser scan of an iced airfoil acquired in the Icing Research Tunnel at NASA Glenn. For both roughness patterns, surfaces were constructed using plastic hemispherical elements, plastic conical elements, and aluminum conical elements. Infrared surface thermometry data from tests run in the VIST were used to calculate area averaged heat transfer coefficient values. The values from the roughness surfaces were compared to the smooth control surface, showing convective enhancement as high as 400% in some cases. The data gathered during this study will ultimately be used to improve the physical modeling in LEWICE or other ice accretion codes and produce predictions of in-flight ice accretion on aircraft surfaces with greater confidence.

Icing Physics↗

Surface passivation of InP solar cells with InAlAs layers

The efficiency of indium phosphide solar cells is limited by high values of surface recombination. The effect of a lattice-matched In(0.52)Al(0.48)As window layer material for InP solar cells, using the numerical code PC-1D is investigated. It was found that the use of InAlAs layer significantly enhances the p(+)n cell efficiency, while no appreciable improvement is seen for n(+)p cells. The conduction band energy discontinuity at the heterojunction helps in improving the surface recombination. An optimally designed InP cell efficiency improves from 15.4 percent to 23 percent AMO for a 10 nm thick InAlAs layer. The efficiency improvement reduces with increase in InAlAs layer thickness, due to light absorption in the window layer.

Jain, Raj K.↗

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Low-gravity electrodeposition and growth of polymer thin films with large third-order optical nonlinearities by electrochemical processes for devices: Thiophene-based polymers

It has been proposed that NLO thin film properties may be improved by low-gravity processing. Strong candidates for NLO thin film applications are the polythiophenes. Polymeric thiophenes are attractive materials due to their ease of preparation, stability, and high X(exp 3). A simple and convenient method for preparation of polythiophenes is electrochemical oxidation. We will apply some of our experience and lessons learned in low-gravity metal, metal/cermet electrode position to improve the quality of polythiophene(s) thin films. In low gravity electrode position of Ni at a high rate on an Au substrate often results in the production of an x-ray non diffracting surface. Cobalt metal deposition does not give this result nor does Ni when deposited similarly on a glassy carbon substrate. Co/Ni alloy composition produced during electrode position is strongly dependent upon the amount of convection. Code position of neutral inert cermets with metals is influenced significantly by the presence of gravity and the size of the cermets. Tracks left in the 1-g surfaces by unsuccessful particle occlusion indicate suspension of the large particles is not the only reason for poor volume percentages of the larger particles in the deposits. All size particles are more homogeneously distributed in the deposits in low-gravity electrocodeposition than in 1-g. Low gravity gives larger volume percentages for the larger particles in the deposits, while 1-g gives larger volume percentages for the smaller particles. Intermediate size particles give mixed results. The experimental cells were constructed with flat electrode end plates such that 1-g bench reference electrode positions could be carried out at various orientations with respect to gravity. A series of bench studies using similar designed cells are suggested so that convection modification can be applied to electrochemical thin film preparation. Convection effects can then be coupled with other parameter variations in current, temperature, concentration, solvent, electrolyte, and anode substrates to optimize the properties of NLO thin films before resorting to low-gravity processing.

Penn, Benjamin G.↗