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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Unsteady blade pressure measurements for the SR-7A propeller at cruise conditions

The unsteady blade surface pressures were measured on the SR-7A propeller. The freestream Mach no., inflow angle, and advance ratio were varied while measurements were made at nine blade stations. At a freestream Mach no. of 0.8, the data in terms of unsteady pressure coefficient vs. azimuth angle are compared to an unsteady 3-D Euler solution, yielding very encouraging results. The code predicts the shape (phase) of the waveform very well, while the magnitude is over-predicted in many cases. At tunnel Mach nos. below 0.6, an unusually large response on the suction surface at 0.15 chord and 0.88 radius was observed. The behavior of this response suggests the presence of a leading edge vortex. The midchord measuring stations on the suction surface exhibit a response that leads the forcing function while most other locations show a phase lag.

Heidelberg, L. J.↗
NASA NTRS1991

Deterministic multi-zone ice accretion modeling

The focus here is on a deterministic model of the surface roughness transition behavior of glaze ice. The initial smooth/rough transition location, bead formation, and the propagation of the transition location are analyzed. Based on the hypothesis that the smooth/rough transition location coincides with the laminar/turbulent boundary layer transition location, a multizone model is implemented in the LEWICE code. In order to verify the effectiveness of the model, ice accretion predictions for simple cylinders calculated by the multizone LEWICE are compared to experimental ice shapes. The glaze ice shapes are found to be sensitive to the laminar surface roughness and bead thickness parameters controlling the transition location, while the ice shapes are found to be insensitive to the turbulent surface roughness.

Yamaguchi, K.↗
NASA NTRS1998

A Fast Code for Jupiter Atmospheric Entry

A fast code was developed to calculate the forebody heating environment and heat shielding that is required for Jupiter atmospheric entry probes. A carbon phenolic heat shield material was assumed and, since computational efficiency was a major goal, analytic expressions were used, primarily, to calculate the heating, ablation and the required insulation. The code was verified by comparison with flight measurements from the Galileo probe's entry; the calculation required 3.5 sec of CPU time on a work station. The computed surface recessions from ablation were compared with the flight values at six body stations. The average, absolute, predicted difference in the recession was 12.5% too high. The forebody's mass loss was overpredicted by 5.5% and the heat shield mass was calculated to be 15% less than the probe's actual heat shield. However, the calculated heat shield mass did not include contingencies for the various uncertainties that must be considered in the design of probes. Therefore, the agreement with the Galileo probe's values was considered satisfactory, especially in view of the code's fast running time and the methods' approximations.

Tauber, Michael E.↗
NASA NTRS1990

Unsteady blade pressure measurements for the SR-7A propeller at cruise conditions

The unsteady blade surface pressures were measured on the SR-7A propeller. The freestream Mach number, inflow angle, and advance ratio were varied while measurements were made at nine blade stations. At a freestream Mach number of 0.8, the data in terms of unsteady pressure coefficient versus azimuth angle are compared to an unsteady three-dimensional Euler solution, yielding very encouraging results. The code predicts the shape (phase) of the waveform very well, while the magnitude is over-predicted in many cases. At tunnel Mach numbers below 0.6, an unusually large response on the suction surface at 0.15 chord and 0.88 radius was observed. The behavior of this response suggests the presence of a leading-edge vortex.

Heidelberg, L. J.↗
DOE OSTI2022

Community Data Mining Approach for Surface Complexation Database Development

This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. Here the framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Pressure-distribution measurements on a transonic low-aspect ratio wing

Experimental surface pressure distributions and oil flow photographs are presented for a 0.90 m semispan model of NASA/Lockheed Wing C, a generic transonic, supercritical, low aspect ratio, highly 3-dimensional configuration. This wing was tested at the design angle of attack of 5 deg over a Mach number range from 0.25 to 0.96, and a Reynolds number range from 3.4 x 1,000,000 to 10 x 1,000,000. Pressures were measured with both the tunnel floor and ceiling suction slots open for most of the tests but taped closed for some tests to simulate solid walls. A comparison is made with the measured pressures from a small model in high Reynolds number facility and with predicted pressures using two three dimesional, transonic full potential flow wing codes: design code FLO22 (nonconservative) and TWING code (conservative). At the given design condition, a small region of flow separation occurred. At a Mach number of 0.82 the flow was unseparated and the surface flow angles were less than 10 deg, indicating that the boundary layer flow was not 3-D. Evidence indicate that wings that are optimized for mild shock waves and mild pressure recovery gradients generally have small 3-D boundary layer flow at design conditions for unseparated flow.

Keener, E. R.↗
NASA NTRS1989

The measurement of boundary layers on a compressor blade in cascade. Volume 1: Experimental technique, analysis and results

Measurements were made of the boundary layers and wakes about a highly loaded, double-circular-arc compressor blade in cascade. These laser Doppler velocimetry measurements have yielded a very detailed and precise data base with which to test the application of viscous computational codes to turbomachinery. In order to test the computational codes at off-design conditions, the data were acquired at a chord Reynolds number of 500,000 and at three incidence angles. Moreover, these measurements have supplied some physical insight into these very complex flows. Although some natural transition is evident, laminar boundary layers usually detach and subsequently reattach as either fully or intermittently turbulent boundary layers. These transitional separation bubbles play an important role in the development of most of the boundary layers and wakes measured in this cascade and the modeling or computing of these bubbles should prove to be the key aspect in computing the entire cascade flow field. In addition, the nonequilibrium turbulent boundary layers on these highly loaded blades always have some region of separation near the trailing edge of the suction surface. These separated flows, as well as the subsequent near wakes, show no similarity and should prove to be a challenging test for the viscous computational codes.

Zierke, William C.↗
NASA NTRS1997

Enhanced Control Effector Designs for Airbreathing Transatmospheric Vehicles

A study was conducted to evaluate the potential effectiveness of a moveable cowl-trailing-edge design for airbreathing hypersonic single-stage-to-orbit (SSTO) configurations, which can be extended or deflected from the nominal position in order to provide additional pitching moment capability. This additional pitching moment capability may reduce the necessary deflection angle from conventional control surfaces and the associated trim drag penalty. Calculations for a generic SSTO configuration with baseline and modified cowl trailing edge geometries were performed at Mach 6 and 10 initially with an engineering analysis code in order to examine several design parametrics. In order to more accurately model geometries and flow physics, 2-D viscous computational fluid dynamics (CFD) predictions were obtained. FInally, a limited set of 3D CFD predictions were obtained at Mach 6 in order to show the effects of modeling 3D flow fields as well as the full 3D vehicle geometry. Comparisons of aftbody surface pressures and force and moment predictions show differences between initial predictions and 2-D CFD solutions due to geometry modeling and calculation method differences. The 3-D CFD predictions confirm the trends observed in the 2-D solutions and provide additional information on 3D effects. These analyses shows that cowl-trailing-edge extensions were effective in providing additional (nose-down) pitching moment increments as well as increased thrust compared to the baseline geometry. These effects reduce the control surface deflection angle required to trim the vehicle and the associated trim drag. Cowl-trailingedge deflections were not as effective, generating more nose-up pitching moment and decreased thrust compared to non-deflected cases.

Cockrell, Charles E.↗
DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
NASA NTRS1985

User's manual for three dimensional boundary layer (BL3-D) code

An assessment has been made of the applicability of a 3-D boundary layer analysis to the calculation of heat transfer, total pressure losses, and streamline flow patterns on the surface of both stationary and rotating turbine passages. In support of this effort, an analysis has been developed to calculate a general nonorthogonal surface coordinate system for arbitrary 3-D surfaces and also to calculate the boundary layer edge conditions for compressible flow using the surface Euler equations and experimental data to calibrate the method, calculations are presented for the pressure endwall, and suction surfaces of a stationary cascade and for the pressure surface of a rotating turbine blade. The results strongly indicate that the 3-D boundary layer analysis can give good predictions of the flow field, loss, and heat transfer on the pressure, suction, and endwall surface of a gas turbine passage.

Anderson, O. L.↗
NASA NTRS1990

Study of vortex breakdown of F-106B by Euler code

The 'Three-dimensional Euler Aerodynamic Method' (TEAM) is presently applied to the F-106B at subsonic speed, in order to examine the relationship between off- and on-surface flow features at angles-of-attack sufficiently great for the occurrence of vortex breakdown. Although TEAM's flow separation is triggered by numerical dissipation, the general trend of vortex-breakdown effect on computed lift characteristics is similar to extant wind tunnel results.

Pao, Jenn Louh↗
NASA NTRS1992

MOM3D method of moments code theory manual

MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.

Shaeffer, John F.↗
NASA NTRS1991

Deterministic multi-zone ice accretion modeling

The study focuses on a deterministic model of the surface roughness transition behavior of glaze ice and analyzes the initial smooth/rough transition location, bead formation, and the propagation of the transition location. Based on a hypothesis that the smooth/rough transition location coincides with the laminar/turbulent boundary-layer transition location, a multizone model is implemented in the LEWICE code. In order to verify the effectiveness of the model, ice accretion predictions for simple cylinders calculated by the multizone LEWICE are compared to experimental ice shapes. The glaze ice shapes are found to be sensitive to the laminar surface roughness and bead thickness parameters controlling the transition location, while the ice shapes are found to be insensitive to the turbulent surface roughness.

Yamaguchi, K.↗
NASA NTRS2005

Advanced Small Perturbation Potential Flow Theory for Unsteady Aerodynamic and Aeroelastic Analyses

An advanced small perturbation (ASP) potential flow theory has been developed to improve upon the classical transonic small perturbation (TSP) theories that have been used in various computer codes. These computer codes are typically used for unsteady aerodynamic and aeroelastic analyses in the nonlinear transonic flight regime. The codes exploit the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP theory was developed methodically by first determining the essential elements required to produce full-potential-like solutions with a small perturbation approach on the requisite Cartesian grid. This level of accuracy required a higher-order streamwise mass flux and a mass conserving surface boundary condition. The ASP theory was further developed by determining the essential elements required to produce results that agreed well with Euler solutions. This level of accuracy required mass conserving entropy and vorticity effects, and second-order terms in the trailing wake boundary condition. Finally, an integral boundary layer procedure, applicable to both attached and shock-induced separated flows, was incorporated for viscous effects. The resulting ASP potential flow theory, including entropy, vorticity, and viscous effects, is shown to be mathematically more appropriate and computationally more accurate than the classical TSP theories. The formulaic details of the ASP theory are described fully and the improvements are demonstrated through careful comparisons with accepted alternative results and experimental data. The new theory has been used as the basis for a new computer code called ASP3D (Advanced Small Perturbation - 3D), which also is briefly described with representative results.

Batina, John T.↗
NASA NTRS1987

An aerothermodynamic study of slender conical vehicles

A numerical study was performed to assess the applicability of some current techniques which can be used for aerothermal predictions over slender spherically blunted cones. Predictions using a viscous-shock-layer method and several engineering approaches were compared with experimental results from flight and ground-based tests, with each other, and with other detailed results. Good agreement was obtained in comparisons with laminar and turbulent heating data from the Reentry F flight vehicle and with the wind-tunnel data. In particular, the viscous-shock-layer method was shown to yield excellent comparisons and should be useful in providing detailed flowfield and surface values for slender blunted cones. Additional predictions were obtained with these methods for two 5-deg half-angle cones with different nose radii to illustrate the effects of nose bluntness and angle of attack on drag and heat transfer. These results demonstrate the benefit of nose blunting with respect to heating and drag for laminar and transitional flow at zero-lift conditions and the benefit of heating reduction at angle of attack. Detailed comparison of the engineering code predictions with the viscous-shock-layer results for these additional cases generally showed good agreement except for the laminar prodictions at angle of attack on the forward cone surface. Nonequilibrium calculations at 0-deg angle of attack showed that substantial benefits of low surface catalysis existed only in the nose-dominated region.

Thompson, R. A.↗
NASA NTRS1998

Locations of Halide Ions in Tetragonal Lysozyme Crystals

Anions play an important role in the crystallization of lysozyme, and are known to bind to the crystalline protein. Previous studies employing X-ray crystallography had found one chloride ion binding site in the tetragonal crystal form of the protein and four nitrate ion binding sites in the monoclinic form. Studies using other approaches have reported more chloride ion binding sites, but their locations were not known. Knowing the precise location of these anions is also useful in determining the correct electrostatic fields surrounding the protein. In the first part of this study the anion positions in the tetragonal form were determined from the difference Fourier map obtained from the lysozyme crystals grown in bromide and chloride solutions under identical conditions. The anion locations were then obtained from standard crystallographic methods and five possible anion binding sites were found in this manner. The sole chloride ion binding site found in previous studies was confirmed. The remaining four sites were new ones for tetragonal lysozyme crystals. However, three of these new sites and the previously found one corresponded to the four unique binding sites found for nitrate ions in monoclinic crystals. This suggests that most of the anion binding sites in lysozyme remain unchanged, even when different anions and different crystal forms of lysozyme are employed. It is unlikely that there are many more anions in the tetragonal lysozyme crystal structure. Assuming osmotic equilibrium it can be shown that there are at most three more anions in the crystal channels. Some of the new anion binding sites found in this study were, as expected, in pockets containing basic residues. However, some of them were near neutral, but polar, residues. Thus, the study also showed the importance of uncharged, but polar groups, on the protein surface in determining its electrostatic field. This was important for the second part of this study where the electrostatic field surrounding the protein was accurately determined. This was achieved by solving the linearized version of the Poisson-Boltzmann equation for the protein in solution. The solution was computed employing the commercial code Delphi which uses a finite difference technique. This has recently become available as a module in the general protein visualization code Insight II. Partial charges were assigned to the polar groups of lysozyme for the calculations done here. The calculations showed the complexity of the electrostatic field surrounding the protein. Although most of the region near the protein surface had a positive field strength, the active site cleft was negatively charged and this was projected a considerable distance. This might explain the occurrence of "head-to-side" interactions in the formation of lysozyme aggregates in solution. Pockets of high positive field strength were also found in the vicinity of the anion locations obtained from the crystallographic part of this study, confirming the validity of these calculations. This study clearly shows not only the importance of determining the counterion locations in protein crystals and the electrostatic fields surrounding the protein, but also the advantage of performing them together.

Lim, Kap↗