SEARCH · Engineering Papers
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.
Direct numerical simulations of a temporally developing, droplet laden mixing layer undergoing transition to mixing turbulence are conducted. The formulation includes complete two-way couplings of mass, momentum, and energy.
Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.
Numerical simulation of global atmospheric circulation for predicting performance of balloons in Eole experiment
Numerical simulation of compressible viscous flow fields is performed for a transonic transport configuration. A single structured grid system is constructed using analytical transformations such as conformal mapping, shearing/twisting/rotating/clustering/stretching transformations. The Reynolds-averaged, thin-layer Navier-Stokes equations are solved on a supercomputer, FACOM VP-400, using the LU-ADI factorization method.
The resolution requirements for direct numerical simulations of transition to turbulence are investigated. A reliable resolution criterion is determined from the results of several detailed simulations of channel and boundary-layer transition.
A 1-2/2 dimensional hybrid numerical simulation code is used to study the heating of cold H(+) ions and heavy ions by electromagnetic ion-cyclotron waves (ICWs) in the ring current region of the equatorial magnetosphere. Consideration is given to a plasma consisting of electrons, hot H(+) ions, and cold heavy ions in which the ICWs are driven by the temperature anisotropy of the hot protons. For large-amplitude ICWs, it is found that the cold H(+) ions are preferentially heated over the heavy ions although the cold H(+) ions are heated by a three-step process.
Direct numerical simulations of a compressible wake are performed to determine the sensitivity of large-scale motions to initial phase lags between the most unstable two-dimensional wave and a pair of oblique waves, and to determine the effect of Mach number on the development of three-dimensional motions. The vorticity dynamics of the flow are studied, and it is determined that, depending upon the initial phase lag, the interaction of the streamwise vorticity with the spanwise vortices may result in the formation of closed vortex loops, or in the removal of streamwise vorticity from the saddle region between the spanwise vortices. Due to lower growth rates experienced at high Mach numbers, the roll-up process is delayed, and the resultant nonlinear motions are much more diffuse.
The resolution requirements for direct numerical simulations of transition to turbulence are investigated. A reliable resolution criterion is determined from the results of several detailed simulations of channel and boundary-layer transition.
Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.
Direct numerical simulation (DNS) is used to investigate a time-dependent turbulent wake evolving in a stably stratified background. A large initial Froude number is chosen to allow the wake to become fully turbulent and axisymmetric before stratification affects the spreading rate of the mean defect. The uncertainty introduced by the finite sample size associated with gathering statistics from a simulation of a time-dependent flow is reduced, compared to earlier simulations of this flow. The DNS reveals the buoyancy-induced changes to the turbulence structure, as well as to the mean-defect history and the terms in the mean-momentum and turbulence-kinetic-energy budgets, that characterize the various states of this flow - namely the three-dimensional (essentially unstratified), non-equilibrium (or 'wake-collapse') and quasi-two-dimensional (or 'two-component') regimes observed elsewhere for wakes embedded in both weakly and strongly stratified backgrounds. The wake-collapse regime is not accompanied by transfer (or 'reconversion') of the potential energy of the turbulence to the kinetic energy of the turbulence, implying that this is not an essential feature of stratified-wake dynamics. The dependence upon Reynolds number of the duration of the wake-collapse period is demonstrated, and the effect of the details of the initial/near-field conditions of the wake on its subsequent development is examined.
In the present study numerical simulations of epitaxial growth of gallium nitride in Metal Organic Vapor Phase Epitaxy reactor AIX-200/4RF-S is addressed. Epitaxial growth means crystal growth that progresses while inheriting the laminar structure and the orientation of substrate crystals. One of the technological problems is to obtain homogeneous growth rate over the main deposit area. Since there are many agents influencing reaction on crystal area such as temperature, pressure, gas flow or reactor geometry, it is difficult to design optimal process. According to the fact that it's impossible to determine experimentally the exact distribution of heat and mass transfer inside the reactor during crystal growth, modeling is the only solution to understand the process precisely. Numerical simulations allow to understand the epitaxial process by calculation of heat and mass transfer distribution during growth of gallium nitride. Including chemical reactions in numerical model allows to calculate the growth rate of the substrate and estimate the optimal process conditions for obtaining the most homogeneous product.
The tip clearance flows of transonic compressor rotors have a significant impact on rotor and stage performance. Although numerical simulations of these flows are quite sophisticated, they are seldom verified through rigorous comparisons of numerical and measured data because, in high-speed machines, measurements acquired in sufficient detail to be useful are rare. Researchers at the NASA Glenn Research Center at Lewis Field compared measured tip clearance flow details (e.g., trajectory and radial extent) of the NASA Rotor 35 with results obtained from a numerical simulation. Previous investigations had focused on capturing the detailed development of the jetlike flow leaking through the clearance gap between the rotating blade tip and the stationary compressor shroud. However, we discovered that the simulation accuracy depends primarily on capturing the detailed development of a wall-bounded shear layer formed by the relative motion between the leakage jet and the shroud.
A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.
The magneto-static behaviour of soft magnetic composites (SMCs) is investigated using tomography based direct numerical simulation. The microgeometry crucially affects the magnetic properties of the composite since a geometry dependent demagnetizing field is established inside the composite, which lowers the magnetic permeability. We determine the magnetic field information inside the SMC using direct numerical simulation of the magnetic field based on high resolution micro-computed tomography data of the SMC's microstructure as well as artificially generated data made of statistically homogeneous systems of identical fully penetrable spheres and prolate spheroids. Quasi-static electromagnetic behaviour and linear material response are assumed. The 3D magnetostatic Maxwell equations are solved using Whitney finite elements. Simulations show that clustering and percolation behaviour determine the demagnetizing factor of SMCs rather than the particle shape. The demagnetizing factor correlates with the slope of a 2-point probability function at its origin, which is related to the specific surface area of the SMC. Comparison with experimental results indicates that the relatively low permeability of SMCs cannot be explained by demagnetizing effects alone and suggests that the permeability of SMC particles has to be orders of magnitude smaller than the bulk permeability of the particle material.
An evolutionary approach proposed in [1, 2] combining the achievements of traditional macroscopic theory of solid mechanics and basic ideas of nonlinear dynamics is applied in a numerical simulation of present-day tectonic plates motion and seismic process in Central Asia. Relative values of strength parameters of rigid blocks with respect to the soft zones were characterized by the δ parameter that was varied in the numerical experiments within δ = 1.1–1.8 for different groups of the zonal-block divisibility. In general, the numerical simulations of tectonic block motion and accompanying seismic process in the model geomedium indicate that the numerical solutions of the solid mechanics equations characterize its deformation as a typical behavior of a nonlinear dynamic system under conditions of self-organized criticality.
Here, this paper reports the systematic parameter study of a tristable nonlinear electromagnetic energy harvester for ambient low-frequency vibration. Numerical simulations and experimental investigations are performed on the harvester which consists of a cantilever beam, a tip coil, two tip magnets and two external side magnets. The external side magnets are deployed symmetrically along a concave surface parallel to the trajectory of the cantilever tip with a controllable distance so that the magnetic orientation of the tip magnets are matched with that of the side magnets. Therefore, instead of the ternary position parameters (d, h, α), a binary parameters pair (d 0 , d) is used to characterize the position of the side magnets and the performance of the energy harvester. The magnetic force and magnetic field on the cantilever tip therefore depend on the relative distance in the tip displacement direction between the tip magnets and side magnets, but is independent of the position of the side magnets on the concave surface. The magnetic force (field)-distance relationship is measured experimentally and curve fitted to obtain explicit expressions, in order to characterize the magnetic force (field) when the side magnets are placed at varied positions along the concave surface. Numerical simulation is, then, performed to predict the electromagnetic voltage output and the bandwidth of the energy harvester. The simulation results coincided with the measured data. Significant broadband response is obtained experimentally and the maximum RMS power output is 40.2 mW at 0.45 g of excitation. Furthermore, the proposed structure showcasing the matched magnetic orientation is characterized by the binary parameters pair (d 0 , d) and the systematic parametric approach could contribute to the design and study of nonlinear broadband energy harvesters.
The Probability of Detection (PoD) curve method has emerged as an important tool for the assessment of the performance of NDE techniques, a topic of particular interest to the nuclear industry where inspection qualification is very important. The conventional experimental means of generating PoD curves though, can be expensive, requiring large data sets (covering defects and test conditions), and equipment and operator time. Several methods of achieving faster estimates for PoD curves using physics-based modelling have been developed to address this problem. Numerical modelling techniques are also attractive, especially given the ever-increasing computational power available to scientists today. Here we develop procedures for obtaining PoD curves, assisted by numerical simulation and based on Bayesian statistics. Numerical simulations are performed using Finite Element analysis for factors that are assumed to be independent, random and normally distributed. PoD curves so generated are compared with experiments on austenitic stainless steel (SS) plates with artificially created notches. We examine issues affecting the PoD curve generation process including codes, standards, distribution of defect parameters and the choice of the noise threshold. We also study the assumption of normal distribution for signal response parameters and consider strategies for dealing with data that may be more complex or sparse to justify this. These topics are addressed and illustrated through the example case of generation of PoD curves for pulse-echo ultrasonic inspection of vertical surface-breaking cracks in SS plates.
The knowledge of the radionuclide content of radioactive waste is of utmost importance for safety and waste management reasons. Numerical simulations are used at EDF-CIDEN to anticipate the dismantling and the radioactive waste management. The activation scheme by neutron flux developed at EDF-CIDEN comprises four steps. First of all, a 3 dimensional description of the reactor, the chemical composition of each structure, and the neutron sources at nominal power rating conditions are used as the input data to simulate the neutron propagation and to generate a 3 dimensional multi energy group neutron flux distribution in nominal power rating conditions. Secondly, this output data is used with the chemical compositions (including impurities) and the history of the operating conditions to calculate the radiological inventory for each structure. Thirdly, according to the radioactive inventory of each component or sub-component, a waste classification can be made; the classification criteria are based on the levels of radioactivity and radiotoxicity of 143 radionuclides. Finally, comparisons between the calculated radioactive inventories and the results of radio-chemical analyses can be made. These comparisons are made with regards to the activation of 'standard' chemical elements (meaning radionuclides being perfectly modeled with no uncertainty about the chemical concentrations). The results, linked with 'standard' chemical elements, can be used to validate the calculation scheme, including the neutron propagation, the neutron activation and all the other associated simplified assumptions (geometry, history of irradiation, etc). For steel, activation of 'Fe' is retained, and for stainless steel, activation of 'Fe' 'Ni' are retained. Overall, the feedback from Chooz A (Pressurized Water Reactor) shows that numerical simulations produce a slight overestimation of 'Fe-55' and 'Ni-63'. The results linked with minor chemical elements or impurities allowed us to validate the use of the average measured compositions. (authors)