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,513 records · Page 84

Interpretable machine learning models classify minerals via spectroscopy

Developing methods to identify mineral species confidently and rapidly from Raman spectral analysis is critical to numerous fields. Traditionally, analysis relies on pattern matching the Raman spectrum of an unknown dataset with a supporting library of well-characterized spectral data, which may prove difficult for environmental samples that are poorly crystalline or phase mixtures. Here, we developed interpretable machine learning models that can classify uranium minerals by secondary oxyanion chemistry and other physicochemical properties based solely on Raman spectra. This new ML method produces a mineral profile of physical and chemical properties for an unknown sample and can rapidly classify or identify unknown minerals from Raman data, without the need for an exact pattern match in a spectral library. Training models are validated by 1. Strong correlation of high confidence model regions with published spectroscopic assignments and 2. Correct classification of a mineral not present in training data. Training data are from the Compendium of Uranium Raman and Infrared Experimental Spectra and available crystallographic information files within the open-source Smart Spectral Matching scientific framework. Physically meaningful classifier models can rapidly identify key structural and chemical information about unknown uranium minerals and the overall methodology is broadly applicable for mineral phases.

Machine learning↗

Towards Effective Clustering Techniques for the Analysis of Electric Power Grids

Clustering is an important data analysis technique with numerous applications in the analysis of electric power grids. Standard clustering techniques are oblivious to the rich structural and dynamic information available for power grids. Therefore, by exploiting the inherent topological and electrical structure in the power grid data, we propose new methods for clustering with applications to model reduction, locational marginal pricing, phasor measurement unit (PMU or synchrophasor) placement, and power system protection. We focus our attention on model reduction for analysis based on time-series information from synchrophasor measurement devices, and spectral techniques for clustering. By comparing different clustering techniques on two instances of realistic power grids we show that the solutions are related and therefore one could leverage that relationship for a computational advantage. Thus, by contrasting different clustering techniques we make a case for exploiting structure inherent in the data with implications for several domains including power systems.

clustering, k-means, power grid, graphs↗

Development and Validation of High Performance Unshrouded Centrifugal Impeller

The feasibility of using a two-stage unshrouded impeller turbopump to replace the current three-stage reusable launch vehicle engine shrouded impeller hydrogen pump has been evaluated from the standpoint of turbopump weight reduction and overall payload improvement. These advantages are a by-product of the higher tip speeds that an unshrouded impeller can sustain. The issues associated with the effect of unshrouded impeller tip clearance on pump efficiency and head have been evaluated with one-dimensional tools and full three-dimensional rotordynamic fluid reaction forces and coefficients have been established through time dependent computational fluid dynamics (CFD) simulation of the whole 360 degree impeller with different rotor eccentricities and whirling ratios. Unlike the shrouded impeller, the unshrouded impeller forces are evaluated as the sum of the pressure forces on the blade and the pressure forces on the hub using the CFD results. The turbopump axial thrust control has been optimized by adjusting the first stage impeller backend wear ring seal diameter and diverting the second stage backend balance piston flow to the proper location. The structural integrity associated with the high tip speed has been checked by analyzing a 3D-Finite Element Model at maximum design conditions (6% higher than the design speed). This impeller was fabricated and tested in the NASA Marshall Space Flight Center water-test rig. The experimental data will be compared with the analytical predictions and presented in another paper. The experimental data provides validation data for the numerical design and analysis methodology. The validated numerical methodology can be used to help design different unshrouded impeller configurations.

Chen, Wei-Chung↗

Regional Scale Meteorological Analysis and Prediction Using GPS Occultation and EOS Data

The main objective of the research under this award is to improve regional meteorological analysis and prediction for traditionally data limited regions, particularly over the Southern Ocean and Antarctica, using the remote sensing observations from current and upcoming GPS radio occultation missions and the EOS instrument suite. The major components of this project are: 1.Develop and improve the methods for retrieving temperature, moisture, and pressure profiles from GPS radio occultation data and EOS radiometer data. 2. Develop and improve a regional scale data assimilation system (MM5 4DVAR). 3. Perform case studies involving data analysis and numerical modeling to investigate the impact of different data for regional meteorological analysis and the importance of data assimilation for regional meteorological simulation over the Antarctic region. 4. Apply the findings and improvements from the above studies to weather forecasting experiments. 5. In the third year of the award we made significant progress toward the remaining goals of the project. The work included carefully evaluating the performance of an atmospheric mesoscale model, the Polar MM5 in Antarctic applications and improving the upper boundary condition.

Bromwich, David H.↗

Lack of robustness and accuracy of many numerical schemes for phase-field simulations

In this paper, we study the stability, accuracy and convergence behavior of various numerical schemes for phase-field modeling through a simple ODE model. Both theoretical analysis and numerical experiments are carried out on this ODE model to demonstrate the limitation of most numerical schemes that have been used in practice. One main conclusion is that the first-order fully implicit scheme is the only robust algorithm for phase-field simulations while all other schemes (that have been analyzed) may have convergence issue if the time step size is not exceedingly small. More specifically, by rigorous analysis in most cases, we have the following conclusions: (i) The first-order fully implicit scheme converges to the correct steady state solution for all time step sizes. In the case of multiple solutions, one of the solution branches always converges to the correct steady state solution. (ii) The first-order convex splitting scheme, which is equivalent to the first-order fully implicit scheme with a different time scaling, always converges to the correct steady state solution but may seriously lack numerical accuracy for transient solutions. (iii) For the second-order fully implicit and convex splitting schemes, for any time step size $δt > 0$, there exists an initial condition $u_0$ , with $|u_0| > 1$, such that the numerical solution converges to the wrong steady state solution. (iv) For $|u_0|$ ≤ 1, all second-order schemes studied in this paper converge to the correct steady state solution although severe numerical oscillations occur for most of them if the time step size is not sufficiently small. (v) An unconditionally energy-stable scheme (such as the modified Crank–Nicolson scheme) is not necessarily better than a conditionally energy-stable scheme (such as the Crank–Nicolson scheme). Finally, most, if not all, of the above conclusions are expected to be true for more general Allen–Cahn and other phase-field models.

97 MATHEMATICS AND COMPUTING↗

Kinetic Alfvén wave turbulence and formation of localized structures

This work presents non-linear interaction of magnetosonic wave with kinetic Alfvén wave for intermediate β-plasma (m{sub e}/m{sub i}≪β≪1). A set of dimensionless equations have been developed for analysis by considering ponderomotive force due to pump kinetic Alfvén wave in the dynamics of magnetosonic wave. Stability analysis has been done to study modulational instability or linear growth rate. Further, numerical simulation has been carried out to study the nonlinear stage of instability and resulting power spectrum applicable to solar wind around 1 AU. Due to the nonlinearity, background density of magnetosonic wave gets modified which results in localization of kinetic Alfvén wave. From the obtained results, we observed that spectral index follows k{sup −3.0}, consistent with observation received by Cluster spacecraft for the solar wind around 1 AU. The result shows the steepening of power spectrum which may be responsible for heating and acceleration of plasma particles in solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

Radhadrishnan, Krishnan↗

Improved analysis of electron penetration and numerical procedures for space radiation shielding

Electron penetration calculational techniques are reviewed with regard to their suitability for shield analysis in future space operations. Methods based on the transmission factors of Mar are discussed and a correction term for low-energy electrons, which results in slightly conservative shield estimates, is derived. This modified Mar's method provides estimates of the dose for electrons that penetrate through shields of arbitrary elemental material with an atomic number greater than four. A complete computer algorithm is included.

Wilson, J. W.↗

Ignition and structure of a laminar diffusion flame in the field of a vortex

The distortion of flames in flows with vortical motion is examined via asymptotic analysis and numerical simulation. The model consists of a constant density, one step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a vortex. The evolution in time of the temperature and mass fraction fields is followed. Emphasis is placed on the ignition time and location as a function of vortex Reynolds number and initial temperature differences of the reacting species. The study brings out the influence of the vortex on the chemical reaction. In all phases, good agreement is observed between asymptotic analysis and the full numerical solution of the model equations.

Macaraeg, Michele G.↗

Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods

Numerical simulation of an electrodynamic system in empty space requires the implementation of open boundary conditions (BC) to terminate the solution of Maxwell’s equations on the boundaries of the computational domain. The Perfectly Matched Layer (PML) has become the method of choice for open BC with wave equations, as it is straightforward and relatively easy to implement, and offers very efficient and user-adjustable absorption rates. PMLs are most often employed with the Finite-Difference Time-Domain (FDTD) algorithm, which in its most common implementation offers second-order accuracy in space and time on Cartesian grids. Yet, simulations (including some class of electromagnetic Particle-In-Cell simulations) that require higher precision may resort to higher-order Maxwell solvers employing extended finite-difference stencils, or even to pseudo-spectral Maxwell solvers, for which a general, versatile and efficient formulation of the PML has been missing so far. Here we propose a novel “two-step” formulation of the PML that is simple, very versatile and can be used as is with any Maxwell solver. In particular, it is applicable to a large class of Maxwell solvers including the arbitrary-order Pseudo-Spectral Analytical Time-Domain (PSATD) solver, which offers arbitrarily low numerical dispersion when increasing solver order and becomes dispersion-free at infinite order. Analysis and numerical simulations demonstrate that the new formulation is as efficient as the standard PML formulation, both for the FDTD and the PSATD implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fully explicit nonlinear optics model in a particle-in-cell framework

A numerical technique which incorporates the nonlinear optics of anisotropic crystals into a particle-in-cell framework is described. The model is useful for simulating interactions between crystals, ultra-short laser pulses, intense relativistic electron bunches, plasmas, or any combination thereof. The frequency content of the incident and scattered radiation is limited only by the resolution of the spatial and temporal grid. A numerical stability analysis indicates that the Courant condition is more stringent than in the vacuum case. Numerical experiments are carried out illustrating the electro-optic effect, soliton propagation, and the generation of fields in a crystal by a relativistic electron bunch.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recent developments in thermal analysis of large space structures

A numerical procedure for analysis of shadowed space heating of sparse structures, SSQ, is discussed. The SSQ program avoids inordinate computational complexity by confining attention to a single elemental location on a structural member of interest throughout an entire orbital period, proceeding then to similar treatment of individual alternate locations. The procedure considers a spacecraft in circular orbit and assumes fixed-Earth orientation of the spacecraft. Shadow orientation and interval duration, merged shadows, and computation of solar heat flux and thermal response are addressed. The output options of the SSQ FORTRAN 5 program and its efficiency are discussed. Application of the system to the analysis of a parabolic expandable truss antenna is considered.

Oneill, R. F.↗

Demonstration of Probabilistic Sensitivity Analyses Tools on the Structural Response of a Representative Inflatable Space Structure

This work provides an initial step toward demonstrating a probabilistic numerical simulation capability to support trade studies and the development of a certification plan for inflatable space habitats. This study concentrates on interpreting the results from probabilistic analysis and numerical simulation tools to identify parameter sensitivities for a novel inflatable airlock concept, specifically the Non‐Axisymmetric Inflatable Pressure Structure (NAIPS) that was designed and tested under NASA's Minimalistic Advanced Softgoods Hatch (MASH) Program. A brief overview of the finite element model is provided along with the probabilistic sensitivity analysis approach. The sensitivity studies required a model that was numerically stable and efficient enough that hundreds of simulations could be completed in the allotted time. Therefore, the existing full model was simplified by: extracting a quarter symmetry section of the dome; focusing on a single inflation pressure; and replacing the non‐linear material stress‐strain curves with linear, isotropic materials defined by elastic moduli. Responses of interest include the sensitivity of various structural component loads to material properties, cord lengths, inflation pressure and friction. Multiple sensitivity studies were completed and three are reported here. The first study focused on utilizing wide input parameter ranges to provide an opportunity to assess numerical robustness. The next two studies narrowed the parameter ranges to enable focus on understanding uncertainty at a fixed operating condition. The completion of the sensitivity studies improved understanding of the interdependence of multiple inputs on the responses. In addition, numerical stability of the simulations over wide parameter ranges, shows the feasibility of incorporating uncertainty‐based methods in the design and certification of inflatable space habitats. With the experience and trust gained, it is anticipated that these same methods will be applied to nonlinear, orthotropic models in the future.

Lyle, Karen H.↗

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Elishakoff, Isaac↗

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗

Comparative analysis of waveguide plasmon – polariton refractometers based on excitation of surface, symmetric, and antisymmetric plasmon modes

We present the results of a numerical study of three waveguide refractometer configurations based on surface plasmon resonance: a refractometer without a buffer layer based on exciting a surface plasmon mode, and refractometers with a buffer layer based on exciting symmetric and antisymmetric plasmon modes. It is found that the best metrological characteristics are inherent in the refractometer based on the symmetric plasmon mode: the cumulative advantage in spectral sensitivity and width of the resonance dip attains two orders of magnitude as compared to the refractometer based on antisymmetric plasmon mode and one order of magnitude as compared to the refractometer without a buffer layer. For measuring local variations of the refractive index in a thin layer the advantage of using the symmetric plasmon mode is less expressed and amounts to nearly 20 and 6 times, respectively. It is shown that higher metrological characteristics of the refractometer based on symmetric plasmon mode are achieved at the expense of its compactness: the length of the sensitive segment in such a refractometer may exceed that of the refractometers without the buffer layer and those based on the antisymmetric mode by three orders of magnitude. The minimal length of the sensitive segment is possible in the refractometer based on the antisymmetric plasmon mode, which, however, possesses the worst metrological characteristics. As to the fabrication simplicity, the scheme without the buffer layer seems most attractive. This scheme occupies the intermediate position between the refractometers based on symmetric and antisymmetric plasmon modes both in metrological parameters and in the length of the sensitive segment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Surface-enhanced Raman scattering from finite arrays of gold nano-patches

We experimentally investigate the surface-enhanced Raman scattering (SERS) response of a 2D-periodic array of square gold nano-patches, functionalized by means of a conjugated, rigid thiol. We measure a Raman signal enhancement up to 200 times more intense compared to other plasmon-based nanostructures functionalized with the same molecule, and show that the enhancement is not strictly correlated to the presence of plasmonic resonances. The agreement between experimental and theoretical results reveals the importance of a full-wave analysis based on the inclusion of the actual scattering cross section of the molecule. The proposed numerical approach may serve not only as a tool to predict the enhancement of Raman signal scattered from strongly resonant nanostructure but also as an effective instrument to engineer SERS platforms that target specific molecules.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗