Search NASA⌕ Search

SEARCH · Search NASA

Results for “control volume finite element”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Water-cooled Helmholtz coil for Ultra-low Magnetic Field Testing

Ultra-low magnetic field testing requires precise and finely calibrated instrumentation interacting with very marginal magnetic fields in highly controlled environments. We designed a large Helmholtz coil system for implementation around a 17.825-inch outer diameter Helium-3 Cryogenic fridge, to serve as both a field-zeroing device and to enable higher resolution on magnetic field sweeps during experiments. We present finite volume modeling and finite element analysis results which indicate the device could safely run at 20 [A] in perpetuity, producing a magnetic field of 65 [mT], without the surface facing the Cryogenic fridge exceeding 50°C. The overall system parts expense is $\$$4,882.76.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

International Congress on Recent Developments in Air- and Structure-Borne Sound and Vibration, 2nd, Auburn Univ., AL, Mar. 4-6, 1992, Proceedings. Vols. 1-3

The papers presented in this volume cover a variety of topics, including sound intensity, structural intensity, modal analysis and synthesis, statistical energy analysis and energy methods, passive and active damping, and boundary element methods. Attention is also given to diagnostics and condition monitoring, material characterization and nondestructive evaluation, active noise and vibration control, sound radiation and scattering, and finite element analysis.

Crocker, Malcolm J.↗
DOE OSTI2020

Shape Optimization for Control and Isolation of Structural Vibrations in Aerospace and Defense Applications

Among the main challenges in shape optimization is the coupling of Finite Element Method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multi-physics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. There are two mathematically equivalent approaches to computing the shape derivative: the volume method, and the boundary method. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. Prior implementations of shape derivatives at Sandia have been based on the volume method. We introduce the strip method, which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. The development of the strip method also offers us the opportunity to share some lessons learned about implementing the volume method and boundary method, to show shape optimization results on problems of interest, and to begin addressing the other main challenges at hand: constraints on optimized shapes, and their interplay with optimization algorithms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

A Summary of the Development of a Nominal Land Landing Airbag Impact Attenuation System for the Orion Crew Module

Airborne Systems North America (formally Irvin Aerospace Inc) has developed an Airbag Landing System for the Orion Crew Module of the Crew Exploration Vehicle. This work is in support of the NASA Langley Research Center Landing System Advanced Development Project. Orion is part of the Constellation Program to send human explorers back to the moon, and then onwards to Mars and other destinations in the Solar System. A component of the Vision for Space Exploration, Orion is being developed to also enable access to space following the retirement of the Space Shuttle in the next decade. This paper documents the development of a conceptual design, fabrication of prototype assemblies, component level testing and two generations of airbag landing system testing. The airbag system has been designed and analyzed using the transient dynamic finite element code LS-DYNA(RegisteredTradeMark). The landing system consists of six airbag assemblies; each assembly comprising a primary impact venting airbag and a non-venting anti-bottoming airbag. The anti-bottoming airbag provides ground clearance following the initial impact attenuation sequence. Incorporated into each primary impact airbag is an active vent that allows the entrapped gas to exit the control volume. The size of the vent is tailored to control the flow-rate of the exiting gas. An internal shaping structure is utilized to control the shape of the primary or main airbags prior to ground impact; this significantly improves stroke efficiency and performance.

Tutt, Ben↗
NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗
NASA NTRS1991

Slave finite element for non-linear analysis of engine structures. Volume 2: Programmer's manual and user's manual

The programming aspects of SFENES are described in the User's Manual. The information presented is provided for the installation programmer. It is sufficient to fully describe the general program logic and required peripheral storage. All element generated data is stored externally to reduce required memory allocation. A separate section is devoted to the description of these files thereby permitting the optimization of Input/Output (I/O) time through efficient buffer descriptions. Individual subroutine descriptions are presented along with the complete Fortran source listings. A short description of the major control, computation, and I/O phases is included to aid in obtaining an overall familiarity with the program's components. Finally, a discussion of the suggested overlay structure which allows the program to execute with a reasonable amount of memory allocation is presented.

Witkop, D. L.↗
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
NASA NTRS1988

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗
NASA NTRS1980

Steady, oscillatory, and unsteady subsonic Aerodynamics, production version 1.1 (SOUSSA-P1.1). Volume 2: User/programmer manual

A user/programmer manual for the computer program SOUSSA P 1.1 is presented. The program was designed to provide accurate and efficient evaluation of steady and unsteady loads on aircraft having arbitrary shapes and motions, including structural deformations. These design goals were in part achieved through the incorporation of the data handling capabilities of the SPAR finite element Structural Analysis computer program. As a further result, SOUSSA P possesses an extensive checkpoint/ restart facility. The programmer's portion of this manual includes overlay/subroutine hierarchy, logical flow of control, definition of SOUSSA P 1.1 FORTRAN variables, and definition of SOUSSA P 1.1 subroutines. Purpose of the SOUSSA P 1.1 modules, input data to the program, output of the program, hardware/software requirements, error detection and reporting capabilities, job control statements, a summary of the procedure for running the program and two test cases including input and output and listings are described in the user oriented portion of the manual.

Smolka, S. A.↗
NASA NTRSDate not supplied

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗
NASA NTRS2010

Acoustic-Modal Testing of the Ares I Launch Abort System Attitude Control Motor Valve

The Attitude Control Motor (ACM) is being developed for use in the Launch Abort System (LAS) of NASA's Ares I launch vehicle. The ACM consists of a small solid rocket motor and eight actuated pintle valves that directionally allocate.thrust_- 1t.has-been- predicted-that significant unsteady. pressure.fluctuations.will.exist. inside the-valves during operation. The dominant frequencies of these oscillations correspond to the lowest several acoustic natural frequencies of the individual valves. An acoustic finite element model of the fluid volume inside the valve has been critical to the prediction of these frequencies and their associated mode shapes. This work describes an effort to experimentally validate the acoustic finite model of the valve with an acoustic modal test. The modal test involved instrumenting a flight-like valve with six microphones and then exciting the enclosed air with a loudspeaker. The loudspeaker was configured to deliver broadband noise at relatively high sound pressure levels. The aquired microphone signals were post-processed and compared to results generated from the acoustic finite element model. Initial comparisons between the test data and the model results revealed that additional model refinement was necessary. Specifically, the model was updated to implement a complex impedance boundary condition at the entrance to the valve supply tube. This boundary condition models the frequency-dependent impedance that an acoustic wave will encounter as it reaches the end of the supply tube. Upon invoking this boundary condition, significantly improved agreement between the test data and the model was realized.

Davis, R. Benjamin↗
NASA NTRS2019

Nanopore-Templated Silver Nanoparticle Arrays Photopolymerized in Zero-Mode Wavelengths

In situ fabrication of nanostructures within a solid-polymer electrolyte confined to subwavelength-diameter nanoapertures is a promising approach for producing nanomaterials for nanophotonic and chemical sensing applications. The solid-polymerelectrolyte can be patterned by lithographic photopolymerization of poly(ethyleneglycol) diacrylate (PEGDA)-based silver cation (Ag+)-containing polyelectrolyte. Here, we present a new method for fabricating nanopore-templated Ag nanoparticle (AgNP)arrays by in situ photopolymerization using a zero-mode waveguide (ZMW) array to simultaneously template embedded AgNPs and control the spatial distribution of the optical field used for photopolymerization. The approach starts with an array of nanopores fabricated by sequential layer-by-layer deposition and focused ion beam milling. These structures have an optically transparent bottom, allowing access of the optical radiation to the attoliter-volume ZMW region to photopolymerize a PEGDA monomer solution containing AgNPs and Ag+. The electric field intensity distribution is calculated for various ZMW optical cladding layer thicknesses using finite-element simulations, closely following the light-blocking efficiency of the optical cladding layer. The fidelity of the polyelectrolyte nanopillar pattern was optimized with respect to experimental conditions, including the presence or absence of Ag+ and AgNPs and the concentrations of PEGDA and Ag+. The self-templated approach for photo patterning high-resolution photolabile polyelectrolyte nanostructures directly within a ZMW array could lead to a new class of metamaterials formed by embedding metal nanoparticles within a dielectric in a well-defined spatial array.

Zero-mode waveguide↗
NASA NTRS1981

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗
NASA NTRS1976

Optimum vibrating beams with stress and deflection constraints

The fundamental frequency of vibration of an Euler-Bernoulli or a Timoshenko beam of a specified constant volume is maximized subject to the constraint that under a prescribed loading the maximum stress or maximum deflection at any point along the beam axis will not exceed a specified value. In contrast with the inequality constraint which controls the minimum cross-section, the present inequality constraints lead to more meaningful designs. The inequality constraint on stresses is as easily implemented as the minimum cross-section constraint but the inequality constraint on deflection uses a treatment which is an extension of the matrix partitioning technique of prescribing displacements in finite element analysis.

Kamat, M. P.↗
DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Space shuttle main engine definition (phase B). Volume 5: Valves and interconnects

The steady state thermodynamic cycle balance of the single preburner staged combustion engine, coupled with dynamic transient analyses, dictated in detail the location and requirements for each valve defined in this volume. Valve configuration selections were influenced by overall engine and vehicle system weight and failure mode determinations. Modulating valve actuators are external to the valve and are line replaceable. Development and satisfactory demonstration of a high pressure dynamic shaft seal has made this configuration practical. Pneumatic motor driven actuators that use engine pumped hydrogen gas as the working fluid are used. The helium control system is proposed as a module containing a cluster of solenoid actuated valves. The separable couplings and flanges are designed to assure minimum leakage with minimum coupling weight. The deflection of the seal surface in the flange is defined by finite element analysis that has been confirmed with test data. The seal design proposed has passed preliminary pressure cycling and thermal cycling tests.

Schultz, D. F.↗
NASA NTRS2020

Bioinspired multilayered cellular composites with enhanced energy absorption and shape recovery

Inspired by the multiscale configuration of the microstructure of cork, the paper describes the design, 3D printing, and evaluation of a new type of multilayered cellular composite (MCC) structure composed of hard brittle and soft flexible phases. The mechanical behavior of 3D printed MCC structures have been investigated both experimentally and numerically. The experiments show that the MCC structure absorbs four times the amount of energy of a conventional cellular configuration under compressive strains up to 70%. Finite element simulations and 2D digital image correlation (DIC) also show that the multilayered architecture provides a more uniform strain distribution and higher stress transfer efficiency, with a resulting progressive failure mode rather than a catastrophic one. Cyclic loading tests demonstrate that the MCC structure also possesses exceptional shape recoverability under compressive deformations up to 40%. These remarkable performance characteristics result from synergies between the properties of the two constituent materials and the chosen multilayered cellular microstructure. The soft phase, in particular, plays a pivotal role in absorbing elastic energy during loading and then releasing the stored energy while unloading. The volume fraction of the soft phase is also essential to control energy absorption and the transition of failure modes. The deformation mechanisms demonstrated here are robust and applicable to other architected cellular materials across multiple length scales and suggest new ways to design lightweight and high-resilience structural materials.

bioinspired composites↗
NASA NTRS2001

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗