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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A general-purpose computer program for studying ultrasonic beam patterns generated with acoustic lenses

A 3-D computer model of a piston radiator with lenses for focusing and defocusing is presented. To achieve high-resolution imaging, the frequency of the transmitted and received ultrasound must be as high as 10 MHz. Current ultrasonic transducers produce an extremely narrow beam at these high frequencies and thus are not appropriate for imaging schemes such as synthetic-aperture focus techniques (SAFT). Consequently, a numerical analysis program has been developed to determine field intensity patterns that are radiated from ultrasonic transducers with lenses. Lens shapes are described and the field intensities are numerically predicted and compared with experimental results.

Roberti, Dino↗
DOE OSTI2023

Effects of additively manufactured surface roughness on small diameter plenums with multiple side discharges for high temperature gas turbine blade cooling applications

Combined heat and power (CHP) applications have significant environmental and economic benefits that are consistent with the goals of the U.S. Department of Energy (DOE). One area that is currently being studied includes the potential benefits of CHP turbine operation at higher turbine inlet temperatures. Internal cooling concepts enabled by additive manufacturing (AM) are of primary interest. Here, the effectiveness of internal cooling is hindered by many factors such as velocity distribution of the cooling air to the hot surface considering impingement cooling. To simulate cooling air exiting from a series of orifices for internal cooling in an airfoil, a straight smooth wall tubing with multiple side discharging orifices is used and compared to additively manufactured tubing (Ti6Al4V Grade 23) with orifice size and spacing as well as inner and outer diameters identical to the smooth wall tubing. Similar to flow discharging form perforated pipes, the flow discharged from individual orifices along the tubes in this study is found to be nonuniformly distributed, and the horizontal (axial direction) momentum can be observed from the experimental data. Discharge velocities have been measured with two-dimensional (2D) particle imaging velocimetry (PIV) and single element hot wire anemometry. Numerical analysis has also been conducted to predict the velocity distributions along the orifices in the smooth wall and additive manufacturing (AM) tubing, which are inherited with surface roughness. Numerical and measured results in this study are compared, presented, and discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1995

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗
DOE OSTI2021

Quantitative analysis of Raman spectral parameters for carbon fibers: practical considerations and connection to mechanical properties

Although the literature on the Raman spectra of carbon fibers is vast, no consistent, robust predictive relationship between mechanical properties of carbon fibers and spectral parameters exists. This shortcoming is due to the use of numerous fitting functions to evaluate Raman spectra of carbon fibers and the inconsistencies in establishing the best fitting models in a statistically robust fashion. To address this gap, we present a comprehensive work on the Raman spectra of carbon fibers that combines a vast library of experimental data with a robust numerical analysis and a statistical evaluation of a wide range of suggested fitting models. This manuscript begins with a brief review of the commonly applied fitting models. Then, the Raman spectra of 32 commercially available polyacrylonitrile-based carbon fibers collected at excitation wavelengths 532, 633, and 785 nm are presented and the best fit for all fibers is evaluated based on several statistical criteria in conjunction with numerical calculations and physical arguments. The results suggest that high-performance fibers must be fit with at least five peaks, whereas high-tensile modulus fibers are best fit with at least six distinct peaks. Finally, we employ simultaneous fitting of the Raman spectra of specific fibers and wavelengths and demonstrate that strong correlations exist between mechanical properties and the D1 peak position and shape across the range of evaluated mechanical properties. Finally, we suggest straightforward improvements in fitting analysis procedures that can be implemented to increase coherency in the understanding of the underlying carbon fiber microstructure intuited from Raman spectroscopy.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Thermal, Fluid and Neutronic Analysis of an LEU Nuclear Thermal Propulsion Core

This paper describes the use of detailed multidisciplinary fluid/thermal/ structural/neutronic simulations to predict performance of the nuclear fuel elements of a Nuclear Thermal Propulsion rocket reactor. To achieve maximum performance, a rocket reactor's fuel must operate near thermal hydraulic, structural and neutronic limits where multidisciplinary interactions are important. Yet physical testing is expensive, time- consuming and risky. Lower-fidelity correlations (heat transfer) and simulations have always existed for design, and one role of detailed numerical analysis is to confirm correlation validity and accuracy. For complex and subtle issues, detailed numerical simulations may prove their value. The paper gives examples of both of these situations. Limitations of the methods and potential extensions will be explored.

Rocket Engine Design↗
DOE OSTI2025

Exploring orbital angular momentum and spin-orbit correlations for gluons at the Electron-Ion Collider

In our previous work [S. Bhattacharya , ], we introduced a pioneering observable aimed at experimentally detecting the orbital angular momentum (OAM) of gluons. Our focus was on the longitudinal double spin asymmetry observed in exclusive dijet production during electron-proton scattering. We demonstrated the sensitivity of the cos ϕ angular correlation between the scattered electron and proton as a probe for gluon OAM at small x and its intricate interplay with gluon helicity. This current work provides a comprehensive exposition, diving further into the aforementioned calculation with added elaboration and in-depth analysis. We reveal that, in addition to the gluon OAM, one also gains access to the spin-orbit correlation of gluons. We supplement our work with a detailed numerical analysis of our observables for the kinematics of the Electron-Ion Collider. In addition to dijet production, we also consider the recently proposed semi-inclusive diffractive deep inelastic scattering process, which potentially offers experimental advantages over dijet measurements. Finally, we investigate quark-channel contributions to these processes and find an unexpected breakdown of collinear factorization. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗
DOE OSTI2024

Solving the homogeneous Bethe-Salpeter equation with a quantum annealer

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other nonsymmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the GEVP by means of a quantum annealer, i.e., to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the simulated-annealing package and the D-Wave . The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1998

Buckling and Damage Resistance of Transversely-Loaded Composite Shells

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

Wardle, Brian L.↗
NASA NTRS1988

Mathematical and statistical analysis

The goal of the mathematical and statistical analysis component of RICIS is to research, develop, and evaluate mathematical and statistical techniques for aerospace technology applications. Specific research areas of interest include modeling, simulation, experiment design, reliability assessment, and numerical analysis.

Houston, A. Glen↗
NASA NTRS1993

Monte Carlo and experimental studies of nozzle flow in a low-power hydrogen arcjet

The flow in the nozzle of a small hydrogen arcjet intended for low-thrust propulsion is studied without arc ignition using numerical and experimental techniques. The flow conditions in the thruster indicate that low-density nonequilibrium effects will significantly influence the flow. Therefore, the numerical analysis is undertaken using a Monte Carlo approach. The experimental studies employ laser Raman scattering. Comparisons of the measured and computed results are made for total number density, rotational temperature, and for the number density of the first rotational level. The numerical results are found to be quite sensitive to the rotational relaxation rate, and a strong degree of thermal nonequilibrium is observed at the exit plane of the nozzle. Successful comparisons between experiment and analysis permit estimation of the rotational relaxation rate for hydrogen.

Boyd, Iain D.↗
DOE OSTI2014

LDRD HPC4Energy Wrapup Report - LDRD 12-ERD-074

High-performance computing and simulation has the potential to optimize production, distribution, and conversion of energy. Although a number of concepts have been discussed, a comprehensive research project to establish and quantify the effectiveness of computing and simulation at scale to core energy problems has not been conducted. We propose to perform the basic research to adapt existing high-performance computing tools and simulation approaches to two selected classes of problems common across the energy sector. The first, applying uncertainty quantification and contingency analysis techniques to energy optimization, allows us to assess the effectiveness of LLNL core competencies to problems such as grid optimization and building-system efficiency. The second, applying adaptive meshing and numerical analysis techniques to physical problems at fine scale, could allow immediate impacts in key areas such as efficient combustion and fracture and spallation. By creating an integrated project team with the necessary expertise, we can efficiently address these issues, delivering both near-term results as well as quantifying developments needed to address future energy challenges.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Quantum Control Landscapes Beyond the Dipole Approximation: Controllability, Singular Controls, and Resources

We investigate the control landscapes of closed n-level quantum systems beyond the dipole approximation by including a polarizability term in the Hamiltonian. The latter term is quadratic in the control field. Theoretical analysis of singular controls is presented, which are candidates for producing landscape traps. The results for considering the presence of singular controls are compared to their counterparts in the dipole approximation (i.e., without polarizability). A numerical analysis of the existence of traps in control landscapes for generating unitary transformations beyond the dipole approximation is made upon including the polarizability term. An extensive exploration of these control landscapes is achieved by creating many random Hamiltonians which include terms linear and quadratic in a single control field. The discovered singular controls are all found not to be local optima. This result extends a great body of recent work on typical landscapes of quantum systems where the dipole approximation is made. We further investigate the relationship between the magnitude of the polarizability and the fluence of the control resulting from optimization. It is also shown that including a polarizability term in an otherwise uncontrollable dipole coupled system removes traps from the corresponding control landscape by restoring controllability. We numerically assess the effect of a polarizability term on a known example of a particular three-level Λ-system with a second order trap in its control landscape. It is found that the addition of the polarizability removes the trap from the landscape. The general practical control implications of these simulations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Quinoa, Version 0.2

Quinoa is a set of computational tools that enables research and numerical analysis in fluid dynamics. At this time it remains a test-bed to experiment with various algorithms using fully asynchronous runtime systems. Currently, Quinoa consists of the following tools: (1) Walker, a numerical integrator for systems of stochastic differential equations in time. It is a mathematical tool to analyze and design the behavior of stochastic differential equations. It allows the estimation of arbitrary coupled statistics and probability density functions and is currently used for the design of statistical moment approximations for multiple mixing materials in variable-density turbulence. (2) Inciter, an overdecomposition-aware finite element field solver for partial differential equations using 3D unstructured grids. Inciter is used to research asynchronous mesh-based algorithms and to experiment with coupling asynchronous to bulk-synchronous parallel code. Two planned new features of Inciter, compared to the previous release (LA-CC-16-015), to be implemented in 2017, are (a) a simple Navier-Stokes solver for ideal single-material compressible gases, and (b) solution-adaptive mesh refinement (AMR), which enables dynamically concentrating compute resources to regions with interesting physics. Using the NS-AMR problem we plan to explore how to scale such high-load-imbalance simulations, representative of large production multiphysics codes, to very large problems on very large computers using an asynchronous runtime system. (3) RNGTest, a test harness to subject random number generators to stringent statistical tests enabling quantitative ranking with respect to their quality and computational cost. (4) UnitTest, a unit test harness, running hundreds of tests per second, capable of testing serial, synchronous, and asynchronous functions. (5) MeshConv, a mesh file converter that can be used to convert 3D tetrahedron meshes from and to either of the following formats: Gmsh, http://www.geuz.org/gmsh, Netgen, http://sourceforge.net/apps/mediawiki/netgen-mesher, ExodusII, http://sourceforge.net/projects/exodusii, HyperMesh, http://www.altairhyperworks.com/product/HyperMesh.

Bakosi, Jozsef↗
NASA NTRS1973

Finite deformation of elasto-plastic solids

A theoretical basis is established for analysis of finite deformation of metals. The observation that finite deformation of such elastoplastic materials may be viewed as a process rather than an event leads to derivation of a complete initial and boundary value problem distinguished by its quasilinear nature. This feature of the formulation motivates adoption of an incremental approach to numerical problem solving. Numerical solution capability is established for problems of plane stress and plane strain. The validity of the theory and numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation for which analytic solutions are available. Subsequently the analysis is employed for the investigation of necking in flat metal tensile bars. The results of this investigation provide the first full numerical solutions for tensile necking in plane stress and plane strain. In addition a basis is provided for assessment of the validity of stress-strain relations inferred from tensile test data.

Osias, J. R.↗
DOE OSTI2015

Tokamak elongation – how much is too much? Part 1. Theory

In this and the accompanying paper, the problem of the maximally achievable elongation κ in a tokamak is investigated. The work represents an extension of many earlier studies, which were often focused on determining κ limits due to (i) natural elongation in a simple applied pure vertical field or (ii) axisymmetric stability in the presence of a perfectly conducting wall. The extension investigated here includes the effect of the vertical stability feedback system which actually sets the maximum practical elongation limit in a real experiment. A basic resistive wall stability parameter, γτ w , is introduced to model the feedback system which although simple in appearance actually captures the essence of the feedback system. Elongation limits in the presence of feedback are then determined by calculating the maximum κ against n=0 resistive wall modes for fixed γτ w . The results are obtained by means of a general formulation culminating in a variational principle which is particularly amenable to numerical analysis. The principle is valid for arbitrary profiles but simplifies significantly for the Solov’ev profiles, effectively reducing the 2-D stability problem into a 1-D problem. Finally, the accompanying paper provides the numerical results and leads to a sharp answer of ‘how much elongation is too much’?

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Geologic mapping using integrated AIRSAR, AVIRIS, and TIMS data

The multi-sensor aircraft campaign called the 'Geologic Remote Sensing Field Experiment' (GRSFE), conducted during 1989 in the southwestern United States, collected multiple airborne remote sensing data sets and associated field and laboratory measurements. The GRSFE airborne data sets used in this study include the airborne Synthetic Aperture Radar (AIRSAR), the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS), and the Thermal Infrared Multispectral Scanner (TIMS). Each sensor's unique characteristics were used for this study in a combined analysis scheme for geologic mapping. AIRSAR was used to map structures and landforms, AVIRIS was used to map mineralogy, and TIMS was used to map lithology. Visual data integration using IHS transforms and combined numerical analysis using derived geophysical and geologic parameters with 'multispectral' techniques resulted in improved geologic mapping over that possible using each data set individually.

Kruse, Fred A.↗
NASA NTRS1950

The effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings

An analysis is presented of the effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings with partial-span constant-percent-chord ailerons extending inboard from the wing tip. The geometric variables considered are aspect ratio, taper ratio, aileron span, and aileron chord. The shape of the wing-torsional-stiffness curve is assumed and the twisting moment is considered to result solely from the pressure distribution caused by aileron deflection, so that the necessity of using a successive-approximation method is avoided. Because of the complexity of the equations resulting from the analysis, numerical calculations from the equations are presented in a series of figures. A computational form is provided to be used in conjunction with these figures so that calculations can be made without reference to the analysis.

Tucker, Warren A↗