SEARCH · Engineering Papers

Results for “Numerical approximation & 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗

NASA NTRSDate not supplied

Quality Control Methods of Tower Data at Kennedy Space Center’s Launch Complex 39-B and the USAF Cape Canaveral 500 Foot Tower

The National Aeronautics and Space Administration (NASA) has long used meteorological data from weather towers located at Kennedy Space Center (KSC) and the United States Air Force’s Eastern Range(ER) in support of their various launch vehicles and climate studies. Some of the most valuable data is gathered at the Launch Complex 39-B (LC39-B) and the 500 foot tower located approximately two miles west of LC39-B and three miles north of the Vehicle Assembly Building (VAB). The data from LC39-B gives NASA Space Launch System (SLS) engineers valuable insight into the weather the vehicle can expect at the pad, while the 500 foot tower has a longer period of record and has been used by the Shuttle and SLS programs. However, numerous data quality (DQ) control issues have arisen when utilizing this data for analysis. In addition, instrumentation tends to deteriorate faster than normal due to the corrosive nature of the high salt content in the air. NASA’s Natural Environments branch has developed several QC databases of these towers, but these studies were for fixed periods of records. EV44 identified a need for a continually updated and QC’ed database of tower data not only to provide along term QC’ed database for vehicle and climate analyses, but also to provide the capability to investigate recent weather events such as downbursts and other high wind events. This study follows the methods developed in prior Natural Environment QC tower databases, which includes variable specific QC thresholds and checks used to generate QC flags. However, this study developed new techniques such as the development and design of the graphical user interface (GUI) for manual verification, and the creation of the final monthly QC files.

Quality Control↗

NASA NTRSDate not supplied

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗

DOE OSTI2022

The EGS Collab Project – Stimulations at Two Depths

The EGS Collab project, supported by the US Department of Energy, is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory to address challenges in implementing enhanced geothermal systems (EGS). Data and observations from the field tests are compared to simulations to understand processes and build confidence in numerical modeling of the processes. We have completed Experiment 1 (of 3), which examined hydraulic fracturing in a well-characterized underground fractured phyllite test bed at a depth of approximately 1.5 km at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. Testbed characterization included fracture mapping, borehole acoustic and optical televiewers, full waveform sonic, conductivity, resistivity, temperature, campaign p- and s-wave investigations and electrical resistance tomography. Borehole geophysical techniques including passive seismic, continuous active source seismic monitoring, electrical resistance tomography, fiber-based distributed strain, distributed temperature, and distributed acoustic monitoring, were used to carefully monitor stimulation events and flow tests. More than a dozen stimulations and nearly one year of flow tests were performed. Quality data and detailed observations were collected and analyzed during stimulation and water flow tests using ambient temperature and chilled water. We achieved adaptive control of the tests using real-time monitoring and rapid dissemination of data and near-real-time simulation. More detailed numerical simulation was performed to answer key experimental design questions, forecast fracture propagation trajectories and extents, and analyze and evaluate results. Data are freely available from the Geothermal Data Repository. Experiment 2 examines the potential for hydraulic shearing in amphibolite at a depth of about 1.25 km at SURF. This site has a different set of stress and fracture conditions than Experiment 1. The Experiment 2 testbed consists of nine subhorizontal boreholes configured in two fans of two boreholes which surround the testbed and contain grouted-in electrical resistance tomography, seismic sensors, active seismic sources and distributed fiber sensors. A “five-spot” set of test wells that extends from a custom mined alcove includes an injection well and four production/monitoring wells. The testbed was characterized geophysically and hydrologically, and three stimulations have been performed using the Step-Rate Injection Method for Fracture In-Situ Properties (SIMFIP) tool to measure strains, and a new strain quantifying tool (downhole robotic strain analysis tool -DORSA) was deployed in a monitoring hole during stimulation. Real-time data were broadcast during stimulations to allow real-time response to arising issues.

EGS Collab, Enhanced Geothermal Systems, EGS, fiel↗

NASA NTRS2001

Quasi 1-D Analysis of a Circular, Compressible, Turbulent Jet Laden with Water Droplets

Recent experimental studies indicate that presence of small amount of liquid droplets reduces the Overall Sound Pressure Level (OASPL) of a jet. Present study is aimed at numerically investigating the effect of liquid particles on the overall flow quantities of a heated, compressible round jet. The jet is assumed perfectly expanded. A quasi-1D model was developed for this purpose which uses area-averaged quantities that satisfy integral conservation equations. Special attention is given to represent the early development region since it is acoustically important. Approximate velocity and temperature profiles were assumed in this region to evaluate entrainment rate. Experimental correlations were used to obtain spreading rate of shear layer. The base flow thus obtained is then laden with water droplets at the exit of the nozzle. Mass, momentum and energy coupling between the two phases is represented using empirical relations. Droplet size and mass loading are varied to observe their effect on flow variables.

Source record↗

DOE OSTI2020

Atmospheric ice nuclei concentration measurements over a high altitude-station in the Western Ghats, India

Mixed-phase and ice clouds play a significant role in modulating the Indian summer monsoon rainfall, and modeling studies suggest that inaccurate representation of ice nucleating particles (INP) is one of the factors, which have resulted in a deficiency in the prediction of cloud properties and hydrological cycle. We do not understand the ice nucleation mechanisms to represent the ice formation in the model, and therefore more INP measurements at various supercooled temperatures are needed. This study reports for the first time real-time INP concentration at a high altitude station in the Western Ghats region of India, and these results will be used to constrain the ice nucleation parameterizations in the model. The variation of INP for any particular cloud condition can have strong microphysical responses of mixed-phase precipitation formation processes. Also, to improve the representation of heterogeneous ice nucleation, INP measurements of various types of aerosols commonly observed over the Indian region are needed. This study presents, the INP concentrations and daily concentration of non-refractory chemical composition (organics, sulfate, nitrate, ammonium, and chloride), cloud condensation nuclei (CCN), and aerosol size distribution. Spectrometer for Ice Nuclei (SPIN) was operated to measure INP concentration at the mountain site from August to December 2018. INP measurements were performed at three different temperatures (-25, -30, and - 34 °C) in the immersion freezing mode relevant for mixed-phase cloud conditions. The average INP concentrations approximately varied from 0.18 to 12.4 L-1, 0.39 to 24 L-1, 1.1 to 40.2 L-1 at -25, -30, and - 34 °C, respectively. These concentrations are within the concentrations reported in the literature from different regions across the globe. The air mass back trajectory analysis showed that continental air mass contains more INP compared to maritime air mass. The information on INP concentration and chemical characteristics of aerosol will help to improve heterogeneous ice nucleation parameterization in numerical models.

Kumar, Anil V.↗

DOE OSTI2020

A Relaxed PV Bus Model in Linear Power Flow

The emerging multimodal active energy resource units make it necessary to introduce constant voltage amplitude buses (PV bus) in the operation of modern distribution systems. To enable the analysis of PV bus using linear power flow models, this letter proposes a relaxed PV bus model for solving power flow related problems using an indirect modeling approach. The developed PV bus model can be applied to linear power flow models in both rectangular and polar coordinates. Besides, the approximation error of the relaxed PV bus model can also be controlled and quantified. Lastly, the proposed relaxed PV bus model is validated using different test systems. According to the numerical studies, the maximum error on the PV bus approximation in all case studies is only 0.61%, and the maximum error of the supporting power at the PV buses is 6.5%.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1976

Technique for atmospheric rate chemistry calculations

The possibility that predictions of atmospheric photochemistry/transport models are sensitive to uncertainties in reaction rates and other inputs stresses the need for rapid numerical integration schemes in rate photochemistry problems. Reducing the computational burden has a major merit in facilitating sensitivity studies to assess the effect of uncertainties on predicted ozone diminutions from NOx (NO + NO2) in the exhaust plume of SST engines. The paper discusses the validity of an algorithmic approach to integration of rate chemistry problems in combustion, developed by Rubel and Baronti for an approximate calculation of the production rate of the i-th chemical species involved. An analysis of two projected SST engines confirms the validity of the proposed algorithm. Because of the relative arithmetical simplicity, it may be easier to treat diffusion rate chemistry calculations using the Rubel and Baronti approximation than would be possible by other approaches.

Matloff, G. L.↗

NASA NTRS1984

Entropy condition satisfying approximations for the full potential equations of transonic flow

A class of conservative difference approximations for the steady full potential equation was presented. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. Rigorous proof indicated that these new schemes satisfied a new discrete entropy inequality, which ruled out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in the analysis is the construction of a new entropy inequality for the full potential equation itself. Results of some numerical experiments using the new schemes are presented.

Osher, S.↗

NASA NTRS1985

Entropy condition satisfying approximations for the full potential equation of transonic flow

Osher, S.↗

NASA NTRS1988

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗

DOE OSTI2021

Validity of Machine Learning in the Quantitative Analysis of Complex Scanning Near-Field Optical Microscopy Signals Using Simulated Data

Scattering-type scanning near-field optical microscope (s-SNOM) is a modern technique for subdiffractional optical imaging and spectroscopy. Over the past two decades, tremendous efforts have been devoted to modeling complex tip-sample interactions in s-SNOM, aimed at understanding the electrodynamics of materials at the nanoscale. However, due to complexities in analytical methods and the limited computation power for fully numerical simulations, compromises must be made to facilitate the modeling of tip-sample interaction, such as using quasistatic approximation or unrealistic tip geometries. Here, we apply a variety of widely utilized machine-learning methods, including k nearest neighbor and feedforward neural network etc. to study the phase-resolved spectroscopic near-field response. With only a small set of training data, which is simulated using the finite-dipole model, we demonstrate that the relation between the experimental near-field signal and sample optical constant can be one to one mapped without the need for tip modeling: for a given material with a moderate dielectric function, its complex near-field spectrum can be accurately determined within the mid-IR spectral range, and vice versa. Our preliminary study sets the stage for future exploration using real experimental data. Our method is beneficial for processing the increasing amount of data accumulated across many research groups and especially useful for user facilities such as synchrotron-based national laboratories where a large amount of data is generated on a daily basis.

36 MATERIALS SCIENCE↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

Vijayaraghavan, D.↗

NASA NTRS1972

Study of low Reynolds number nozzle flows, including radial pressure gradients

An analysis is presented of the laminar, axisymmetric flow in a nozzle, including both axial and radial variations of the pressure. The system of equations derived is believed to contain all of the terms necessary for describing the flow through a relatively sharp throat (i.e., one for which the longitudinal radius of curvature of the throat is comparable to, or less than, the transverse radius). A finite difference approximation of these equations is described, together with a computer program for finding numerical solutions. An instability was found in the starting solution; a series of attempts to eliminate this instability is described.

Rae, W. J.↗

NASA NTRS1989

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗

NASA NTRS1986

A three-dimensional viscous flow analysis for the helicopter tip vortex generation process about square and round tipped blades

The tip vortex flow field occurring in the vicinity of the tip region of a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly effect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results.

Lin, S.-J.↗

NASA NTRS1978

Analytical and numerical investigation of structural response of compliant wall materials

Theoretical analysis of an electrostatically driven wall system for a compliant wall drag reduction program is reported. The electrostatic wall system is capable of producing deflections of many orders greater than the thicknesses and at small wavelengths. An intermediate large response theory was used for structural analysis. The theoretical predictions were compared to bench test results, and good agreement between the two was obtained. The effects of aerodynamic loads and perturbation electric fields on the theoretical solutions were considered. It was shown that for very small wavelengths (approximately 2mm) the aerodynamic effects can be estimated using potential theory without loss of accuracy, and the perturbation electric fields do not affect solutions as long as the deflections are less than one percent of the wavelength. Resonance effects for this type of structure were shown to be fairly small.

Goglia, G. L.↗