Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code Analysis”

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 883 records · Page 49

A Dynamic Calibration Method for Experimental and Analytical Hub Load Comparison

This paper presents the results from an ongoing effort to produce improved correlation between analytical hub force and moment prediction and those measured during wind-tunnel testing on the Aeroelastic Rotor Experimental System (ARES), a conventional rotor testbed commonly used at the Langley Transonic Dynamics Tunnel (TDT). A frequency-dependent transformation between loads at the rotor hub and outputs of the testbed balance is produced from frequency response functions measured during vibration testing of the system. The resulting transformation is used as a dynamic calibration of the balance to transform hub loads predicted by comprehensive analysis into predicted balance outputs. In addition to detailing the transformation process, this paper also presents a set of wind-tunnel test cases, with comparisons between the measured balance outputs and transformed predictions from the comprehensive analysis code CAMRAD II. The modal response of the testbed is discussed and compared to a detailed finite-element model. Results reveal that the modal response of the testbed exhibits a number of characteristics that make accurate dynamic balance predictions challenging, even with the use of the balance transformation.

Kreshock, Andrew R.↗

Pilot Deployment of the LDSD Parachute via a Supersonic Ballute

The Low Density Supersonic Decelerator (LDSD) Project required the use of a pilot system due to the inability to mortar deploy its main supersonic parachute. A mortar deployed 4.4 m diameter supersonic ram-air ballute was selected as the pilot system for its high drag coefficient and stability relative to candidate supersonic parachutes at the targeted operational Mach number of 3. The ballute underwent a significant development program that included the development of a new liquid methanol-based pre-inflation system to assist the ballute inflation process. Both pneumatic and pyrotechnic mortar tests were conducted to verify orderly rigging deployment, bag strip, inflation aid activation, and proper mortar performance. The ballute was iteratively analyzed between fluid and structural analysis codes to obtain aerodynamic and aerothermodynamic estimates as well as estimates of the ballute's structural integrity and shape. The ballute was successfully flown in June 2014 at a Mach number of 2.73 as part of the first LDSD supersonic flight test and performed beyond expectations. Recovery of the ballute indicated that it did not exceed its structural or thermal capabilities. This flight set a historical precedent as it represented the largest ballute to have ever been successfully flown at this Mach number by a NASA entity.

Tanner, Christopher L.↗

Fluid Transient Analysis during Priming of Evacuated Line

Water hammer analysis in pipe lines, in particularly during priming into evacuated lines is important for the design of spacecraft and other in-space application. In the current study, a finite volume network flow analysis code is used for modeling three different geometrical configurations: the first two being straight pipe, one with atmospheric air and other with evacuated line, and the third case is a representation of a complex flow network system. The numerical results show very good agreement qualitatively and quantitatively with measured data available in the literature. The peak pressure and impact time in case of straight pipe priming in evacuated line shows excellent agreement.

Bandyopadhyay, Alak↗

Use of a Hybrid Edge Node-Centroid Node Approach to Thermal Modeling

A recent proposal submitted for an ESA mission required that models be delivered in ESARAD/ESAT AN formats. ThermalDesktop was the preferable analysis code to be used for model development with a conversion done as the final step before delivery. However, due to some differences between the capabilities of the two codes, a unique approach was developed to take advantage of the edge node capability of ThermalDesktop while maintaining the centroid node approach used by ESARAD. In essence, two separate meshes were used: one for conduction and one for radiation. The conduction calculations were eliminated from the radiation surfaces and the capacitance and radiative calculations were eliminated from the conduction surfaces. The resulting conduction surface nodes were coincident with all nodes of the radiation surface and were subsequently merged, while the nodes along the edges remained free. Merging of nodes on the edges of adjacent surfaces provided the conductive links between surfaces. Lastly, all nodes along edges were placed into the subnetwork and the resulting supernetwork included only the nodes associated with radiation surfaces. This approach had both benefits and disadvantages. The use of centroid, surface based radiation reduces the overall size of the radiation network, which is often the most computationally intensive part of the modeling process. Furthermore, using the conduction surfaces and allowing ThermalDesktop to calculate the conduction network can save significant time by not having to manually generate the couplings. Lastly, the resulting GMM/TMM models can be exported to formats which do not support edge nodes. One drawback, however, is the necessity to maintain two sets of surfaces. This requires additional care on the part of the analyst to ensure communication between the conductive and radiative surfaces in the resulting overall network. However, with more frequent use of this technique, the benefits of this approach can far outweigh the additional effort.

Peabody, Hume L.↗

Integration of Information Management System, Workflow and Computational Tools Enabling Multiscale Modeling Within an ICME Paradigm

With the increased emphasis on reducing the cost and time to market of new materials, the need for analytical tools that enable the virtual design and optimization of materials throughout their processing - internal structure - property - performance envelope, along with the capturing and storing of the associated material and model information across its lifecycle, has become critical. This need is also fueled by the demands for higher efficiency in material testing; consistency, quality and traceability of data; product design; engineering analysis; as well as control of access to proprietary or sensitive information. Fortunately, material information management systems and physics-based multiscale modeling methods have kept pace with the growing user demands. Herein, recent efforts to develop a set of Python functions that exchange information between NASA GRC's Integrated multiscale Micromechanics Analysis Code (ImMAC) software toolset and its Integrated Computational Materials Engineering (ICME), Granta MI® database schema is presented. The goal is to enable seamless coupling between both test data and simulation data, which is captured and tracked automatically within Granta MI®, with full model pedigree information. These tools, and this type of linkage, are foundational to realizing the full potential of ICME, in which materials processing, microstructure, properties, and performance are coupled to enable application-driven design and optimization of materials and structures.

multiscale modeling; Micromechanics; Computational↗

Generation of Mars Helicopter Rotor Model for Comprehensive Analyses

The present research is aimed at providing a performance model for the Mars Helicopter (MH), to understand the complexity of the flow, and identify future regions of improvement. The low density of the Martian atmosphere and the relatively small MH rotor, result in very low chord-based Reynolds number flows. The low density and Reynolds numbers reduce the lifting force and lifting efficiency, respectively. The high drag coefficients in subcritical flow, especially for thicker sections, are attributed to laminar separation from the rear of the airfoil. In the absence of test data, efforts have been made to explore these effects using prior very low Reynolds number research efforts. The rotor chord-based Reynolds number range is observed to be subcritical, which makes boundary layer transition unlikely to occur. The state of the two-dimensional rotor boundary layer in hover is approximated by calculating the instability point, laminar separation point, and the transition location to provide understanding of the flow state in the high Mach-low Reynolds number regime. The results are used to investigate the need for turbulence modeling in Computational Fluid Dynamics (CFD) calculations afterwards. The goal is to generate a performance model for the MH rotor for a free wake analysis, because the computational budget for a complete Navier-Stokes solution for a rotating body-fitted rotor is substantial. In this study, a Reynolds-Averaged Navier-Stokes (RANS) based approach is used to generate the airfoil deck using C81Gen with stitched experimental data for very high angles of attack. A full Grid Resolution Study is performed and over 4,500 cases are completed to create the full airfoil deck. The laminar separation locations are predicted within the accuracy of the approximate method when compared with the CFD calculations. The model is presented through airfoil data tables (c81 files) that are used by comprehensive rotor analysis codes such as CAMRADII, or the mid-fidelity CFD solver RotCFD. Finally, the rotor performance is compared with experimental data from the 25ft Space Simulator at the NASA Jet Propulsion Laboratory (JPL) and shows good correlation for the rotor Figure of Merit over the available thrust range.

Model↗

Low Reynolds Number Airfoil Evaluation for the Mars Helicopter Rotor

The present research is aimed at providing a performance model for the Mars Helicopter (MH), to understand the complexity of the flow, and future regions of improvement. The Martian atmosphere's low density and the MH's relatively small rotor result in very low chord-based Reynolds number flows, Rec = O(10(exp 3)-10(exp 4)). The low density and subcritical Reynolds number reduce the lifting force and lifting efficiency, respectively. The high drag coefficients in subcritical flow, especially for thicker sections, are attributed to laminar separation from the rear of the airfoil. The goal is to generate a performance model for the MH rotor for a free wake analysis, since the computational budget for a complete Navier-Stokes solution for a rotating body-fitted rotor is substantial. In this study, a RANS-based approach is used to generate the airfoil deck using OVERFLOW with stitched experimental data for very high angles of attack. The model is presented through airfoil data tables (C81 files) that are used by comprehensive rotor analysis codes such as CAMRADII, or the mid-fidelity CFD solver RotCFD. These codes have proven to provide accurate performance predictions for all rotor operations at only a fraction of the computational expense of three- dimensional body-fitted viscous grids.

Koning, Witold J. K.↗

Fuselage-Induced Velocity Model for LRTA, RTA, and HART Fuselages

In order to make the aerodynamic fuselage-rotor interference effects available to comprehensive rotor codes, a simple analytical model of the fuselage-induced velocities within the volume of rotor blade operation above the fuselage is developed here for the following bodies used in wind tunnel experiments: the Large Rotor Test Apparatus (LRTA), the Rotor Test Apparatus (RTA), and the Higher Harmonic Control Aeroacoustic Rotor Test (HART).While the first two are used in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames, California, the third one is used by DLR in the large low-speed facility of the German-Dutch wind tunnel in the Netherlands. The fuselage-induced velocity model is based on parameter identification of isolated fuselage-induced velocity data (computed by means of computational fluid dynamics, CFD) and is intended to be generic enough to be used for real helicopter fuselages as well. The accuracies obtained in reproducing the CFD data show a remaining average error of less or equal 5 of the peak-to-peak induced velocity range, which is considered sufficient for comprehensive code analysis.

HART↗

Hypersonic Chemically Reacting Boundary-Layer Stability using LASTRAC

Hypersonic boundary layer transition is critical to the design of all hypersonic vehicles due to its effect on the heat transfer into the vehicle surface and potential drag enhancement or reduction during reentry. Boundary layer transition and boundary layer stability analysis under hypersonic conditions has been studied for decades, yet there is ample room for improved accuracy and further investigations into the relevant phenomena. In this work, we present a recent implementation of chemical equilibrium, finite-rate chemistry, and thermochemical nonequilibrium capabilities into LASTRAC, an existing well-established boundary-layer stability analysis code. Verification against existing numerical results in the literature are presented. LASTRAC was previously able to address calorically perfect flows. By using solutions of the Parabolized Stability Equations (PSE) with chemical and thermal nonequilibrium, we are able to investigate the effects of chemical and thermal nonequilibrium on a variety of phenomena including stationary crossflow instability on a swept wing and 2nd mode instabilities over a wedge.

Kline, H. L.↗

Computational Study of the Side-By-Side Urban Air Taxi Concept

High-fidelity computational fluid dynamics simulations of NASA’s Side-by-side air taxi concept have been carried out. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids using high-order accurate schemes, dual-time stepping, and a hybrid turbulence model. The flow solver has been loosely coupled with a helicopter comprehensive analysis code in order to get the trimmed flight solution. The vehicle simulated is a six-passenger side-by-side intermeshing rotor helicopter with hybrid propulsion for air taxi operations, also known as urban air mobility applications. This concept vehicle is intended to focus and guide NASA research activities in support of aircraft development for emerging aviation markets, in particular vertical take-off and landing air taxi operations.

Ventura Diaz, Patricia↗

Estimating Drag and Heating Coefficients for Hollow Reentry Objects in Transitional Flow Using DSMC

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These methods are typically limited to convex, blunt objects (such as spheres) and are applied to other objects such as boxes and cylinders using multiplicative “shape factors” to account for the different behavior. Previous literature has analyzed the aerodynamic and aerothermodynamic properties of flow around sharp-edged objects like boxes and cylinders in transitional flow, though only those objects with solid external boundaries. However, many reentry objects we have encountered in real spacecraft have been hollow (i.e., with the potential to allow flow through them). We present here preliminary results from analyses performed using the NASA Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) on hollow cylinders and boxes (with varying wall thickness-diameter ratio).

Marichalar, Jeremiah J.↗

Computational Study of NASA's Quadrotor Urban Air Taxi Concept

High-fidelity computational fluid dynamics simulations have been carried out in order to analyze NASA's quadrotor urban air taxi concept for urban air mobility, also know as on-demand mobility applications. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. The ow solver has been loosely coupled with a rotorcraft comprehensive analysis code. The vehicle simulated is a six-passenger quadrotor for air taxi operations. A study of power reduction as a function of the rear-rotor to front-rotors vertical separation has been performed, for a quad-rotor without the airframe, in cruise flight conditions. Then, the quadrotor without the airframe has been simulated in hover. The airloads and wake geometries are analyzed. To finish the study the complete quadrotor vehicle is presented. NASA's quadrotor air taxi concept is one of the many concepts being developed by NASA in support of aircraft development for vertical take-o and landing air taxi operations.

Ventura Diaz, Patricia↗

Evolving Material Porosity on an Additive Manufacturing Simulation with the Generalized Method of Cells

The effect of material porosity on final part distortion and residual stresses in a selective laser sintering manufacturing simulation is presented here. A time-dependent thermomechanical model is used with the open-source FEA software CalculiX. Effective homogenized material properties for Inconel 625 are precomputed using NASA’s Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC). The evolving porosity of the material is estimated with each pass of the laser beam during simulation runtime. A comparison with a homogenous model and the evolving model shows that the evolving porous model predicts larger distortions with greater residual stresses.

Silva, Luis F.↗

Multiple Boundary Layer Instability Modes with Nonequilibrium and Wall Temperature Effects Using LASTRAC

Prediction and control of boundary layer transition from laminar to turbulent is important to many flow regimes and vehicle designs, including vehicles operating at hypersonic conditions where nonequilibrium effects may be encountered. Wall cooling is known to affect the instability characteristics of the boundary layer and subsequently the transition location. Design considerations, including material failure and fuel chemistry, require the use of actively cooled walls in hypersonic vehicles, further motivating the study of wall temperature effects on top of the considerations of reducing heat flux, drag, and uncertainty. In this work, we analyze the stability of a boundary layer with chemical and thermal nonequilibrium on a Mach 20, 6 wedge. We investigate the effects of wall temperature on multiple unstable modes individually and on the integrated growth of disturbances along the surface. We use the LAngley Stability and TRansition Analysis Code (LASTRAC) to evaluate boundary layer stability, using capabilities implemented by the authors. Included are results that address chemical nonequilibrium with both thermal equilibrium and nonequilibrium.

Kline, H. L.↗

Effect of Material Porosity on Residual Stress in an Additive Manufacturing Simulation using the Generalized Method of Cells

The effect of material porosity on final part distortion and residual stresses in a selective laser sintering manufacturing simulation is presented here. A time-dependent thermo-mechanical model is used with the open-source FEA software CalculiX. Effective homogenized material properties for Inconel 625 are precomputed using NASA’s Micromechanics Analysis Code with the Generalized Method of Cells. The evolving porosity of the material is estimated with each pass of the laser beam during simulation runtime. A comparison with a homogenous model and the evolving model shows that the model with evolving porosity predicts larger distortions with greater residual stresses.

Additive Manufactuirng↗

Leveraging STARE for Co-aligned Data Locality with netCDF and Python MPI

We have leveraged STARE indexing to package partitioned data chunks from diverse datasets into netCDF files, distributed them on a cluster of 16 lightweight nodes with their placements spatiotemporally co-aligned, and demonstrated a few integrative analyses using netCDF parallel I/O and Python MPI, with single-user performance and scalability comparable to, or even better than, that of a parallel array database management system (ADBMS) such as SciDB. However, records of the node location and STARE index ranges for each data chunk, similar to the chunk maps of SciDB, must be maintained and consulted by the I/O and analysis code for coordinating the analytic operations in parallel, in order to achieve the good performance and scalability.

Kwo-Sen Kuo↗

EnDEVR: An Environment for Data Engineering in VR

Organizations within the public and private sectors are looking to improve their data science operations for research, development, and operational purposes. As this interest in data science grows, so too does interest in tools and programs that facilitating such operations. This paper presents EnDEVR, the Environment for Data Engineering in Virtual Reality, a user-extensible system that allows for the execution of custom data science functions. We describe the user-centered design process for the EnDEVR system, based on the needs of subject matter experts within the aerospace and aviation industries. Fundamental to this design is a method for users to integrate their existing analysis code (that may leverage third-party vendors) as usable objects in the virtual environment; this allows organizations to operate without the price tag and lifetime obligations associated with commercial data science toolkits, a concept known as “vendor lock-in”. We then present our system design with iterative implementations, and conclude with lessons learned. The presented work may provide insight into future user-centered design processes for VR applications, particularly for digital transformation stakeholders to circumvent vendor lock-in.

Human-centered computing↗

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗