Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimensional 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 1,207 records · Page 67

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Towne, Charles E.↗

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

Non-equilibrium one-dimensional two-phase flow in variable area channels

A one-dimensional nonequilibrium flow analysis has been formulated for a one component two phase flow. The flow is considered homogeneous and essentially isothermal. Phase change is assumed to occur at heterogeneous nucleation sites and the growth of the vapor bubbles is governed by heat conduction from the liquid to the bubble. The analysis adjusted for friction is applied to liquid nitrogen flow in a venturi and comparison is made with the NASA experimental results of Simoneau. Good agreement with the experiments is obtained when one assumes the effective activation energy for nucleus formation to be small but nonzero. The computed pressure distributions deviate from the experimental results in the throat region of the venturi in a manner consistent with centrifugal effects not accounted for in the one-dimensional theory. The results are shown to depend not only on cavitation number but on additional dimensionless parameters governing the nonequilibrium production and subsequent growth of nuclei.

Rohatgi, U. S.↗

Vadose and Saturated Zone Flow and Transport Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

The purpose of the fate and transport modeling described in this environmental calculation file (ECF) is to evaluate the impacts to groundwater associated with waste disposal operations at Mixed Waste Low-Level Burial Ground (LLBG) Trenches 31 and 34 to satisfy requirements in DOE O 435.1, Radioactive Waste Management. The model integrates the flow and transport in the vadose zone beneath the active trenches with the saturated zone downgradient of the trenches to predict the radionuclide concentration at the point of assessment (POA). DOE M 435.1-1, Radioactive Waste Management Manual, defines the POA as the point of highest projected dose or concentration beyond a 100 m (328 ft) buffer zone surrounding the disposed waste. The modeling is conducted in accordance with the DOE G 435.1, Implementation Guide for Use with DOE M 435.1-1, performance assessment (PA) guidelines. The modeling involves evaluation of the groundwater concentrations and radionuclide arrival times during the 1,000-year compliance and 10,000-year sensitivity-uncertainty periods per DOE O 435.1 and DOE M 435.1. This analysis does not consider radionuclide release during facility operations, only the post-closure impacts of the radionuclides to the environment. The evaluation of potential radiological dose to groundwater receptors caused by releases from a closed facility containing radioactive waste typically includes the following: (1) Release of radionuclides from that facility (2) Transport of those radionuclides through the environment, and (3) Exposure to humans to environmental concentration levels of those radionuclides The fate and transport three-dimensional (3D) model analysis involves the post-closure impacts to the environment of the technetium-99, iodone-129, and uranium (all isotopes in the waste). The residual inventory estimates include several radionuclides, but technetium-99 is typically responsible for almost all of the beta-gamma dose equivalent associated with groundwater (water resources) protection per 40 CFR 141, “National Primary Drinking Water Regulations” (e.g., see the results in WCH-520, Performance Assessment of Environmental Restoration Disposal Facility, Hanford Site, Washington; hereinafter referred to as the ERDF PA), and iodone-129 can also be a significant dose contributor for some waste (e.g., RPP-RPT-59958, Performance Assessment for the Integrated Disposal Facility, Hanford Site, Washington; hereinafter referred to as the IDF PA). Uranium does not typically factor significantly into the impacts to groundwater, even during the 10,000-year sensitivity-uncertainty period, but always remains of interest as a contaminant. This ECF does not address vadose and saturated zone modeling for Trench 94 of the 200 East Area LLBG. Current information confirms the validity of the low corrosion rate of the naval reactor plant carbon steel (HY-80), and the even lower corrosion rate of the nickel-iron-chromium alloy reactor vessel (Inconel Alloy 600) presented in DOE/EIS-0259, Final Environmental Impact Statement on the Disposal of Decommissioned, Defueled Cruiser, Ohio Class, and Los Angeles Class Naval Reactor Plants. Based on these low corrosion rates, the time to breach the reactor vessel to allow release of radionuclides from the activated metal of the reactor vessel internal structure is at least 10,000 years. This time to breach precludes the need to evaluate the vadose and saturated zone transport of contaminants released from the reactor compartment disposal packages in the 200 East Area LLBGs PA (CP-63826, Waste Release Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Views on the impact of HOST

The Hot Section Technology (HOST) Project, which was initiated by NASA Lewis Research Center in 1980 and concluded in 1987, was aimed at improving advanced aircraft engine hot section durability through better technical understanding and more accurate design analysis capability. The project was a multidisciplinary, multiorganizational, focused research effort that involved 21 organizations and 70 research and technology activities and generated approximately 250 research reports. No major hardware was developed. To evaluate whether HOST had a significant impact on the overall aircraft engine industry in the development of new engines, interviews were conducted with 41 participants in the project to obtain their views. The summarized results of these interviews are presented. Emphasis is placed on results relative to three-dimensional inelastic structural analysis, thermomechanical fatigue testing, constitutive modeling, combustor aerothermal modeling, turbine heat transfer, protective coatings, computer codes, improved engine design capability, reduced engine development costs, and the impacts on technology transfer and the industry-government partnership.

J B Esgar↗

Icarus Application to Dragonfly Heatshield

Icarus, an in-house material response solver developed at NASA Ames, is applied to perform 2-D bondline-temperature analysis on the heatshield surface of the latest Dragonfly entry system design for a mission to Titan. Two axisymmetric domains are used: the near-shoulder region and the entire heatshield. The heatshield consists of a PICA TPS bonded to a layered system of Aluminum honeycomb wrapped by an M55J carbon-fiber facesheet. The backshell, which is also included for more accurate analysis, consists of similar materials layered. Icarus simulations using both domains indicate that the in-plane thermal conductivity property of the M55J facesheet is observed to play a dominant role in bondline temperature. Similarly, the maximum bondline temperature is found around the PICA-tile-interface region instead of the near-shoulder region or the stagnation point based on the current trajectory using the orthotropic properties of M55J. Therefore, these findings indicate the significance of 2-D or higher-dimensional material modeling analysis in fully understanding how bondline temperature behaves across the heatshield and obtaining a basis for TPS design locations. Verification at different radial locations ahead of the shoulder, where thermal conduction is close to 1-D, between Icarus (2-D), Icarus (1-D), and FIAT is performed and good agreement is observed. With the entire heatshield domain, the total mass loss of the aeroshell materials is approximately 3% of their initial mass potentially due to thermal decomposition of phenolic resins inside PICA under the current high-heating environment.

TPS material modeling↗

Supersonic Flow of Chemically Reacting Gas-Particle Mixtures. Volume 2: RAMP - A Computer Code for Analysis of Chemically Reacting Gas-Particle Flows

A computer program written in conjunction with the numerical solution of the flow of chemically reacting gas-particle mixtures was documented. The solution to the set of governing equations was obtained by utilizing the method of characteristics. The equations cast in characteristic form were shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The characteristic directions for the gas-particle system are found to be the conventional gas Mach lines, the gas streamlines and the particle streamlines. The basic mesh construction for the flow solution is along streamlines and normals to the streamlines for axisymmetric or two-dimensional flow. The analysis gives detailed information of the supersonic flow and provides for a continuous solution of the nozzle and exhaust plume flow fields. Boundary conditions for the flow solution are either the nozzle wall or the exhaust plume boundary.

Penny, M. M.↗

Elastic-plastic analysis of a propagating crack under cyclic loading

Development and application of a two-dimensional finite-element analysis to predict crack-closure and crack-opening stresses during specified histories of cyclic loading. An existing finite-element computer program which accounts for elastic-plastic material behavior under cyclic loading was modified to account for changing boundary conditions - crack growth and intermittent contact of crack surfaces. This program was subsequently used to study the crack-closure behavior under constant-amplitude and simple block-program loading.

Newman, J. C., Jr.↗

Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP

A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗

Energetics analysis of the observed and simulated general circulation using three-dimensional normal mode expansions

The energetics characteristics of the observed and simulated general circulation are analyzed using three-dimensional normal mode expansions. The data sets involved are the Goddard Laboratory for Atmospheres (GLA) analysis and simulation data and the Geophysical Fluid Dynamics Laboratory (GFDL) analysis data. The spectral energy properties of the Rossby and gravity modes and energy transformations are presented. Significant influences of model characteristics and the assimilation techniques are observed in the barotropic energy spectrum, particularly for the gravity mode. Energy transformations of the zonal mean field in the GLA analysis and simulation are similar, but distinctly different from that in the GFDL analysis. However, overall, the energy generation in the baroclinic mode is largely balanced by the sink in the barotropic mode. The present study may demonstrate utilities of the three-dimensional normal mode energetics in the analysis of the general circulation.

Tanaka, Hiroshi↗

Experimental observations and finite element analysis of the initiation of fiber microbuckling in notched composite laminates

A better understanding of the factors that affect the semi-circular edge-notched compressive strength is developed, and the associated failure mode(s) of thermoplastic composite laminates with multidirectional stacking sequences are identified. The primary variables in this investigation are the resin nonlinear shear constitutive behavior, stacking sequence (orientation of plies adjacent to the 0 degree plies), resin-rich regions between the 0 degree plies and the off-axis supporting plies, fiber/matrix interfacial bond strength, and initial fiber waviness. Two thermoplastic composite material systems are used in this investigation. The materials are the commercial APC-2 (AS4/PEEK) and a poor interface experimental material, AU4U/PEEK, designed for this investigation. Notched compression specimens are studied at 21, 77, and 132 C. Geometric and material nonlinear two-dimensional finite element analysis is used to model the initiation of fiber microbuckling of both the ideal straight fiber and the more realistic initially wavy fiber. The effects of free surface, fiber constitutive properties, matrix constitutive behavior, initial fiber curvature, and fiber/matrix interfacial bond strength on fiber microbuckling initiation strain levels are considered.

Guynn, E. Gail↗

Coupled multi-disciplinary composites behavior simulation

The capabilities of the computer code CSTEM (Coupled Structural/Thermal/Electro-Magnetic Analysis) are discussed and demonstrated. CSTEM computationally simulates the coupled response of layered multi-material composite structures subjected to simultaneous thermal, structural, vibration, acoustic, and electromagnetic loads and includes the effect of aggressive environments. The composite material behavior and structural response is determined at its various inherent scales: constituents (fiber/matrix), ply, laminate, and structural component. The thermal and mechanical properties of the constituents are considered to be nonlinearly dependent on various parameters such as temperature and moisture. The acoustic and electromagnetic properties also include dependence on vibration and electromagnetic wave frequencies, respectively. The simulation is based on a three dimensional finite element analysis in conjunction with composite mechanics and with structural tailoring codes, and with acoustic and electromagnetic analysis methods. An aircraft engine composite fan blade is selected as a typical structural component to demonstrate the CSTEM capabilities. Results of various coupled multi-disciplinary heat transfer, structural, vibration, acoustic, and electromagnetic analyses for temperature distribution, stress and displacement response, deformed shape, vibration frequencies, mode shapes, acoustic noise, and electromagnetic reflection from the fan blade are discussed for their coupled effects in hot and humid environments. Collectively, these results demonstrate the effectiveness of the CSTEM code in capturing the coupled effects on the various responses of composite structures subjected to simultaneous multiple real-life loads.

Singhal, Surendra N.↗

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries.

Huff, Dennis L.↗

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries.

Huff, Dennis L.↗

Low-cycle fatigue analysis of a cooled copper combustion chamber

A three-dimensional finite element elastoplastic strain analysis was performed for the throat section of regeneratively cooled rocket engine combustion chamber. The analysis included thermal and pressure loads, and the effects of temperature dependent material properties, to determine the strain range corresponding to the engine operating cycle. The strain range was used in conjunction with OFHC copper isothermal fatigue test data to predict engine low-cycle fatigue life. The analysis was performed for chamber configuration and operating conditions corresponding to a hydrogen-oxygen chamber which was fatigue tested to failure at the NASA Lewis Research Center.

Miller, R. W.↗

Nonlinear evolution of interacting oblique waves on two-dimensional shear layers

The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Goldstein, M. E.↗