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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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) [1,2,3]. 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) [4] 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) [5] 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 [6]. 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 [7]. 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.

Thermal Protection Systems↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1994

Land surface temperature measurements for EOS MODIS data

Work accomplished includes: Beta delivery 1 of the MODIS LST product; the first version of MODIS LST ATBD; update of the atmospheric radiative transfer code ATRAD; the development of a new approach look-up table method; and improvement of the TIR spectrometer. Preliminary feasibility analysis of the look-up table approach is presented in terms of showing the effects on the TIR radiance at the top of the atmosphere of the stratospheric and upper atmospheric temperature profiles, the surface emissivity and temperature, the lower atmospheric temperature and water vapor profiles, and the viewing angle.

Wan, Zhengming↗
NASA NTRS2011

Thermal Integration of a Liquid Acquisition Device into a Cryogenic Feed System

Primary objectives of this effort were to define the following: (1) Approaches for quantification of the accumulation of thermal energy within a capillary screen liquid acquisition device (LAD) for a lunar lander upper stage during periods of up to 210 days on the lunar surface, (2) techniques for mitigating heat entrapment, and (3) perform initial testing, data evaluation. The technical effort was divided into the following categories: (1) Detailed thermal modeling of the LAD/feed system interactions using both COMSOL computational fluid device and standard codes, (2) FLOW-3D modeling of bulk liquid to provide interfacing conditions for the LAD thermal modeling, (3) condensation conditioning of capillary screens to stabilize surface tension retention capability, and (4) subscale testing of an integrated LAD/feed system. Substantial progress was achieved in the following technical areas: (1) Thermal modeling and experimental approaches for evaluating integrated cryogen LAD/feed systems, at both the system and component levels, (2) reduced gravity pressure control analyses, (3) analytical modeling and testing for capillary screen conditioning using condensation and wicking, and (4) development of rapid turnaround testing techniques for evaluating LAD/feed system thermal and fluid integration. A comprehensive effort, participants included a diverse cross section of representatives from academia, contractors, and multiple Marshall Space Flight Center organizations.

Hastings, L. J.↗
NASA NTRS2021

Analysis of simultaneous aerosol and ocean glint retrieval using multi-angle observations

Since early 2000, NASA's Multi-angle Imaging SpectroRadiometer (MISR) instrument has been performing remote sensing retrievals of aerosol optical properties from the polar-orbiting Terra spacecraft. A noteworthy aspect of MISR observations over the ocean is that, for much of the Earth, some of the multi-angle views have contributions from solar reflection by the ocean surface (glint, or glitter), while others do not. Aerosol retrieval algorithms often discard these glint-influenced observations because they can overwhelm the signal and are difficult to predict without knowledge of the (wind-speed-driven) ocean surface roughness. However, theoretical studies have shown that multi-angle observations of a location at geometries with and without reflected sun glint can be a rich source of information, sufficient to support simultaneous retrieval of both the aerosol state and the wind speed at the ocean surface. We are in the early stages of creating such an algorithm. In this paper, we describe our assessment of the appropriate level of parameterization for simultaneous aerosol and ocean surface property retrievals using sun glint. For this purpose, we use generalized nonlinear retrieval analysis (GENRA), an information content assessment (ICA) technique employing Bayesian inference, and simulations from the Ahmad–Fraser iterative radiative transfer code. We find that four parameters are suitable: aerosol optical depth (τ), particle size distribution (expressed as the fine mode fraction f of small particles in a bimodal size distribution), surface wind speed (w), and relative humidity (r, to define the aerosol water content and complex refractive index). None of these parameters define ocean optical properties, as we found that the aerosol state could be retrieved with the nine MISR near-infrared views alone, where the ocean body is strongly absorbing in the open ocean. We also found that retrieval capability varies with observation geometry and that as τ increases so does the ability to determine aerosol intensive optical properties (r and f, while it decreases for w). Increases in w decrease the ability to determine the true value of that parameter but have minimal impact on retrieval of aerosol properties. We explored the benefit of excluding the two most extreme MISR view angles for which radiative transfer with the plane-parallel approximation is less certain, but we found no advantage in doing so. Finally, the impact of treating wind speed as a scalar parameter, rather than as a two-parameter directional wind, was tested. While the simpler scalar model does contribute to overall aerosol uncertainty, it is not sufficiently large to justify the addition of another dimension to parameter space. An algorithm designed upon these principles is in development. It will be used to perform an atmospheric correction with MISR for coincident ocean color (OC) observations by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument, also on the NASA Terra spacecraft. Unlike MISR, MODIS is a single-view-angle instrument, but it has a more complete set of spectral channels ideal for determination of optical ocean properties. The atmospheric correction of MODIS OC data can therefore benefit from MISR aerosol retrievals. Furthermore, higher-spatial-resolution data from coincident MISR observations may also improve glint screening.

Kirk Knobelspiesse↗
NASA NTRS1990

Three dimensional full potential method for the aeroelastic modeling of propfans

A three-dimensional time dependent full potential aerodynamic analysis of single rotation propellers has been developed. The primary purpose of the code is to provide a capability of doing propfan aeroelastic analysis in the nonlinear transonic regime. A secondary purpose is to provide a validation of the unsteady lifting surface panel method that has been have developed. Results will be shown for steady state aerodynamic loading, unsteady aerodynamic response to forced aeroelastic deformations, and free aeroelastic response. Comparisons are made to experimental data and corresponding panel code results.

Williams, Marc H.↗
NASA NTRS2014

Quaoar: New, Longitudinaly Resolved, Spectroscopic Characterization of Its Surface

(50000) Quaoar, one of the largest Trans-neptunian objects, is comparable in size to Pluto's moon Charon. However, while Charon's surface is rich almost exclusively in H2O ice, Quaoar's surface characterized by ices of CH4, N2, as well as C2H6, a product of irradiation of CH4 (Dalle Ore et al. 2009). Because of its distance from the Sun, Quaoar is expected to have preserved, to a degree, its original composition, however, its relatively small size did not make it a prime candidate for presence of volatile ices in the study by Schaller and Brown (2007). Furthermore, based on the Brown et al. (2011) study (Brown, Schaller, & Fraser, 2011. A Hypothesis for the Color Diversity of the Kuiper Belt. ApJL, 739, L60) its red coloration points to CH3OH as the ice which, when irradiated, might have produced the red material. We present new visible to near-infrared (0.3-2.48 micrometers) spectro-photometric data obtained with the XSHOOTER (Vernet et al. 2011, A&A, 536A, 105 ) instrument at the VLT-ESO facility at four different longitudes on the surface of Quaoar. The data are complemented by previously published photometric observations obtained in the near-infrared (3.6, 4.5 micrometers) with the Spitzer Space Telescope, which provide an extra set of constraints in the model calculation process in spite of the different observing times that preclude establishing the spatial consistency between the two sets. For each of the four spectra we perform spectral modeling of the entire wavelength range -from 0.3 to 4.5 micrometers- by means of a code based on the Shkuratov radiative transfer formulation of the slab model. We obtain spatially resolved compositional information for the surface of Quaoar supporting the presence of CH4 and C2H6, as previously reported, along with evidence for N2 and NH3OH. The albedo at the two Spitzer bands indicates the likely presence of CO and CO2. CH3OH, predicted on the basis of Quaoar's coloration (Brown et al. 2011), is not found at any of the four longitudes, implying that the presence of this ice is a sufficient, but not necessary condition for reddening of TNO surfaces. Other ices, in particular CH4 (Brunetto et al. 2006), have been shown to be plausible precursors for reddening of TNO surfaces.

Ices↗
DOE OSTI2025

Verification of the Shafranov shift in free-boundary VMEC, DESC, and SPEC calculations with elliptical geometry

Here, given a set of current-carrying filaments that create a magnetic field with rotating elliptical flux surfaces, we present verification calculations using the VMEC, DESC, and SPEC free-boundary magnetohydrostatic equilibrium codes for both the “vacuum” equilibrium, for which there is no plasma pressure and no plasma currents, and for a non-zero pressure, zero toroidal current equilibrium. For the vacuum case, these codes are quantitatively compared to the magnetic field produced by the coils. For the non-zero pressure case, they are compared to each other. As the stepped pressure profile used in SPEC approaches the continuous profile used in VMEC and DESC, the SPEC magnetic axis approaches the VMEC and DESC axes. To within an error given by the difference between the pressure profiles, the VMEC, DESC, and SPEC calculations give the same pressure-induced geometric shift of the magnetic axis, known as the Shafranov shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations with Ly$α$ Forests

Here, we present a measurement of baryonic acoustic oscillations (BAO) from Lyman-α (Lyα) absorption and quasars at an effective redshift $z = 2.33$ using the complete extended Baryonic Oscillation Spectroscopic Survey (eBOSS). The sixteenth and final eBOSS data release (SDSS DR16) contains all data from eBOSS and its predecessor, the Baryonic Oscillation Spectroscopic Survey (BOSS), providing 210,005 quasars with $z_q > 2.10$ that are used to measure Lyα absorption. We measure the BAO scale both in the auto-correlation of Ly$α$ absorption and in its cross correlation with 341,468 quasars with redshift $z_q > 1.77$. Apart from the statistical gain from new quasars and deeper observations, the main improvements over previous work come from more accurate modeling of physical and instrumental correlations and the use of new sets of mock data. Combining the BAO measurement from the auto- and cross-correlation yields the constraints of the two ratios $D_H(z = 2.33)/r_d = 8.99 ± 0.19$ and $D_M(z = 2.33)/r_d = 37.5 ± 1.1$, where the error bars are statistical. These results are within $1.5σ$ of the prediction of the flat-ΛCDM cosmology of Planck (2016). The analysis code, picca, the catalog of the flux-transmission field measurements, and the $Δ\chi^2$ surfaces are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Calculation of metric coefficients for streamline coordinates.

A procedure is given for deriving the equations describing the surface streamline metric in a general form and in a form suitable for incorporation into computer codes to calculate an inviscid flow about three-dimensional bodies. Sample results are included to show the application of the metric in conjunction with an axisymmetric analog to predict heat transfer to a typical space shuttle orbiter.

Rakich, J. V.↗
NASA NTRS1976

Ground location of satellite scanner data

This paper presents simple and accurate mathematical formulation for determining the ground location of remote sensor data. The techniques used are based on elementary concepts of differential geometry and lead to the development of a relation that gives location as a function of surface ellipticity, satellite position, velocity, attitude, and scanner orientation. The formula lends itself to simply computer coding and will hopefully lead to a standardization of the various techniques which have been developed to solve this problem.

Puccinelli, E. F.↗
NASA NTRS1978

The practical impact of elastohydrodynamic lubrication

The use of elastohydrodynamics in the analysis of rolling element bearings is discussed. Relationships for minimum film thickness and tractive force were incorporated into computer codes and used for bearing performance prediction. The lambda parameter (ratio of film thickness to composite surface roughness) was shown to be important in predicting bearing life and failure mode. Results indicate that at values of lambda below 3 failure modes other than the classic subsurface initiated fatigue can occur.

Anderson, W. J.↗
NASA NTRS1980

Destruction of N/2D/ by O2 - A major source of 6300 A dayglow emission

The paper discusses the N(2D) + O2 yields NO + O(1D), for k1 reaction, and examines under different geophysical conditions, the effect of this source of 6300 A surface brightness on the dayglow redline. In calculating the 6300 A volume emission rate, a computer code is used which provides the thermal electron density, photoelectron flux, ion concentrations, electron and ion temperatures and odd nitrogen densities required to evaluate the O(1D) sources and sinks. It is found that under most conditions this source is the dominant mechanism responsible for the production of the 6300 A dayglow above 150 km. The volume emission rate of this source in summer is 2.5 times larger than the sum of all other sources at 180 km, and 1.5 times as large in winter. In addition, when the N(2D) source of O(1D) is taken into account, the rate coefficient for quenching of O(1D) by N2 is estimated to be 4 plus or minus 1 x 10 to the -11th cu cm per sec.

Torr, D. G.↗
NASA NTRS1985

Use of a corrugated surface to enhance radiation tolerance in a GaAs solar cell

The use of a corrugated surface on a GaAs solar cell and its effects on radiation resistance were studied. A compute code was developed to determine the performance of the cell for various geometric parameters. The large optical absorption coefficient of GaAs allows grooves to be only 4-5 micrometers deep. Using accepted material parameters for GaAs solar cells the theoretical performances were compared for various corrugated cells before and after minority carrier diffusion length degradation. The total power output was maximized for both n(+)/p and p(+)/n cells. Optimum values of 1.0-1.5 and 5.0 micrometers for groove and ridge widths respectively were determined.

Leon, Rosa P.↗
NASA NTRS1988

Vibration and flutter characteristics of the SR7L large-scale propfan

An investigation of the vibration characteristics and aeroelastic stability of the SR7L Large-Scale Advanced Propfan was performed using a finite element blade model and an improved aeroelasticity code. Analyses were conducted for different blade pitch angles, blade support conditions, number of blades, rotational speeds, and freestream Mach numbers. A finite element model of the blade was used to determine the blade's vibration behavior and sensitivity to support stiffness. The calculated frequencies and mode shape obtained with this model agreed well with the published experimental data. A computer code recently developed at NASA Lewis Research Center and based on three-dimensional, unsteady, lifting surface aerodynamic theory was used for the aeroelastic analysis to examine the blade's stability at a cruise condition of Mach 0.8 at 1700 rpm. The results showed that the blade is stable for that operating point. However, a flutter condition was predicted if the cruise Mach number was increased to 0.9.

August, Richard↗
NASA NTRS1989

Thermal analysis of engine inlet anti-icing systems

A hot air anti-icing system of a gas turbine engine inlet is analyzed numerically. A three-dimensional potential flow code, which accounts for compressibility effects, is used to determine the flowfield in and around the inlet. A particle trajectory code is developed using a local linearization technique. The trajectory code is used to calculate local water impingement rates. Energy balances are performed on both the surface runback water and the metallic skin to determine their temperature distributions. A variety of test cases are considered in order to validate the various numerical components of the process as well as to demonstrate the procedure.

Keith, Theo G., Jr.↗
NASA NTRS1992

Grid generation for a complex aircraft configuration

The procedure used to create a grid around the F/A-18 fighter aircraft is presented. This work was done for the NASA High Alpha Technology Program. As part of this program, LeRC is numerically and experimentally investigating the flow in the F/A-18 inlet duct at high angles of attack. A grid was needed which could be used to calculate both the external and internal flow around the F/A-18. The grid had to be compatible with the computational fluid dynamics (CFD) codes PARC3D and CFL3D. The programs used to create this grid were I3GVIRGO and GRIDGEN. A surface definition used to create the grid was obtained from McDonnell Aircraft Company (MCAIR) and was composed of numerous files each containing a point definition of a portion of the aircraft. These files were read into the geometry manipulation program I3GVIRGO, where they were modified and grouped into smaller GRIDGEN database files. Next, the block outlines and boundary conditions were specified in the GRIDBLOCK program. The GRIDGEN2D program was used to create the surface grid on the block faces, and GRIDGEN3D was used to create the full 3-D grid.

Bruns, Jim↗