Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Beam-, Target- and Double- Spin Asymmetry Measurements in Deeply Virtual π 0 , π + and π - Meson Production

This thesis presents measurements of spin asymmetries in single pion electroproduction from the EG1-DVCS experiment carried out at the Thomas Jefferson National Accelerator Facility in 2009 with use of the CEBAF Large Acceptance Spectrometer. The experiment made use of a 6 GeV longitudinally polarised electron beam and (dynamically) longitudinally polarised targets of 14 NH 3 and 14 ND 3 (giving access to protons and deuterons), providing the means to measure three spin asymmetries: beam-, target-, and double-spin. Deeply Virtual Meson Production for the π 0 , π + and π - channels were studied for the kinematic range of Q 2 > 1 GeV 2 and W > 2 GeV. The collection of measurements includes the analysis of numerous reconstruction topologies, as well as making use of different target materials. These are the world's first charged pion spin asymmetry measurements in this kinematic regime. The sPlot maximum likelihood fit based background subtraction technique was implemented in the data analysis, showing strong agreement with a traditional and commonly used cuts-based approach. This novel technique was thoroughly tested to develop implementation and execution. The results obtained will provide a reference point for future analyses at the upgraded experimental facility and motivate the use of the sPlot technique. These measurements will provide insight into the structure of nucleons through their interpretation in relation to Generalised Parton Distributions which describe substructure correlations between longitudinal momentum and transverse position of partons. In particular, here the spin asymmetries dependence on t , the squared four-momentum transfer to the target nucleon, carries information on the spatial parton distributions in the transverse plane.

Murdoch, Gavin↗
NASA NTRS1984

Free-Vibration Analysis of Structures

Unified numerical procedure for free-vibration analysis of structures developed and incorporated into EIGSOL computer program. Dynamic response analysis of primary importance in design of wide range of practical structures such as space-craft, buildings, and rotating machineries. Procedure determines natural frequencies and associated modes in structural design.

Gupta, K. K.↗
NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗
NASA NTRS2001

Advanced Methods for Aircraft Engine Thrust and Noise Benefits: Nozzle-Inlet Flow Analysis

Three connected sub-projects were conducted under reported project. Partially, these sub-projects are directed to solving the problems conducted by the HU/FM&AL under two other NASA grants. The fundamental idea uniting these projects is to use untraditional 3D corrugated nozzle designs and additional methods for exhaust jet noise reduction without essential thrust lost and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves forming in propulsion systems; this mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) to use porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of exhaust hot jet and pressure compensation for non-design conditions (so-called continuous ejector with small mass flow rate; and (3) to propose and analyze new effective methods fuel injection into flow stream in air-breathing engines. Note that all these problems were formulated based on detailed descriptions of the main experimental facts observed at NASA Glenn Research Center. Basically, the HU/FM&AL Team has been involved in joint research with the purpose of finding theoretical explanations for experimental facts and the creation of the accurate numerical simulation technique and prediction theory for solutions for current problems in propulsion systems solved by NASA and Navy agencies. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analysis for advanced aircraft and rocket engines. The F&AL Team uses analytical methods, numerical simulations, and possible experimental tests at the Hampton University campus. We will present some management activity and theoretical numerical simulation results obtained by the FM&AL Team in the reporting period in accordance with the schedule of the work.

Morgan, Morris H.↗
DOE OSTI2017

An improved rotated staggered-grid finite-difference method with fourth-order temporal accuracy for elastic-wave modeling in anisotropic media

The rotated staggered-grid (RSG) finite-difference method is a powerful tool for elastic-wave modeling in 2D anisotropic media where the symmetry axes of anisotropy are not aligned with the coordinate axes. We develop an improved RSG scheme with fourth-order temporal accuracy to reduce the numerical dispersion associated with prolonged wave propagation or a large temporal step size. The high-order temporal accuracy is achieved by including high-order temporal derivatives, which can be converted to high-order spatial derivatives to reduce computational cost. Dispersion analysis and numerical tests show that our method exhibits very low temporal dispersion even with a large temporal step size for elastic-wave modeling in complex anisotropic media. Using the same temporal step size, our method is more accurate than the conventional RSG scheme. In conclusion, our improved RSG scheme is therefore suitable for prolonged modeling of elastic-wave propagation in 2D anisotropic media.

58 GEOSCIENCES↗
DOE OSTI2018

The Demand-Side Grid (dsgrid) Model Documentation

Electrical load is the backdrop for all power systems analysis. As a coequal partner in the supply-demand balance that must be maintained on electrical grids at all times and for all timescales, load is a major source of variability and uncertainty in grid operations. Despite this, utility load forecasting is typically done in a top-down manner that lacks the granularity in time, geography, end use, and technology that is needed to explore the potential impact of technological shifts; and many studies of future grid systems employ relatively simple scaling of historical loads, placing more emphasis on supply-side resources such as generation and transmission that are fewer in number and are better understood in terms of cost and performance. The primary purpose of the demand-side grid model (dsgrid) is to fill these gaps by creating comprehensive load data sets at a sufficient temporal, geographic, sectoral, and end-use resolution to enable detailed analyses of current patterns and future projections of end-use load. Furthermore, the dsgrid platform uses a bottom-up methodology that allows highly resolved analysis of 'what-if' scenarios. dsgrid leverages detailed sectoral energy models to provide hourly time series of load by subsector, end use, and county covering a full year. Although dsgrid currently emphasizes electricity load data, its component sector models for residential buildings, commercial buildings, and industry provide information on other fuel use, including natural gas. The data sets can thus be leveraged to support analysis of numerous demand-side technology-driven changes, such as energy efficiency, electrification, and operational flexibility (i.e., demand response). The electricity use data are time-synchronized with solar and wind data sets so as to be suitable for use in power systems analysis. This report documents dsgrid and its initial data set, which covers electricity demand in the contiguous United States for the historical year of 2012. We start with a historical year to enable model calibration and validation, as well as the development of a mathematical description of our model residuals that can be applied to future-year 'what-if' scenarios to remove known errors between our bottom-up modeling and historical load shapes. This work is part of the Electrification Futures Study (EFS), for which the dsgrid team will be developing future load snapshots that describe projected year 2050 under baseline and a range of electrification assumptions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Plasma Facing Components with Capillary Porous System and Liquid Metal Coolant Flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [Khodak and Maingi (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile, as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a Fusion Nuclear Science Facility [Kessel et al. (2019)]: 10T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

High Heat Flux↗
NASA NTRS1996

An aeroelastic analysis of helicopter rotor blades incorporating piezoelectric fiber composite twist actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consist of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for numerical integration. The twist actuation responses for three conceptual full-scale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗
NASA NTRS2018

Probabilistic Sensitivity Analyses Incorporating Surrogate Models for the Structural Response of a Representative Inflatable Space Structure

This study concentrates on results from probabilistic analysis and numerical simulation tools to identify parameter sensitivities and assess surrogate model suitability for a novel inflatable airlock concept, specifically the Non-Axisymmetric Inflatable Pressure Structure (NAIPS) under NASA's Minimalistic Advanced Softgoods Hatch (MASH) Program. The current studies extended the demonstration of probabilistic analysis tools to identify parameter sensitivities by: incorporating nonlinear material load-strain property curves and orthotropic representation of the fabric; and exploring the use of surrogate models to facilitate rapid evaluations as new information becomes available, as well as enable optimizations. The airlock model has been simplified to focus on the end dome section. Brief overviews of the finite element model and associated surrogate models are provided along with the probabilistic sensitivity analysis approaches. Responses of interest include the sensitivities of various loads in the structural elements of the NAIPS to material properties, cord lengths, inflation pressure and friction between softgoods components. The use of surrogate models was explored to facilitate rapid sensitivity evaluations as parameter ranges change. The completion of the sensitivity studies improved understanding of the dependence of load responses to several uncertain parameters and confidence in the ability to use surrogate models to represent the finite element simulations.

Lyle, Karen H.↗
DOE OSTI2016

On Improving Analytical Models of Cosmic Reionization for Matching Numerical Simulations

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large scale statistical properties. These mock catalogs are particularly useful for CMB polarization and 21cm experiments, where large volumes are required to simulate the observed signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON IMPROVING ANALYTICAL MODELS OF COSMIC REIONIZATION FOR MATCHING NUMERICAL SIMULATION

The methods for studying the epoch of cosmic reionization vary from full radiative transfer simulations to purely analytical models. While numerical approaches are computationally expensive and are not suitable for generating many mock catalogs, analytical methods are based on assumptions and approximations. We explore the interconnection between both methods. First, we ask how the analytical framework of excursion set formalism can be used for statistical analysis of numerical simulations and visual representation of the morphology of ionization fronts. Second, we explore the methods of training the analytical model on a given numerical simulation. We present a new code which emerged from this study. Its main application is to match the analytical model with a numerical simulation. Then, it allows one to generate mock reionization catalogs with volumes exceeding the original simulation quickly and computationally inexpensively, meanwhile reproducing large-scale statistical properties. These mock catalogs are particularly useful for cosmic microwave background polarization and 21 cm experiments, where large volumes are required to simulate the observed signal.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Computational analyses of crystal growth

Two important aspects of Hg/Cd/Te crystal growth processes are discussed. First, the thermal field and second, the fluid movement in the melt zone. The thermal analysis includes numerical calculation of axisymmetric heat conduction within the sample. It also includes a three-dimensional radiation model to calculate the radiative heat exchange between the furnace and the crystal as determined by the complex geometry of the furnace and the adiabatic shield. The thermal analysis also includes a crystal conductivity which is dependent on temperature and composition. To tackle the fluid flow aspect of the problem, an attempt was made to use a newly developed incompressible flow code based on the slight compressibility, and hence the finite sound speed, of all real fluids.

Dakhoul, Youssef M.↗
NASA NTRS1989

Buoyancy effects on the vapor condensation rate on a horizontal liquid surface

The results are presented of a numerical study of the effects of buoyancy on the direct condensation of saturated or nearly saturated vapor on a horizontal liquid surface in a cylindrical tank. The liquid motion beneath the liquid-vapor interface is induced by an axisymmetric laminar jet of subcooled liquid. Analysis and numerical results show that the dominant parameter which determines the influence of buoyancy on the condensation rate is the Richardson number. However, the effect of buoyancy on the condensation rate cannot be quantified in terms of the Richardson number alone. The critical value of the Richardson number below which the condensation rate is not significantly reduced depends on the Reynolds number as well as the Prandtl number.

Hasan, Mohammad M.↗
NASA NTRS1990

Buoyancy effects on the vapor condensation rate on a horizontal liquid surface

The results are presented of a numerical study of the effects of buoyancy on the direct condensation of saturated or nearly saturated vapor on a horizontal liquid surface in a cylindrical tank. The liquid motion beneath the liquid-vapor interface is induced by an axisymmetric laminar jet of subcooled liquid. Analysis and numerical results show that the dominant parameter which determines the influence of buoyancy on the condensation rate is the Richardson number. However, the effect of buoyancy on the condensation rate cannot be quantified in terms of the Richardson number alone. The critical value of the Richardson number below which the condensation rate is not significantly reduced depends on the Reynolds number as well as the Prandtl number.

Hasan, Mohammad M.↗
NASA NTRS1972

Theoretical study of corrugated plates: Shear stiffness of a trapezoidally corrugated plate with discrete attachments to a rigid flange at the ends of the corrugations

Analysis and numerical results are presented for the elastic shear stiffness of a corrugated shear web with a certain type of discrete attachments at the ends of the trough lines of the corrugations, namely point attachments to a rigid flange which interferes with the deformations of the end cross sections by preventing downward movement but permitting upward (lifting off) movement. The analysis is based on certain assumed modes of deformation of the cross sections in conjunction with the method of minimum total potential energy and the calculus of variations in order to obtain equations for the manner in which the assumed modes of deformation vary along the length of the corrugation. The numerical results are restricted to the case of equal-width crests and troughs but otherwise apply to a wide variety of geometries. They are in the form of graphs which give the overall shear stiffness as a fraction of the overall shear stiffness that could be obtained by having continuous attachment at the ends of the corrugations.

Hsiao, C.↗
NASA NTRS1996

Aeroelastic Analysis of Helicopter Rotor Blades Incorporating Anisotropic Piezoelectric Twist Actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for dynamics simulation using numerical integration. The twist actuation responses for three conceptual fullscale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗