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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure

The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.

Jameson, Antony↗
NASA NTRS1992

A general numerical model for wave rotor analysis

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Paxson, Daniel W.↗
NASA NTRS1974

A numerical-perturbation method for the nonlinear analysis of structural vibrations

A numerical-perturbation method is proposed for the determination of the nonlinear forced response of structural elements. Purely analytical techniques are capable of determining the response of structural elements having simple geometries and simple variations in thickness and properties, but they are not applicable to elements with complicated structure and boundaries. Numerical techniques are effective in determining the linear response of complicated structures, but they are not optimal for determining the nonlinear response of even simple elements when modal interactions take place due to the complicated nature of the response. Therefore, the optimum is a combined numerical and perturbation technique. The present technique is applied to beams with varying cross sections.

Nayfeh, A. H.↗
DOE OSTI2018

Numerical Modeling and Dynamic Analysis of a Wave-Powered Reverse-Osmosis System

A wave energy converter (WEC) system has the potential to convert the wave energy resource directly into the high-pressure flow that is needed by the desalination system to pump saltwater to the reverse-osmosis membrane and provide the required pressure level to generate freshwater. In this study, a wave-to-water numerical model was developed to investigate the potential use of a wave-powered desalination system (WPDS) for water production. The model was developed by coupling a time-domain radiation-and-diffraction method-based numerical tool (WEC-Sim) for predicting the hydrodynamic performance of WECs with a solution-diffusion model that was used to simulate the reverse-osmosis (RO) process. The objective of this research is to evaluate the WPDS dynamics and the overall efficiency of the system. To evaluate the feasibility of the WPDS, the wave-to-water numerical model was applied to simulate a desalination system that used an oscillating surge WEC device to pump seawater through the system. The hydrodynamics WEC-Sim simulation results for the oscillating surge WEC device were validated against existing experimental data. The RO simulation was verified by comparing the results to those from the Dow Chemical Company’s reverse osmosis system analysis (ROSA) model, which has been widely used to design and simulate RO systems. The wave-to-water model was then used to analyze the WPDS under a range of wave conditions and for a two-WECs-coupled RO system to evaluate the influence of pressure and flow rate fluctuation on the WPDS performance. The results show that the instantaneous energy fluctuation from waves has a significant influence on the responding hydraulic pressure and flow rate, as well as the recovery ratio and, ultimately, the water-production quality. Nevertheless, it is possible to reduce the hydraulic fluctuation for different sea states while maintaining a certain level of freshwater production, and a WEC array that produces water can be a viable, near-term solution to the nation’s water supply. A discussion on the dynamic impact of hydraulic fluctuation on the WPDS performance and potential options to reduce the fluctuation and their trade-offs is also presented.

16 TIDAL AND WAVE POWER↗
DOE OSTI2013

Influence of finite radial geometry on the growth rate of ion-channel free electron laser

The influence of finite radial geometry on the instability of a tenuous relativistic electron beam propagating in an ion-channel in a waveguide is investigated. The instability analysis is based on the linearized Vlasov-Maxwell equations for the perturbation about a self-consistent beam equilibrium. With the help of characteristic method the dispersion relation for the TE-mode is derived and analyzed through the numerical solutions. It is found that the positioning of the beam radius R{sub b} relative to the waveguide radius R{sub c}, and the ion-channel frequency can have a large influence on the maximum growth rate and corresponding wave number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Hot-spot induced mild ignition: Numerical simulation and scaling analysis

Mild ignition is a 'non-ideal' autoignition process observed at some conditions in rapid compression machine, shock tube and flow reactor experiments where flames and/or reaction fronts initiate due to non-uniformities within the test mixture. These flames/fronts can consume the entire mixture, or compression-heat the unburned gas, causing bulk ignition at times earlier than in homogeneous environments. This study numerically investigates mild ignition features using computational fluid dynamics techniques under a scenario where a small hot-spot is employed to initiate the process in a spherically-symmetric geometry. Only subsonic regimes are considered. The intent is to develop greater fundamental understandings of competing phenomena, and formulate a robust framework to quantify the perturbative behavior. Syngas blends (CO/H 2 = 80/20) are employed as the fuel, mixed with air at lean fuel loadings of phi = 0.2 and 0.5, with the initial temperatures and pressures covering 900-1120K, and 1.5-15 bar, respectively. Hot-spots ranging from 50 to 2,000 mu m in diameter are utilized with temperatures elevated from the surroundings by 1 to 90%, i.e., T hot /T 0 =1.01-1.90. Scaling arguments are derived to analyze competing processes of: (a) quenching vs. ignition, (b) extinction vs. propagation, and (c) flame consumption/forced autoignition vs. homogeneous ignition. Here, this work provides new insight into non-uniform autoignition phenomena, and discusses opportunities and challenges associated with mitigating it. In particular, scaling arguments indicate that experimental parameters could have potential to either suppress hot-spot initiated flame formation (e.g., via bath gases with high thermal diffusivity), or avert flame compression of the end gas (e.g., via lean/diluted mixtures, or bath gases with high heat capacity).

42 ENGINEERING↗
NASA NTRS1991

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗
NASA NTRS1992

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup -alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗
DOE OSTI2015

Local random potentials of high differentiability to model the Landscape

We generate random functions locally via a novel generalization of Dyson Brownian motion, such that the functions are in a desired differentiability class C{sup k}, while ensuring that the Hessian is a member of the Gaussian orthogonal ensemble (other ensembles might be chosen if desired). Potentials in such higher differentiability classes (k≥2) are required/desirable to model string theoretical landscapes, for instance to compute cosmological perturbations (e.g., k=2 for the power-spectrum) or to search for minima (e.g., suitable de Sitter vacua for our universe). Since potentials are created locally, numerical studies become feasible even if the dimension of field space is large (D∼100). In addition to the theoretical prescription, we provide some numerical examples to highlight properties of such potentials; concrete cosmological applications will be discussed in companion publications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

A mathematical and numerical framework for the analysis of compressible thermal convection in gases at very high temperatures

The relevance of non-equilibrium phenomena, nonlinear behavior, gravitational effects and fluid compressibility in a wide range of problems related to high-temperature gas-dynamics, especially in thermal, mechanical and nuclear engineering, calls for a concerted approach using the tools of the kinetic theory of gases, statistical physics, quantum mechanics, thermodynamics and mathematical modeling in synergy with advanced numerical strategies for the solution of the Navier–Stokes equations. The reason behind such a need is that in many instances of relevance in this field one witnesses a departure from canonical models and the resulting inadequacy of standard CFD approaches, especially those traditionally used to deal with thermal (buoyancy) convection problems. Starting from microscopic considerations and typical concepts of molecular dynamics, passing through the Boltzmann equation and its known solutions, we show how it is possible to remove past assumptions and elaborate an algorithm capable of targeting the broadest range of applications. Moving beyond the Boussinesq approximation, the Sutherland law and the principle of energy equipartition, the resulting method allows most of the fluid properties (density, viscosity, thermal conductivity, heat capacity and diffusivity, etc.) to be derived in a rational and natural way while keeping empirical contamination to the minimum. Special attention is deserved as well to the well-known pressure issue. With the application of the socalled multiple pressure variables concept and a projection-like numerical approach, difficulties with such a term in the momentum equation are circumvented by allowing the hydrodynamic pressure to decouple from its thermodynamic counterpart. The final result is a flexible and modular framework that on the one hand is able to account for all the molecule (translational, rotational and vibrational) degrees of freedom and their effective excitation, and on the other hand can guarantee adequate interplay between molecular and macroscopic-level entities and processes. Performances are demonstrated by computing some incompressible and compressible benchmark test cases for thermal (gravitational) convection, which are then extended to the high-temperature regime taking advantage of the newly developed features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Lepton number violation as a key to low-scale leptogenesis

We explore the possibility of having a successful leptogenesis through oscillations between new sterile fermion states added to the Standard Model field content in a well motivated framework, naturally giving rise to the required mass splitting between the sterile states through a small total lepton number violation. We propose a framework with only two sterile states forming a pseudo-Dirac state, in which their mass difference as well as the smallness of the neutrino masses are due to two sources of lepton number violation with ΔL=2, corresponding to an Inverse Seesaw framework extended by a Linear Seesaw mass term. We also explore the pure Inverse Seesaw mechanism in its minimal version, requiring at least four new sterile states in order to comply with neutrino data. Our analytical and numerical studies reveal that one can have a successful leptogenesis at the temperature of the electroweak scale through oscillations between the two sterile states with a “natural” origin of the strong degeneracy in their mass spectrum. We also revisit the analytical expression of the baryon asymmetry of the Universe in the weak washout regime of this framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Digital Speckle Photography of Subpixel Displacements of Speckle Structures Based on Analysis of Their Spatial Spectra

We have investigated the capabilities of the method of digital speckle interferometry for determining subpixel displacements of a speckle structure formed by a displaceable or deformable object with a scattering surface. An analysis of spatial spectra of speckle structures makes it possible to perform measurements with a subpixel accuracy and to extend the lower boundary of the range of measurements of displacements of speckle structures to the range of subpixel values. The method is realized on the basis of digital recording of the images of undisplaced and displaced speckle structures, their spatial frequency analysis using numerically specified constant phase shifts, and correlation analysis of spatial spectra of speckle structures. Transformation into the frequency range makes it possible to obtain quantities to be measured with a subpixel accuracy from the shift of the interference-pattern minimum in the diffraction halo by introducing an additional phase shift into the complex spatial spectrum of the speckle structure or from the slope of the linear plot of the function of accumulated phase difference in the field of the complex spatial spectrum of the displaced speckle structure. The capabilities of the method have been investigated in natural experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

A study of the extended-range forecasting problem blocking

Wavenumber frequency spectral analysis of a 90 day winter (Jan. 15 - April 14) wind field simulated by a climate experiment of the GLAS atmospheric circulation model is made using the space time Fourier analysis which is modified with Tukey's numerical spectral analysis. Computations are also made to examine how the model wave disturbances in the wavenumber frequency domain are maintained by nonlinear interactions. Results are compared with observation. It is found that equatorial easterlies do not show up in this climate experiment at 200 mb. The zonal kinetic energy and momentum transport of stationary waves are too small in the model's Northern Hemisphere. The wavenumber and frequency spectra of the model are generally in good agreement with observation. However, some distinct features of the model's spectra are revealed. The wavenumber spectra of kinetic energy show that the eastward moving waves of low wavenumbers have stronger zonal motion while the eastward moving waves of intermediate wavenumbers have larger meridional motion compared with observation. Furthermore, the eastward moving waves show a band of large spectral value in the medium frequency regime.

Chen, T. C.↗
NASA NTRS2018

Achieving Improved Reliability with Failure Analysis

Reliability is the ability of a product to properly function, within specified performance limits, for a specified period of time, under the life cycle application conditions. Failure analysis is a vital tool in the effort to ensure reliability of electronic products and systems throughout their product lifecycle. Today, organizations involved in activities within the electronics supply chain are facing new challenges, not just from complex assembly styles, harsher lifecycle environments, and sophisticated supply chains, but also from customers who are demanding a quicker turn-around. Unfortunately, root cause failure analysis is often performed incompletely, leading to a poor understanding of failure mechanisms and causes and, customer dissatisfaction due to recurring failures. The PDC starts with an introduction to reliability concepts, physics of failure and an overview of failure mechanisms that affect PCBs, PCBAs and components. The PDC then dives into root cause hypothesizing techniques (Pareto, FMEA, fishbone, FTA), non-destructive and destructive analysis and, materials characterization will be discussed. Numerous failure analysis case studies will be used to illustrate the techniques and analysis principles to arrive at the root cause(s) of field failures on printed circuit boards, active components, and assemblies. What Will You Learn: Topics include: Overview of Reliability Concepts Failure mechanisms of electronic products Root cause analysis Failure analysis techniques -Non-destructive techniques (optical, CSAM etc.) -Destructive analysis (DPA, Decap, FIB etc.) -Materials characterization (XRF, EDS, TMA/DSC etc.) Who Will Benefit: Reliability engineers, failure analysis engineers, engineering managers, design engineers, component engineers, quality assurance functions and, personnel involved with reliability activities within their company.

non-destructive techniques↗
NASA NTRS2019

Achieving Improved Reliability with Failure Analysis

Reliability is the ability of a product to properly function, within specified performance limits, for a specified period of time, under the life cycle application conditions. Failure analysis is a vital tool in the effort to ensure reliability of electronic products and systems throughout their product lifecycle. Today, organizations involved in activities within the electronics supply chain are facing new challenges, not just from complex assembly styles, harsher lifecycle environments, and sophisticated supply chains, but also from customers who are demanding a quicker turn-around. Unfortunately, root cause failure analysis is often performed incompletely, leading to a poor understanding of failure mechanisms and causes and, customer dissatisfaction due to recurring failures. The PDC (Professional Development Course) starts with an introduction to reliability concepts, physics of failure and an overview of failure mechanisms that affect PCBs (Printed Circuit Boards), PCBAs (Printed Circuit Board Assembly) and components. The PDC then dives into root cause hypothesizing techniques (Pareto, FMEA (Failure Modes and Effects Analysis), fishbone (Cause-And-Effect Diagram), FTA (Fault Tree Analysis)), non-destructive and destructive analysis and, materials characterization will be discussed. Numerous failure analysis case studies will be used to illustrate the techniques and analysis principles to arrive at the root cause(s) of field failures on printed circuit boards, active components, and assemblies. What Attendees will Learn: Topics include: Overview of Reliability Concepts Failure mechanisms of electronic products Root cause analysis Failure analysis techniques -Non-destructive techniques (optical, CSAM (Confocal Scanning Electron Microscopy) etc.) -Destructive analysis (DPA (Destructive Physical Analysis), Decap (Decapsulation), FIB (Focused Ion Beam) etc.) -Materials characterization (XRF (X-Ray Fluorescence) , EDS (Error Detection Sequential), TMA/DSC (Thermal Mechanical Analysis/Differential Scanning Calorimetry) etc.)

PCB quality↗
NASA NTRS2020

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗
DOE OSTI2015

Effects of thermal fluctuations on thermal inflation

The mechanism of thermal inflation, a relatively short period of accelerated expansion after primordial inflation, is a desirable ingredient for a certain class of particle physics models if they are not to be in contention with the cosmology of the early Universe. Though thermal inflation is most simply described in terms of a thermal effective potential, a thermal environment also gives rise to thermal fluctuations that must be taken into account. We numerically study the effects of these thermal fluctuations using lattice simulations. We conclude that though they do not ruin the thermal inflation scenario, the phase transition at the end of thermal inflation proceeds through phase mixing and is therefore not accompanied by the formations of bubbles nor appreciable amplitude of gravitational waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Investigation of the electron trajectories and gain regimes of the whistler pumped free-electron laser

A free-electron laser (FEL) scheme, which employs the whistler wave as a slow electromagnetic wave wiggler, was studied theoretically. Subjected to the transverse fields of whistler wave wiggler, the beam electrons are the source of the energy needed to produce electromagnetic radiation. The strength and the period of the wiggler field depend on the parameters of the magnetoplasma medium. This configuration has a higher tunability by controlling the plasma density, on top of the {gamma}-tunability of the conventional FELs. The theory of linear gain and electron trajectories was presented and four groups (I, II, III, and IV) of electron orbits were found in the presence of an axial guide magnetic field. Using perturbation analysis, it is found that these groups of orbits were stable except small regions of group I and IV orbits. The function {Phi} which determines the rate of change of axial velocity with beam energy was also derived. In the case in which {Phi}<0 represents a negative-mass regime in which the axial velocity accelerates as the electrons lose energy. Numerical solutions showed that by increasing the cyclotron frequency, the gain for group I and III orbits increased, while a gain decrement was obtained for group II and IV orbits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗