Search NASA⌕ Search

SEARCH · Search NASA

Results for “benchmark and validation”

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 145 records · Page 8

Summary of Shock Wave Turbulent Boundary Layer Interaction Experiments in a Circular Test Section

A series of experiments were performed at Mach 2.5 in a 17 cm diameter circular test section to characterize an impinging/reflected shock wave turbulent boundary layer interaction generated by a cone-cylinder centerbody. The cone-cylinder centerbody generates a conical shock wave that interacts with the naturally occurring boundary layer developing on the test section wall. Three different cone angles were used in the experiment to study unseparated, incipiently separated, and separated interactions. When the cone-cylinder centerbody is positioned on the centerline, a flowfield which is two-dimensional in the mean is generated. Three dimensional interactions were also created by offsetting the cone-cylinder centerbody from the test section centerline. The results are intended to provide benchmark quality datasets for computational fluid dynamics (CFD) validation without the pitfalls inherent in rectangular configurations where corner effects prohibit a truly two-dimensional flow in the mean. The experimental measurements included surface flow visualization, wall static pressure, flowfield Pitot tube pressure, constant-voltage anemometry (CVA) normal hot-wire, and particle image velocimetry (PIV) measurements. The hot-wire measurements were used to calculate mean mass flux and total temperature profiles, mass flux and total temperature turbulence intensities, and the mass flux-total temperature correlation. The PIV measurements provide three-dimensional mean velocity measurements. Agreement between the pressure, hot-wire, and PIV measurements is established in the undisturbed upstream flowfield.

compressible flow↗

Mastering the Cryogenic Frontier: Predicting Sloshing in Cryogenic Propellant Tanks

In recent years, NASA has undertaken the development and validation of computation fluid dynamic (CFD) models for predicting cryogenic two-phase sloshing within propellant tanks. These models have been rigorously anchored against cryogenic ground tests that induce sloshing, as well as limited microgravity tests. The primary focus of these models is to accurately predict the two-phase condensation and, at times, evaporation within the cryogenic propellant. This approach is important for estimating the required autogenous and/or helium pressurant for completing propulsion system burns or ensuring adequate pressure for in-space propellant transfers. Unlike in hypergolic propulsion systems, where the propellant is stored at room temperature, cryogenic propellant is highly volatile. Therefore, ensuring the accurate estimation of pressurant necessitates the precise modeling of two-phase interfacial heat and mass transfer processes. This paper aims to provide an overview of NASA’s continuing efforts in developing and validating these CFD models, with particular emphasis on benchmark cases with sloshing that enhances phase change, tracing their historical development, and demonstrating their predictive accuracy when compared to test results.

low-gravity fluid modeling↗

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗

Comprehensive 3D Multiphysics Model on Electrochemical Recovery of O 2 from Metabolic CO 2 at the International Space Station (ISS)

The International Space Station (ISS) is presently equipped with an elaborate, heavy, and high-power consuming system that recovers approximately 50% of O 2 from metabolic CO 2 as part of the atmospheric revitalization (AR) at the ISS habitat. Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery to reach the self-sufficiency required for long space missions beyond earth’s low orbit. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to not only increase significantly current O 2 recovery efficiency, improving self-sufficiency on AR at the ISS habitat and future long missions, but also reduce the complexity of the system. The authors have developed and deployed a comprehensive 3D multiphysics model that thoroughly replicates the actual configuration and fluid/material domains of the MFECR. The coupled physics in this multiphysics model include multicomponent-multiphase electrochemical-driven reactions, non-ideal mass transport mechanism, free and porous flow, heat transfer, CO 2 solubility on alkaline electrolyte, water condensation on porous medium, and DC electrical current generation along with Joule heating effect. This model is aimed to conduct quantitive benchmark on three different MFECR’s layouts, one without serpentine paths (plain) and two with serpentines leading to four and twelve paths respectively. Once experimental data is generated via a test matrix of 200 tests, the model will be validated to conduct MFECR’s process optimization and revalidate the quantitive benchmark on three different MFECR’s layouts.

Jesus A. Dominguez↗

A benchmark for galactic cosmic ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

A benchmark for galactic cosmic-ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

Investigation of the Flow Physics Driving Stall-Side Flutter in Advanced Forward Swept Fan Designs

Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.

Sanders, Albert J.↗

Buckling and Damage Resistance of Transversely-Loaded Composite Shells

Experimental and numerical work was conducted to better understand composite shell response to transverse loadings which simulate damage-causing impact events. The quasi-static, centered, transverse loading response of laminated graphite/epoxy shells in a [+/-45(sub n)/O(sub n)](sub s) layup having geometric characteristics of a commercial fuselage are studied. The singly-curved composite shell structures are hinged along the straight circumferential edges and are either free or simply supported along the curved axial edges. Key components of the shell response are response instabilities due to limit-point and/or bifurcation buckling. Experimentally, deflection-controlled shell response is characterized via load-deflection data, deformation-shape evolutions, and the resulting damage state. Finite element models are used to study the kinematically nonlinear shell response, including bifurcation, limit-points, and postbuckling. A novel technique is developed for evaluating bifurcation from nonlinear prebuckling states utilizing asymmetric spatial discretization to introduce numerical perturbations. Advantages of the asymmetric meshing technique (AMT) over traditional techniques include efficiency, robustness, ease of application, and solution of the actual (not modified) problems. The AMT is validated by comparison to traditional numerical analysis of a benchmark problem and verified by comparison to experimental data. Applying the technique, bifurcation in a benchmark shell-buckling problem is correctly identified. Excellent agreement between the numerical and experimental results are obtained for a number of composite shells although predictive capability decreases for stiffer (thicker) specimens which is attributed to compliance of the test fixture. Restraining the axial edge (simple support) has the effect of creating a more complex response which involves unstable bifurcation, limit-point buckling, and dynamic collapse. Such shells were noted to bifurcate into asymmetric deformation modes but were undamaged during testing. Shells in this study which were damaged were not observed to bifurcate. Thus, a direct link between bifurcation and atypical damage could not be established although the mechanism (bifurcation) was identified. Recommendations for further work in these related areas are provided and include extensions of the AMT to other shell geometries and structural problems.

Wardle, Brian L.↗

Development of a Linearized Unsteady Euler Analysis with Application to Wake/Blade-Row Interactions

A three-dimensional, linearized, Euler analysis is being developed to provide a comprehensive and efficient unsteady aerodynamic analysis for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. The mathematical models needed to describe nonlinear and linearized, inviscid, unsteady flows through a blade row operating within a cylindrical annular duct are presented in this report. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to far-field eigen analyses, is also described. The linearized aerodynamic and numerical models have been implemented into the three-dimensional unsteady flow code, LINFLUX. This code is applied herein to predict unsteady subsonic flows driven by wake or vortical excitations. The intent is to validate the LINFLUX analysis via numerical results for simple benchmark unsteady flows and to demonstrate this analysis via application to a realistic wake/blade-row interaction. Detailed numerical results for a three-dimensional version of the 10th Standard Cascade and a fan exit guide vane indicate that LINFLUX is becoming a reliable and useful unsteady aerodynamic prediction capability that can be applied, in the future, to assess the three-dimensional flow physics important to blade-row, aeroacoustic and aeroelastic responses.

Verdon, Joseph M.↗

Validation of CFD/Heat Transfer Software for Turbine Blade Analysis

I am an intern in the Turbine Branch of the Turbomachinery and Propulsion Systems Division. The division is primarily concerned with experimental and computational methods of calculating heat transfer effects of turbine blades during operation in jet engines and land-based power systems. These include modeling flow in internal cooling passages and film cooling, as well as calculating heat flux and peak temperatures to ensure safe and efficient operation. The branch is research-oriented, emphasizing the development of tools that may be used by gas turbine designers in industry. The branch has been developing a computational fluid dynamics (CFD) and heat transfer code called GlennHT to achieve the computational end of this analysis. The code was originally written in FORTRAN 77 and run on Silicon Graphics machines. However the code has been rewritten and compiled in FORTRAN 90 to take advantage of more modem computer memory systems. In addition the branch has made a switch in system architectures from SGI's to Linux PC's. The newly modified code therefore needs to be tested and validated. This is the primary goal of my internship. To validate the GlennHT code, it must be run using benchmark fluid mechanics and heat transfer test cases, for which there are either analytical solutions or widely accepted experimental data. From the solutions generated by the code, comparisons can be made to the correct solutions to establish the accuracy of the code. To design and create these test cases, there are many steps and programs that must be used. Before a test case can be run, pre-processing steps must be accomplished. These include generating a grid to describe the geometry, using a software package called GridPro. Also various files required by the GlennHT code must be created including a boundary condition file, a file for multi-processor computing, and a file to describe problem and algorithm parameters. A good deal of this internship will be to become familiar with these programs and the structure of the GlennHT code. Additional information is included in the original extended abstract.

Kiefer, Walter D.↗

Benchmark Wall Heat Flux Data for a GO2/GH2 Single Element Injector

Computational Fluid Dynamics (CFD) is becoming an important component of injector design in the rocket industry. Injector designers who use CFD in the design process need to understand the accuracy level of the particular code being used for certain aspects of the design. This paper presents a recent effort to acquire benchmark quality data to be used for CFD code validation. Detailed chamber wall temperature and heat flux data was acquired for a gaseous oxygen, gaseous hydrogen single element shear coaxial injector in a 1.5 inch diameter copper heat sink chamber at Penn State University. The data was taken using both coaxial and water cooled heat flux gauges. Tests were run using hot gases generated from both fuel and oxidizer preburners. Tests were conducted over a chamber pressure range of 300 to 750 psia. Data analysis and uncertainty information will also be presented.

Pal, Sibtosh↗

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The influence of bulk turbulence and boundary layer thickness on the secondary flow development in a square, 90 degree turning duct was investigated. A three-dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity at six cross-planes in the duct. The results from this investigation, with entrance boundary layer thickness of 20 percent, were compared with the thin boundary layer results documented in NASA CR-174811. The axial velocity profiles, cross-flow velocities, and turbulence intensities were compared and evaluated with regard to the influence of bulk turbulence intensity and boundary layer thickness, and the influence was significant. The results of this investigation expand the 90 degree curved duct experimental data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, R. A.↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, Roger A.↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Shohadaee, Ahmad A.↗

Experimental and computational surface and flow-field results for an all-body hypersonic aircraft

Personnel from NASA Ames Research Center presented a paper on establishing a benchmark experimental data base for generic hypersonic vehicle shape for validation and/or calibration of advanced computational fluid dynamics computer codes. The need for this capability is based on a requirement for extensive hypersonic data to fully validate CFD codes to be used for NASP and other hypersonic vehicles. The use of wind tunnel models in the Ames 3.5-ft Hypersonic Wind Tunnel to obtain pertinent surface and flow-field data over a broad range of test conditions is described.

Lockman, William K.↗

Development of Constraint Force Equation Methodology for Application to Multi-Body Dynamics Including Launch Vehicle Stage Seperation

The objective of this report is to develop and implement a physics based method for analysis and simulation of multi-body dynamics including launch vehicle stage separation. The constraint force equation (CFE) methodology discussed in this report provides such a framework for modeling constraint forces and moments acting at joints when the vehicles are still connected. Several stand-alone test cases involving various types of joints were developed to validate the CFE methodology. The results were compared with ADAMS(Registered Trademark) and Autolev, two different industry standard benchmark codes for multi-body dynamic analysis and simulations. However, these two codes are not designed for aerospace flight trajectory simulations. After this validation exercise, the CFE algorithm was implemented in Program to Optimize Simulated Trajectories II (POST2) to provide a capability to simulate end-to-end trajectories of launch vehicles including stage separation. The POST2/CFE methodology was applied to the STS-1 Space Shuttle solid rocket booster (SRB) separation and Hyper-X Research Vehicle (HXRV) separation from the Pegasus booster as a further test and validation for its application to launch vehicle stage separation problems. Finally, to demonstrate end-to-end simulation capability, POST2/CFE was applied to the ascent, orbit insertion, and booster return of a reusable two-stage-to-orbit (TSTO) vehicle concept. With these validation exercises, POST2/CFE software can be used for performing conceptual level end-to-end simulations, including launch vehicle stage separation, for problems similar to those discussed in this report.

Pamadi, Bandu N.↗

Benchmarking the QUAD4/TRIA3 element

The QUAD4 and TRIA3 elements are the primary plate/shell elements in NASTRAN. These elements enable the user to analyze thin plate/shell structures for membrane, bending and shear phenomena. They are also very new elements in the NASTRAN library. These elements are extremely versatile and constitute a substantially enhanced analysis capability in NASTRAN. However, with the versatility comes the burden of understanding a myriad of modeling implications and their effect on accuracy and analysis quality. The validity of many aspects of these elements were established through a series of benchmark problem results and comparison with those available in the literature and obtained from other programs like MSC/NASTRAN and CSAR/NASTRAN. Never-the-less such a comparison is never complete because of the new and creative use of these elements in complex modeling situations. One of the important features of QUAD4 and TRIA3 elements is the offset capability which allows the midsurface of the plate to be noncoincident with the surface of the grid points. None of the previous elements, with the exception of bar (beam), has this capability. The offset capability played a crucial role in the design of QUAD4 and TRIA3 elements. It allowed modeling layered composites, laminated plates and sandwich plates with the metal and composite face sheets. Even though the basic implementation of the offset capability is found to be sound in the previous applications, there is some uncertainty in relatively simple applications. The main purpose of this paper is to test the integrity of the offset capability and provide guidelines for its effective use. For the purpose of simplicity, references in this paper to the QUAD4 element will also include the TRIA3 element.

Pitrof, Stephen M.↗