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 415 records · Page 23

LEWICE/E - An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

A Summary of Validation Results for LEWICE 2.0

A research project is underway at NASA Lewis to produce a computer code which can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from version 2.0 of this code, which is called LEWICE. This version differs from previous releases due to its robustness and its ability to reproduce results accurately for different point spacing, and time step criteria across general computing platforms. It also differs in the extensive amount of effort undertaken to compare the results in a quantifiable manner against the database of ice shapes which have been generated in the NASA Lewis Icing, Research Tunnel (IRT), The complete set of data used for this comparison is available in a recent contractor report . The result of this comparison shows that the difference between the predicted ice shape from LEWICE 2.0 and the average of the experimental data is 7.2% while the variability of the experimental data is 2.5%.

Wright, William B.↗

Validation Results for LEWICE 2.0

A research project is underway at NASA Lewis to produce a computer code which can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from version 2.0 of this code, which is called LEWICE. This version differs from previous releases due to its robustness and its ability to reproduce results accurately for different spacing and time step criteria across computing platform. It also differs in the extensive amount of effort undertaken to compare the results in a quantified manner against the database of ice shapes which have been generated in the NASA Lewis Icing Research Tunnel (IRT). The results of the shape comparisons are analyzed to determine the range of meteorological conditions under which LEWICE 2.0 is within the experimental repeatability. This comparison shows that the average variation of LEWICE 2.0 from the experimental data is 7.2% while the overall variability of the experimental data is 2.5%.

Wright, William B.↗

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deployable reflector antenna performance optimization using automated surface correction and array-feed compensation

Methods for increasing the electromagnetic (EM) performance of reflectors with rough surfaces were tested and evaluated. First, one quadrant of the 15-meter hoop-column antenna was retrofitted with computer-driven and controlled motors to allow automated adjustment of the reflector surface. The surface errors, measured with metric photogrammetry, were used in a previously verified computer code to calculate control motor adjustments. With this system, a rough antenna surface (rms of approximately 0.180 inch) was corrected in two iterations to approximately the structural surface smoothness limit of 0.060 inch rms. The antenna pattern and gain improved significantly as a result of these surface adjustments. The EM performance was evaluated with a computer program for distorted reflector antennas which had been previously verified with experimental data. Next, the effects of the surface distortions were compensated for in computer simulations by superimposing excitation from an array feed to maximize antenna performance relative to an undistorted reflector. Results showed that a 61-element array could produce EM performance improvements equal to surface adjustments. When both mechanical surface adjustment and feed compensation techniques were applied, the equivalent operating frequency increased from approximately 6 to 18 GHz.

Schroeder, Lyle C.↗

Code Description for "Brief Communication: Monitoring snow depth using small, cheap, and easy-to-deploy ground surface temperature sensors"

Temporally continuous snow depth estimates are vital for understanding changing snow patterns and impacts on permafrost in the Arctic. We train a random forest machine learning model to predict snow depth from variability in ground surface temperature. To our knowledge, this is the first time that small ground surface temperature sensors have been used to estimate snow depth. The model performs well at sites where the model was trained and at pan-arctic evaluation sites (RMSE <= 0.15 m). Small temperature sensors are cheap and easy-to-deploy, so this technique enables spatially distributed and temporally continuous snowpack monitoring to an extent previously infeasible. The model is flexible and can be applied to datasets retroactively to retrieve snow depth estimates at additional sites. This code package includes a *.joblib file of the trained random forest model and a *.ipynb file showing how to clean input data, train the random forest model, and apply the model.

Bachand, Claire↗

Cockpit weather information needs

The primary objective is to develop an advanced pilot weather interface for the flight deck and to measure its utilization and effectiveness in pilot reroute decision processes, weather situation awareness, and weather monitoring. Identical graphical weather displays for the dispatcher, air traffic control (ATC), and pilot crew should also enhance the dialogue capabilities for reroute decisions. By utilizing a broadcast data link for surface observations, forecasts, radar summaries, lightning strikes, and weather alerts, onboard weather computing facilities construct graphical displays, historical weather displays, color textual displays, and other tools to assist the pilot crew. Since the weather data is continually being received and stored by the airborne system, the pilot crew has instantaneous access to the latest information. This information is color coded to distinguish degrees of category for surface observations, ceiling and visibilities, and ground radar summaries. Automatic weather monitoring and pilot crew alerting is accomplished by the airborne computing facilities. When a new weather information is received, the displays are instantaneously changed to reflect the new information. Also, when a new surface or special observation for the intended destination is received, the pilot crew is informed so that information can be studied at the pilot's discretion. The pilot crew is also immediately alerted when a severe weather notice, AIRMET or SIGMET, is received. The cockpit weather display shares a multicolor eight inch cathode ray tube and overlaid touch panel with a pilot crew data link interface. Touch sensitive buttons and areas are used for pilot selection of graphical and data link displays. Time critical ATC messages are presented in a small window that overlays other displays so that immediate pilot alerting and action can be taken. Predeparture and reroute clearances are displayed on the graphical weather system so pilot review of weather along the route can be accomplished prior to pilot acceptance of the clearance. An ongoing multiphase test series is planned for testing and modifying the graphical weather system. Preliminary data shows that the nine test subjects considered the graphical presentation to be much better than their current weather information source for situation awareness, flight safety, and reroute decision making.

Scanlon, Charles H.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

Gear optimization

The use of formal numerical optimization methods for the design of gears is investigated. To achieve this, computer codes were developed for the analysis of spur gears and spiral bevel gears. These codes calculate the life, dynamic load, bending strength, surface durability, gear weight and size, and various geometric parameters. It is necessary to calculate all such important responses because they all represent competing requirements in the design process. The codes developed here were written in subroutine form and coupled to the COPES/ADS general purpose optimization program. This code allows the user to define the optimization problem at the time of program execution. Typical design variables include face width, number of teeth and diametral pitch. The user is free to choose any calculated response as the design objective to minimize or maximize and may impose lower and upper bounds on any calculated responses. Typical examples include life maximization with limits on dynamic load, stress, weight, etc. or minimization of weight subject to limits on life, dynamic load, etc. The research codes were written in modular form for easy expansion and so that they could be combined to create a multiple reduction optimization capability in future.

Vanderplaats, G. N.↗

Design of a natural laminar flow wing for a transonic corporate transport

Two- and three-dimensional computational methods were used to design a wing for a transonic corporate transport that had significant runs of laminar flow on both upper and lower wing surfaces at the cruise condition. The airfoil was derived from a low-speed laminar-flow section. The contour was systematically modified based on results from a two-dimensional transonic code to give favorable pressure gradients to 50-percent chord on the upper surface and 65-percent chord on the lower surface. Three-dimensional transonic codes were used to determine the wing twist and to evaluate the aerodynamic characteristics of the complete configuration at various flight conditions. A leading-edge modification for improved stall characteristics was also designed. The theoretical pressure distributions for the final airfoil correlated well with results from wind-tunnel tests.

Campbell, R. L.↗

Surface Modeling and Grid Generation of Orbital Sciences X34 Vehicle

The surface modeling and grid generation requirements, motivations, and methods used to develop Computational Fluid Dynamic volume grids for the X34-Phase 1 are presented. The requirements set forth by the Aerothermodynamics Branch at the NASA Langley Research Center serve as the basis for the final techniques used in the construction of all volume grids, including grids for parametric studies of the X34. The Integrated Computer Engineering and Manufacturing code for Computational Fluid Dynamics (ICEM/CFD), the Grid Generation code (GRIDGEN), the Three-Dimensional Multi-block Advanced Grid Generation System (3DMAGGS) code, and Volume Grid Manipulator (VGM) code are used to enable the necessary surface modeling, surface grid generation, volume grid generation, and grid alterations, respectively. All volume grids generated for the X34, as outlined in this paper, were used for CFD simulations within the Aerothermodynamics Branch.

Alter, Stephen J.↗

Application of engineering codes to the assessment of spacecraft charging induced hazards

One of the major hurdles to developing a better understanding of the environmentally induced interactions in spacecraft systems (spacecraft charging) is the lack of simple tools to analyze the behavior of the spacecraft in a sub-storm environment. The available tool, NASA Charging Analyzer Program (NASCAP), is expensive to use and treats only half of the problem - surface effects only. In response to this need, we have developed, under our IRAD Program, a set of engineering tools that operate on personal computers and addresses both the surface and bulk charging aspects of the interactions. These engineering codes utilize circuit analysis techniques and model simplifications to allow rapid computations. The codes were compared to both NASCAP predictions and to bulk charging data with excellent agreement for both phenomena. The models used in our engineering codes are shown, and in the engineering codes, the models are defined in terms of plates or triangles in an unconfined space. The predictions are based upon current balances to each plate. The code has routines to allow interactions between the plates to simulate 'barrier' effects. The agreement for the predicted differential surface voltages for the model is given. In conducting the comparison of engineering codes to NASCAP predictions, it was found that some inconsistencies existed in the NASCAP code that are not obvious to the user. The use of the grid imposes constraints that voltage gradients exist between the grid points. This allows the existence of voltage gradients in conductors and influences surface resistivity predictions. These points are discussed along with the predictions for treating these inconsistencies properly. The model predictions for a variable intensity substorm is also presented.

Stevens, N. J.↗

Radar transponder antenna pattern analysis for the space shuttle

In order to improve tracking capability, radar transponder antennas will soon be mounted on the Shuttle solid rocket boosters (SRB). These four antennas, each being identical cavity-backed helices operating at 5.765 GHz, will be mounted near the top of the SRB's, adjacent to the intertank portion of the external tank. The purpose is to calculate the roll-plane pattern (the plane perpendicular to the SRB axes and containing the antennas) in the presence of this complex electromagnetic environment. The large electrical size of this problem mandates an optical (asymptotic) approach. Development of a specific code for this application is beyond the scope of a summer fellowship; thus a general purpose code, the Numerical Electromagnetics Code - Basic Scattering Code, was chosen as the computational tool. This code is based on the modern Geometrical Theory of Diffraction, and allows computation of scattering of bodies composed of canonical problems such as plates and elliptic cylinders. Apertures mounted on a curved surface (the SRB) cannot be accomplished by the code, so an antenna model consisting of wires excited by a method of moments current input was devised that approximated the actual performance of the antennas. The improvised antenna model matched well with measurements taken at the MSFC range. The SRB's, the external tank, and the shuttle nose were modeled as circular cylinders, and the code was able to produce what is thought to be a reasonable roll-plane pattern.

Radcliff, Roger↗

Mechanical Erosion Modeling of TPS Materials

The goal of this work is to predict the mechanical response of TPS materials, and specifically, to determine if there is additional surface recession in the heat shield’s surface as a result of mechanical erosion due to the mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within the PATO material response code, enabling it to model the potential mechanical erosion in three steps: first, having the effective mechanical properties as function of temperature, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed with the DPLR hypersonic CFD code, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the TPS material due to surface chemistry and mechanical erosion.

Stress Analysis↗

Organ Dose Assessment and Evaluation of Cancer Risk on Mars Surface

Organ specific fluence spectra and doses for large solar particle events (SPE) and galactic cosmic rays (GCR) at various levels of solar activity are simulated on the surface of Mars using the HZETRN/QMSFRG computer code and the 2010 version of the Badhwar and O Neill GCR model. The NASA JSC propensity model of SPE fluence and occurrence is used to consider upper bounds on SPE fluence for increasing mission lengths. To account for the radiation transmission through the Mars atmosphere, a vertical distribution of Mars atmospheric thickness is calculated from the temperature and pressure data of Mars Global Surveyor. To describe the spherically distributed atmospheric distance on the Mars surface at each elevation, the directional cosine distribution is implemented. The resultant directional shielding by Mars atmosphere at each elevation is then coupled with vehicle and body shielding for organ dose estimates. Finally, cancer risks for astronauts exploring Mars can be assessed by applying the NASA Space Radiation Cancer Risk 2010 model with the resultant organ dose estimates. Variations of organ doses and cancer risk quantities on the surface of Mars, which are due to a 16-km elevation range between the Tharsis Montes and the Hellas impact basin, are visualized on the global topography of Mars measured by the Mars Orbiter Laser Altimeter. It is found that cancer incidence risks are about 2-fold higher than mortality risks with a disproportionate increase in skin and thyroid cancers for male and female astronauts and in breast cancer for female astronauts. The number of safe days, defined by the upper 95% percent confidence level to be below cancer limits, on Mars is analyzed for several Mars mission design scenarios.

Kim, Myung-Hee↗

User's manual for the one-dimensional hypersonic experimental aero-thermodynamic (1DHEAT) data reduction code

A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite-volume heat transfer model are implemented in this code. The analytical solution is based on a one-dimensional, semi-infinite wall thickness model with the approximation of constant substrate thermal properties, which is empirically corrected for the effects of variable thermal properties. The finite-volume solution is based on a one-dimensional, implicit discretization. The finite-volume model directly incorporates the effects of variable substrate thermal properties and does not require the semi-finite wall thickness approximation used in the analytical model. This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either method. This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in this report.

Hollis, Brian R.↗