SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A study is made of the effects of variation in the lamination and geometric parameters of multilayered composite cylinders on the accuracy of the static and vibrational responses predicted by eight modeling approaches, based on two-dimensional shear-deformation shell theories. The standard of comparison is taken to be the exact three-dimensional elasticity solutions, and the quantities compared include both the gross response characteristics and detailed, through-the-thickness, distributions of displacements, stresses and strain energy densities. Based on the numerical studies conducted, a predictor-corrector approach, used in conjunction with the first-order shear-deformation theory (with five displacement parameters in the predictor phase), appears to be the most effective among the eight modeling approaches considered. For multilayered orthotropic cylinders the response quantities obtained by the predictor-corrector approach are shown to be in close agreement with the exact three-dimensional elasticity solutions for a wide range of lamination and geometric parameters.
Molten salt reactor designs such as the FHR must consider the safety implications of freezing, requiring the development of thermal hydraulics models that account for freezing behavior. Pipe models that couple heat transfer, fluid dynamics, and phase change physics will require verification and validation. A static freezing study is being conducted to provide an initial set of validation data with LiF-BeF{sub 2} (FLiBe). The study will also be performed with LiF-NaF-KF (FLiNaK) in order to quantify the distortion that arises from using FLiNaK as a non-toxic surrogate salt for thermal-hydraulic studies. In order to develop the numerical model, accurate inputs for the thermal hydraulic properties of the two salts of interest need to be measured near the freezing point. For this reason, a series of three experiments including Differential Scanning Calorimetry (DSC) measurements, the Density Measurement Experiment, and the Static Freezing Experiment are being conducted to determine these thermo-physical properties at temperatures near the freezing point. The Static Freezing Experiment is also designed to determine the degree of supercooling that FLiBe and FLiNaK experience. Supercooling is a temperature depression below the freezing point that the liquid experiences before freezing initiates. Increased supercooling in a liquid flowing through a pipe leads to increased dendritic freezing and an increased likelihood of pipe blockage and freeze-shut. Based on previous research on supercooling in water and non-fluoride salts, supercooling is likely affected most strongly by cooling rate and sample size (as well as salt purity). No data has been generated previously on supercooling in fluoride salts such as FLiBe and FLiNaK. The results of this set of DSC, Density Measurement and Static Freezing experiments will be used in an ongoing fluoride salt freezing study being conducted to develop pipe freezing models for molten salts to be used in safety analysis. The thermal hydraulic properties measured in this study will serve as inputs for all future models. With the supercooling behavior of both salts quantified, it can be assessed whether supercooling must be included in future thermal hydraulic models. The validated Static Freezing Experiment model can then serve as the framework for a pipe freezing model and eventually a 1D system code.
A coherent bistatic vegetation scattering model, based on a Monte Carlo simulation, is being developed to simulate polarimetric bi-static reflectometry at VHF/UHF-bands (240-270 MHz). The model is aimed to assess the value of geostationary satellite signals of opportunity to enable estimation of the Earth's biomass and root-zone soil moisture. An expression for bistatic scattering from a vegetation canopy is derived for the practical case of a ground-based/low altitude platforms with passive receivers overlooking vegetation. Using analytical wave theory in conjunction with distorted Born approximation (DBA), the transmit and receive antennas effects (i.e., polarization, orientation, height, etc.) are explicitly accounted for. Both the coherency nature of the model (joint phase and amplitude information) and the explicit account of system parameters (antenna, altitude, polarization, etc) enable one to perform various beamforming techniques to evaluate realistic deployment configurations. In this paper, several test scenarios will be presented and the results will be evaluated for feasibility for future biomass and root-zone soil moisture application using geostationary communication satellite signals of opportunity at low frequencies.
The performance of a three-dimensional continuum damage mechanics (CDM) model for intralaminar damage, coupled with a cohesive zone model for delamination, is examined for damage prediction in cross-ply and quasi-isotropic laminated plates subjected to quasi-static transverse indentation. Details of the CDM model based on the deformation gradient decomposition (DGD) method and the finite element (FE) modeling strategy are presented and discussed. Critical features of the damage model and FE modeling approach are assessed by comparing analysis predictions to experimental data for cross-ply and quasi-isotropic panels in which damage is dominated by complex interactions between matrix cracks and delaminations.
Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.
This report provides an overview on the verification and validation (V&V) of the Thermal Energy Distribution System (TEDS) model developed in the Modelica process modeling ecosystem using experimental data. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Energy Systems Laboratory (ESL) at Idaho National Laboratory (INL). The model was then used during the preconstruction phase of the experimental effort to inform experimental design (e.g., insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit: a 200kW Chromalox heater; a single-tank packed-bed thermal energy storage system filled with 0.125-inch alumina (Al2O3) beads; an ethylene-glycol-to-Therminol-66 heat exchanger; system piping; five control valves; and all associated temperature, pressure, and volumetric flow sensors. Using the Institute of Electrical and Electronics Engineers (IEEE) V&V methodologies, considered the gold standard in the engineering field, the model was verified using a combination of static analysis, spatial convergence, and regression tests. Then using dynamic time warping (DTW) initial runs to validate and tune the TEDS model versus the experiment were conducted. This tuning method was accomplished using the INL Risk Analysis Virtual ENvironment (RAVEN) software package. Tuning is required to account for physical phenomena that are less understood within the empirical heat transfer correlations. Through the commencement of this work, a systems-level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model was utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Then, initial experimental startup and operational data were used to demonstrate the validation and tuning methodology. This process demonstrates the classical two-step approach of a model informing experimental design followed by the experiment validation and tuning the model.
This paper presents an automation process for structural sizing of subsonic and supersonic aircraft concepts under static aeroelastic constraints. The automation process starts with an OpenVSP geometry and ends with a PATRAN plot of a NASTRAN solution for static aeroelastic analysis or optimization. ModelCenter is used to integrate all analysis codes with easy-to-use interfaces. Automation tools are developed to streamline the setup process and avoid user errors. Fuel is distributed by solving an optimization problem to match the center of gravity of aircraft at a specified flight condition. Fuel weights are also automatically attached to the structural model as point masses. All other weights used in FLOPS mission analysis (excluding fuselage and wing structural weights) are automatically attached to or smeared on the structural model. For any given OpenVSP geometry and FLOPS analysis data, a static aeroelastic sizing model for NASTRAN analysis can be generated in a couple of hours. The empirical fuselage and wing structural weights from FLOPS are replaced by structural panel weights from the sized finite-element model. Three supersonic and two subsonic aircraft concepts are used to demonstrate the automation process as a physics-based weight estimation tool for aircraft conceptual design.
In this work, we present a 3D Phase Field Dislocation Dynamics (PFDD) model for body-centered cubic (BCC) metals. The model formulation is extended to account for the dependence of the Peierls barrier on the line-character of the dislocation. Simulations of the expansion of a dislocation loop belonging to {110} < 111 > the slip system are presented with direct comparison to Molecular Statics (MS) simulations. The extended PFDD model is able to capture the salient features of dislocation loop expansion predicted by MS simulations. In conclusion, the model is also applied to simulate the motion of a straight screw dislocation through kink-pair motion.
Investigation of space shuttle orbiter configurations 2A(VL70-000089B) and 3(VL70-000139B) was performed at the Langley Research Center Unitary Plan Wind Tunnel (UPWT) from June 1, 1973, to June 15, 1973, for 60 test hours. The primary test objectives were to obtain stability and control characteristics for Configurations 2A and 3 and an alternate forebody used with Configuration 3. In addition, hinge moments were measured on the elevons and rudder for Configuration 2A only. The configurations were tested at Mach numbers 2.5, 3.9 and 4.6. Pitch runs were made at angles of attack from -4 to 44 deg and sideslip angles from -4 to +6. Static pressures were measured at the fuselage base for use in force data correction.
A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.
Several finite-element models are applied to the linear static, stability, and vibration analysis of laminated composite plates and shells. The study is based on linear shallow-shell theory, with the effects of shear deformation, anisotropic material behavior, and bending-extensional coupling included. Both stiffness (displacement) and mixed finite-element models are considered. Discussion is focused on the effects of shear deformation and anisotropic material behavior on the accuracy and convergence of different finite-element models. Numerical studies are presented which show the effects of increasing the order of the approximating polynomials, adding internal degrees of freedom, and using derivatives of generalized displacements as nodal parameters.
The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.
The thermal models used for analyzing dendritic web growth and calculating the thermal stress were reexamined to establish the validity limits imposed by the assumptions of the models. Also, the effects of thermal conduction through the gas phase were evaluated and found to be small. New growth designs, both static and dynamic, were generated using the modeling results. Residual stress effects in dendritic web were examined. In the laboratory, new techniques for the control of temperature distributions in three dimensions were developed. A new maximum undeformed web width of 5.8 cm was achieved. A 58% increase in growth velocity of 150 micrometers thickness was achieved with dynamic hardware. The area throughput goals for transient growth of 30 and 35 sq cm/min were exceeded.