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 397 records · Page 22

Introduction of the ASP3D Computer Program for Unsteady Aerodynamic and Aeroelastic Analyses

A new computer program has been developed called ASP3D (Advanced Small Perturbation 3D), which solves the small perturbation potential flow equation in an advanced form including mass-consistent surface and trailing wake boundary conditions, and entropy, vorticity, and viscous effects. The purpose of the program is for unsteady aerodynamic and aeroelastic analyses, especially in the nonlinear transonic flight regime. The program exploits 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 ASP3D code is the result of a decade of developmental work on improvements to the small perturbation formulation, performed while the author was employed as a Senior Research Scientist in the Configuration Aerodynamics Branch at the NASA Langley Research Center. The ASP3D code is a significant improvement to the state-of-the-art for transonic aeroelastic analyses over the CAP-TSD code (Computational Aeroelasticity Program Transonic Small Disturbance), which was developed principally by the author in the mid-1980s. The author is in a unique position as the developer of both computer programs to compare, contrast, and ultimately make conclusions regarding the underlying formulations and utility of each code. The paper describes the salient features of the ASP3D code including the rationale for improvements in comparison with CAP-TSD. Numerous results are presented to demonstrate the ASP3D capability. The general conclusion is that the new ASP3D capability is superior to the older CAP-TSD code because of the myriad improvements developed and incorporated.

Batina, John T.↗

Three Years of Global Positioning System Experience on International Space Station

The International Space Station global positioning system (GPS) receiver was activated in April 2002. Since that time, numerous software anomalies surfaced that had to be worked around. Some of the software problems required waivers, such as the time function, while others required extensive operator intervention, such as numerous power cycles. Eventually enough anomalies surfaced that the three pieces of code included in the GPS unit have been re-written and the GPS units upgraded. The technical aspects of the problems are discussed, as well as the underlying causes that led to the delivery of a product that has had so many problems. The technical aspects of the problems included physical phenomena that were not well understood, such as the affect that the ionosphere would have on the GPS measurements. The underlying causes were traced to inappropriate use of legacy software, changing requirements, inadequate software processes, unrealistic schedules, incorrect contract type, and unclear ownership responsibilities..

Gomez, Susan↗

Three Years of Global Positioning System Experience on International Space Station

The International Space Station global positioning systems (GPS) receiver was activated in April 2002. Since that time, numerous software anomalies surfaced that had to be worked around. Some of the software problems required waivers, such as the time function, while others required extensive operator intervention, such as numerous power cycles. Eventually, enough anomalies surfaced that the three pieces of code included in the GPS unit have been re-written and the GPS units were upgraded. The technical aspects of the problems are discussed, as well as the underlying causes that led to the delivery of a product that has had numerous problems. The technical aspects of the problems included physical phenomena that were not well understood, such as the affect that the ionosphere would have on the GPS measurements. The underlying causes were traced to inappropriate use of legacy software, changing requirements, inadequate software processes, unrealistic schedules, incorrect contract type, and unclear ownership responsibilities.

Gomez, Susan↗

Simulation of a Periodic Jet in a Crossflow with a RANS Solver Using an Unstructured Grid

A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method originally developed by Anderson[1, 2], and is currently supported by members of the Fast Adaptive Aerospace Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configurations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. A detailed description of the FUN3D code can be found in the references given above. The code is under continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two methods for computing gradients, several turbulence models, and several solution methodologies; all in varying states of development. Only the most robust and reliable components, based on experiences with steady aerodynamic simulations, were employed in this work. As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured meshes using a vertex-centered approach. The inviscid terms are evaluated by a flux-difference splitting formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss gradient evaluations are used for viscous and turbulence modeling terms. The discrete spatial operator is combined with a backward time operator which is then solved iteratively using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at each physical time step. For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to model the action of the diaphragm. This boundary condition is described below in more detail. The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid spacing in on the surfaces and in the inviscid regions are indirectly controlled by specifying sources. The viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the spatial variation of the first step off the surface, growth rates, and the termination criterion by providing small problem dependent subroutines.

Atkins, H. L.↗

Comparison of DSMC and CFD Solutions of Fire II Including Radiative Heating

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. These flows may also contain significant radiative heating. To prepare for these missions, NASA is developing the capability to simulate rarefied, ionized flows and to then calculate the resulting radiative heating to the vehicle's surface. In this study, the DSMC codes DAC and DS2V are used to obtain charge-neutral ionization solutions. NASA s direct simulation Monte Carlo code DAC is currently being updated to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced Quantum-Kinetic chemistry model, and to include electronic energy levels as an additional internal energy mode. The Fire II flight test is used in this study to assess these new capabilities. The 1634 second data point was chosen for comparisons to be made in order to include comparisons to computational fluid dynamics solutions. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid. It is shown that there can be quite a bit of variability in the vibrational temperature inferred from DSMC solutions and that, from how radiative heating is computed, the electronic temperature is much better suited for radiative calculations. To include the radiative portion of heating, the flow-field solutions are post-processed by the non-equilibrium radiation code HARA. Acceptable agreement between CFD and DSMC flow field solutions is demonstrated and the progress of the updates to DAC, along with an appropriate radiative heating solution, are discussed. In addition, future plans to generate more high fidelity radiative heat transfer solutions are discussed.

Liechty, Derek S.↗

GEOS-5 Seasonal Forecast System: ENSO Prediction Skill and Bias

The GEOS-5 AOGCM known as S2S-1.0 has been in service from June 2012 through January 2018 (Borovikov et al. 2017). The atmospheric component of S2S-1.0 is Fortuna-2.5, the same that was used for the Modern-Era Retrospective Analysis for Research and Applications (MERRA), but with adjusted parameterization of moist processes and turbulence. The ocean component is the Modular Ocean Model version 4 (MOM4). The sea ice component is the Community Ice CodE, version 4 (CICE). The land surface model is a catchment-based hydrological model coupled to the multi-layer snow model. The AGCM uses a Cartesian grid with a 1 deg × 1.25 deg horizontal resolution and 72 hybrid vertical levels with the upper most level at 0.01 hPa. OGCM nominal resolution of the tripolar grid is 1/2 deg, with a meridional equatorial refinement to 1/4 deg. In the coupled model initialization, selected atmospheric variables are constrained with MERRA. The Goddard Earth Observing System integrated Ocean Data Assimilation System (GEOS-iODAS) is used for both ocean state and sea ice initialization. SST, T and S profiles and sea ice concentration were assimilated.

GEOS-↗

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗

LDEF environment modeling updates

An updated gas dynamics model for gas interactions around the LDEF is presented that includes improved scattering algorithms. The primary improvement is more accurate predictions of surface fluxes in the wake region. The code used is the Integrated Spacecraft Environments Model (ISEM). Additionally, initial results of a detailed ISEM prediction model of the Solar Array Passive LDEF Experiment (SAMPLE), A0171, is presented. This model includes details of the A0171 geometry and outgassing characteristics of the many surfaces on the experiment. The detailed model includes the multiple scattering that exists between the ambient atmosphere, LDEF outgassing, and atomic oxygen erosion products. Predictions are made for gas densities, surface fluxes and deposition at three different time periods of the LDEF mission.

Gordon, Tim↗

Effect of surface catalysis on heating to ceramic coated thermal protection systems for transatmospheric vehicles

This paper describes the effect of surface catalysis on the heat transfer rate to the heat shield of a typical Transatmospheric Vehicle (TAV) during ascent and atmospheric entry. Surface kinetics and optical properties obtained from arc-jet tests on candidate thermal protection systems (coated metals) were used in a reacting boundary layer code to estimate the heating distribution along the surface of a TAV. Thermochemical stability of the coatings is described in terms of reduction in emittance and loss of opacifiers from the coatings during the arc-jet tests.

Stewart, David A.↗

Validation of the Hypersolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

Validation of the HyperSolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

Single Airfoil Gust Response Problem: Category 3, Problem 1

An unsteady aerodynamic code, called GUST3D (ref. 3), has been developed to solve equation (8) for flows with periodic vortical disturbances. The code uses a frequency-domain approach with second-order central differences and a pressure radiation condition in the far field. GUST3D requires as input certain mean flow quantities which are calculated separately by a potential flow solver. The solver calculates the mean ow using a Gothert's Rule approximation (ref. 3). On the airfoil surface, it uses the solution calculated by the potential code FLO36 (ref. 4). Figures 1-2 show the mean pressure along the airfoil surface for the two airfoil geometries. In Figures 3 - 8, we present the RMS pressure on the airfoil surface. Each figure shows three GUST3D solutions (calculated on grids with different far-field boundary locations). Three solutions are shown to provide some indication of the numerical uncertainty in the results. Figures 9 - 13 present the acoustic intensity. We again show three solutions per case. Note that no results are presented for the k1 = k2 = 2.0 loaded airfoil case, as an acceptable solution could not be obtained. A few comments need to be made about the results shown. First, since the last Workshop, the GUST3D code has been substantially upgraded. This includes implementing a more accurate far-field boundary condition (ref. 5) and developing improved gridding capabilities. This is the reason for any differences that may exist between the present results and results from the last Workshop. Second, the intensity results on the circle R = 4C were obtained using a Kirchoff method (ref. 6). The Kirchoff surface was the circle R = 2C. Finally, the GUST3D code is most accurate for low reduced frequencies. A new domain decomposition approach (ref. 7) has been developed to improve accuracy. Both the single domain and domain decomposition approaches were used in generating the present results.

Scott, James R.↗

Los Alamos RAGE Simulations of the HAIV Mission Concept

The mitigation of potentially hazardous objects (PHOs) can be accomplished by a variety of methods including kinetic impactors, gravity tractors and several nuclear explosion options. Depending on the available lead time prior to Earth impact, non- nuclear options can be very effective at altering a PHOs orbit. However if the warning time is short nuclear options are generally deemed most effective at mitigating the hazard. The NIAC mission concept for a nuclear mission has been presented at several meetings, including the last PDC (2013).We use the adaptive mesh hydrocode RAGE to perform detailed simulations of this Hypervelocity Asteroid Intercept Vehicle (HAIV) mission concept. We use the RAGE code to simulate the crater formation by the kinetic impactor as well as the explosion and energy coupling from the follower nuclear explosive device (NED) timed to detonate below the original surface to enhance the energy coupling. The RAGE code has been well validated for a wide variety of applications. A parametric study will be shown of the energy and momentum transfer to the target 100 m diameter object: 1) the HAIV mission as planned; 2) a surface explosion and 3) a subsurface (contained) explosion; both 2) and 3) use the same source energy as 1).Preliminary RAGE simulations show that the kinetic impactor will carve out a surface crater on the object and the subsequent NED explosion at the bottom of the crater transfers energy and momentum to the target effectively moving it off its Earth crossing orbit. Figure 1 shows the initial (simplified) RAGE 2D setup geometry for this study. Figure 2 shows the crater created by the kinetic impactor and Figure 3 shows the time sequence of the energy transfer to the target by the NED.

Asteroid Mitigation↗

Stellarator coil design using cubic splines for improved access on the outboard side

In recent years many efforts have been undertaken to simplify coil designs for stellarators due to the difficulties in fabricating non-planar coils. Additionally, the FOCUS code removes the need for a winding surface and represents the coils as arbitrary curves in 3D. In the following work, the implementation of a spline representation for the coils in FOCUS is described, along with the implementation of a new engineering constraint to design coils with a straighter outer section. The new capabilities of the code are shown as an example on HSX, NCSX, and a prototype quasi-axisymmetric reactor-sized stellarator. The flexibility granted by splines along with the new constraint will allow for stellarator coil designs with improved accessibility and simplified maintenance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Guide to AERO2S and WINGDES Computer Codes for Prediction and Minimization of Drag Due to Lift

The computer codes, AER02S and WINGDES, are now widely used for the analysis and design of airplane lifting surfaces under conditions that tend to induce flow separation. These codes have undergone continued development to provide additional capabilities since the introduction of the original versions over a decade ago. This code development has been reported in a variety of publications (NASA technical papers, NASA contractor reports, and society journals). Some modifications have not been publicized at all. Users of these codes have suggested the desirability of combining in a single document the descriptions of the code development, an outline of the features of each code, and suggestions for effective code usage. This report is intended to supply that need.

Carlson, Harry W.↗

Radar Reflectivity of the Martian Polar Regions

Radar experiments provide a unique method of probing the surfaces and subsurfaces of planetary bodies. Information on surface and subsurface structures and properties can be extracted from radar data. There is a well developed history of radar investigations of the planet Mars, beginning with the first reports of variations of scattering properties as a function of martian longitude in the mid-1960's. Because of the rapid rotation of Mars, the standard technique of delay-doppler mapping cannot be used as effectively as it is on Mercury and Venus, making it a more difficult radar target. Techniques have been developed to overcome this difficulty, including an inversion technique to combine many doppler-only (CW) experiments from different viewing geometries into a map of surface reflectivity, and new random long-code techniques. These techniques still suffer from ambiguities, however. A third technique which does not suffer from these ambiguities is the combination of the powerful Goldstone transmitter with the VLA as the receiving instrument to create a combined radar imaging instrument. We have used this combined radar instrument to image the surface of Mars in 3.5-cm radar reflectivity during the 1988, 1992/93, and 1999 oppositions. During the 1988 experiments, the residual south polar ice cap (RSPIC) was the brightest radar reflector on the planet - intrinsically brighter than even the Tharsis lava flows. This was quite an unexpected result. In contrast, during the 1992/93 experiments, the residual north polar ice cap (RNPIC) was not nearly so bright, and in fact showed no enhancement at all, This was puzzling, given the 1988 results for the RSPIC. We attributed the lack of a radar reflection enhancement to a combination of three effects: the geometry was different; the season was different; and the intrinsic scattering from the residual ice caps was different. The 1999 experiments provided a chance to test the relative importance of these three effects, since the RNPIC would be visible in a geometry very similar to the RSPIC in 1988, and the season would also be very similar. Preliminary reduction of the data taken in 1999 shows that the RNPIC did in fact show a radar reflectivity enhancement, but that it was still not as bright as the RSPIC in 1988. Details regarding the radar reflectivity of the residual ice caps, as well as the surrounding polar layered terrains (in light of our radar reflectivity data) will be discussed. Additional information is contained in the original extended abstract.

Butler, B. J.↗

GPU-friendly surface model for Monte-Carlo detector simulations

The demands for Monte-Carlo simulation are drastically increasing with the Large Hadron Collider’s high-luminosity upgrade, and are expected to exceed the currently available compute resources. At the same time, modern high-performance computing has adopted powerful hardware accelerators, particularly GPUs. The AdePT and Celeritas projects aim to address the demanding computational needs by leveraging these heterogeneous computing architectures. While both have successfully ported realistic detector simulations to GPUs using the VecGeom library, the complexity of geometry modeling emerged as a bottleneck. Thread divergence and high register usage were degrading the GPU performance. Therefore, a new, GPU-friendly surface-based model has been introduced in the VecGeom library that decomposes the divergent code of the 3D primitive solids into simpler and more balanced surface algorithms. In this work, we present the latest developments, focusing on the additions required to efficiently model complex setups like the CMS Phase-2 geometry. This includes memory reduction techniques, and adding accelerating structures for faster traversal.

Diederichs, Severin [CERN]↗

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.↗