Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 757 records · Page 42

Crush Test Abuse Stand

The purpose of this system is to simulate an internal short on battery cells by causing deformation (a crushing force) in a cell without penetration. This is performed by activating a hydraulic cylinder on one side of a blast wall with a hydraulic pump located on the other. The operator can control the rate of the crush by monitoring a local pressure gauge connected to the hydraulic cylinder or a load cell digital display located at the hydraulic pump control area. The internal short simulated would be considered a worst-case scenario of a manufacturer fs defect. This is a catastrophic failure of a cell and could be a very destructive event. Fully charged cells are to have an internal short simulated at the center of the length of the cell (away from terminals). The crush can be performed with a .- to 1-in. (.0.6- to 2.5-cm) rod placed crossways to the cell axis, causing deformation of the cell without penetration. The OCV (open-circuit voltage) and temperature of the cells, as well as the pressure and crushing force, are recorded during the operation. Occurrence of an internal short accompanied by any visible physical changes such as venting, fires, or explosions is reported. Typical analytical data examined after the test would be plots of voltage, temperature, and pressure or force versus time. The rate of crushing force can be increased or decreased based on how fast the operator pumps the hydraulic pump. The size of cylinder used to compress the battery cell can be easily changed by adding larger or smaller fittings onto the end of the hydraulic cylinder based on the battery/cell size being tested. The cell is crushed remotely and videotaped, allowing the operator to closely monitor the situation from a safe distance.

Collins, Jacob↗

Numerical study of axial turbulent flow over long cylinders

The effects of transverse curvature are investigated by means of direct numerical simulations of turbulent axial flow over cylinders. Two cases of Reynolds number of about 3400 and layer-thickness-to-cylinder-radius ratios of 5 and 11 were simulated. All essential turbulence scales were resolved in both calculations, and a large number of turbulence statistics were computed. The results are compared with the plane channel results of Kim et al. (1987) and with experiments. With transverse curvature the skin friction coefficient increases and the turbulence statistics, when scaled with wall units, are lower than in the plane channel. The momentum equation provides a scaling that collapses the cylinder statistics, and allows the results to be interpreted in light of the plane channel flow. The azimuthal and radial length scales of the structures in the flow are of the order of the cylinder diameter. Boomerang-shaped structures with large spanwise length scales were observed in the flow.

Neves, J. C.↗

The Effect of Solar Wind Turbulence on Parallel and Oblique Firehose Instabilities

We consider the firehose instability coexisting with the omnipresent ambient solar wind turbulence. The characteristic temporal and spatial scales of the turbulence are comparable to those of the instability. Therefore, turbulence may violate the common assumption of a uniform and stationary background used to describe instabilities and make the properties of the instabilities different. To investigate this effect, we perform threedimensional hybrid simulations with particle-in-cell ions and a quasi-neutralizing electron fluid. We find that the turbulence significantly reduces the growth rates and saturation levels of both instabilities. Comparing the cases with and without turbulence, the former results in a higher temperature anisotropy in the asymptotic marginally stable state at large times. In the former case, the distribution function averaged over the simulation box is also closer to the initial one.

S A Markovskii↗

Theory and Simulation of Real and Ideal Magnetohydrodynamic Turbulence

Incompressible, homogeneous magnetohydrodynamic (MHD) turbulence consists of fluctuating vorticity and magnetic fields, which are represented in terms of their Fourier coefficients. Here, a set of five Fourier spectral transform method numerical simulations of two-dimensional (2-D) MHD turbulence on a 512(sup 2) grid is described. Each simulation is a numerically realized dynamical system consisting of Fourier modes associated with wave vectors k, with integer components, such that k = |k| less than or equal to k(sub max). The simulation set consists of one ideal (non-dissipative) case and four real (dissipative) cases. All five runs had equivalent initial conditions. The dimensions of the dynamical systems associated with these cases are the numbers of independent real and imaginary parts of the Fourier modes. The ideal simulation has a dimension of 366104, while each real simulation has a dimension of 411712. The real runs vary in magnetic Prandtl number P(sub M), with P(sub M) is a member of {0.1, 0.25, 1, 4}. In the results presented here, all runs have been taken to a simulation time of t = 25. Although ideal and real Fourier spectra are quite different at high k, they are similar at low values of k. Their low k behavior indicates the existence of broken symmetry and coherent structure in real MHD turbulence, similar to what exists in ideal MHD turbulence. The value of PM strongly affects the ratio of kinetic to magnetic energy and energy dissipation (which is mostly ohmic). The relevance of these results to 3-D Navier-Stokes and MHD turbulence is discussed.

Shebalin, John V.↗

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 32(exp 3) simulation (Run A) and a 64(exp 3) simulation (Run B). Other than the grid size, the primary differences between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. ('Compressible motion' is that part of the velocity which has zero curl, but non-zero divergence). The numerical results indicate that gravitational collapse can proceed in a variety of ways. In Run A, collapse led to an elongated tube-like structure, while in Run B, collapse led to a flatter, disk-like structure.

Shebalin, J. V.↗

Transfer-AE: A novel autoencoder-based impact detection model for structural digital twin

Accurately detecting the location and intensity of impacts is crucial for ensuring structural safety. Currently, AI-based structural impact detection methods are widely used for their excellent detection accuracy. However, their generalization capability is limited by the scenarios present in the training data. Many complex and dangerous impact scenarios are difficult to conduct real-world experiments on to collect sufficient samples. To capture all impact scenarios and fully leverage the advantages of AI-based detection technologies, advanced methods involve combining real-world structural monitoring data with corresponding numerical models to construct digital twins. These methods continuously refine the created numerical models with limited real-world data and provide diverse impact scenarios through numerical model simulations. However, there are inevitable differences between digital models and physical models that are challenging to correct through mechanical means. This discrepancy in data distribution between the two models significantly hinders the application of digital twin technology in impact/event identification tasks. To address this challenge, this study proposes a novel model based on autoencoders, named Transfer-AE. Transfer-AE encodes the common features of digital twins in the latent space to bridge the uncertainty gap at a macro scale between numerical models and physical models and synchronously fits the magnitude and location of the impact load in the decoder. This enables consistent detection results for the same impact event, whether the sample comes from the numerical model or the physical model. Transfer-AE includes two operating modes: Mode 1 has a fixed computational complexity with stable inference speed, but the training cost and difficulty increase with data distribution. Mode 2's computational complexity increases with data distribution, but it has a fixed training cost and speed. In both cases involving the geodesic dome structure simulating a deep space habitat and the IASC-ASCE benchmark structure, Transfer-AE demonstrated the best performance in impact localization and quantification tasks compared to mainstream domain-adaptive transfer models.

Chengjia Han↗

Nonlinear dynamic simulation of single- and multi-spool core engines

In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.

Schobeiri, T.↗

Computer and laboratory simulation of interactions between spacecraft surfaces and charged-particle environments

Cases where the charged-particle environment acts on the spacecraft (e.g., spacecraft charging phenomena) and cases where a system on the spacecraft causes the interaction (e.g., high voltage space power systems) are considered. Both categories were studied in ground simulation facilities to understand the processes involved and to measure the pertinent parameters. Computer simulations are based on the NASA Charging Analyzer Program (NASCAP) code. Analytical models are developed in this code and verified against the experimental data. Extrapolation from the small test samples to space conditions are made with this code. Typical results from laboratory and computer simulations are presented for both types of interactions. Extrapolations from these simulations to performance in space environments are discussed.

Stevens, N. J.↗

Analysis of TIMS performance subjected to simulated wind blast

The results of the performance of the Thermal Infrared Multispectral Scanner (TIMS) when it is subjected to various wind conditions in the laboratory are described. Various wind conditions were simulated using a 24 inch fan or combinations of air jet streams blowing toward either or both of the blackbody surfaces. The fan was used to simulate a large volume of air flow at moderate speeds (up to 30 mph). The small diameter air jets were used to probe TIMS system response in reaction to localized wind perturbations. The maximum nozzle speed of the air jet was 60 mph. A range of wind directions and speeds were set up in the laboratory during the test. The majority of the wind tests were conducted under ambient conditions with the room temperature fluctuating no more than 2 C. The temperature of the high speed air jet was determined to be within 1 C of the room temperature. TIMS response was recorded on analog tape. Additional thermistor readouts of the blackbody temperatures and thermocouple readout of the ambient temperature were recorded manually to be compared with the housekeeping data recorded on the tape. Additional tests were conducted under conditions of elevated and cooled room temperatures. The room temperature was varied between 19.5 to 25.5 C in these tests. The calibration parameters needed for quantitative analysis of TIMS data were first plotted on a scanline-by-scanline basis. These parameters are the low and high blackbody temperature readings as recorded by the TIMS and their corresponding digitized count values. Using these values, the system transfer equations were calculated. This equation allows us to compute the flux for any video count by computing the slope and intercept of the straight line that relates the flux to the digital count. The actual video of the target (the lab floor in this case) was then compared with a simulated target. This simulated target was assumed to be a blackbody at emissivity of .95 degrees and the temperature was assumed to be at ambient temperature as recorded by the TIMS for each scanline. Using the slope and the intercept the flux corresponding to this target was converted into digital counts. The counts were observed to have a strong correlation with the actual video as recorded by the TIMS. The attached graphs describe the performance of the TIMS when compressed air is blown at each one of the blackbodies at different speeds. The effect of blowing a fan and changing the room temperature is also being analyzed. Results indicate that the TIMS system responds to variation in wind speed in real time and maintains the capability to produce accurate temperatures on a scan line basis.

Jaggi, S.↗

Large-Scale Tropospheric Transport in the Chemistry-Climate Model Initiative (CCMI) Simulations

Understanding and modeling the large-scale transport of trace gases and aerosols is important for interpreting past (and projecting future) changes in atmospheric composition. Here we show that there are large differences in the global-scale atmospheric transport properties among the models participating in the IGAC SPARC Chemistry–Climate Model Initiative (CCMI). Specifically, we find up to 40% differences in the transport timescales connecting the Northern Hemisphere (NH) midlatitude surface to the Arctic and to Southern Hemisphere high latitudes, where the mean age ranges between 1.7 and 2.6 years. We show that these differences are related to large differences in vertical transport among the simulations, in particular to differences in parameterized convection over the oceans. While stronger convection over NH midlatitudes is associated with slower transport to the Arctic, stronger convection in the tropics and subtropics is associated with faster interhemispheric transport. We also show that the differences among simulations constrained with fields derived from the same reanalysis products are as large as (and in some cases larger than) the differences among free-running simulations, most likely due to larger differences in parameterized convection. Our results indicate that care must be taken when using simulations constrained with analyzed winds to interpret the influence of meteorology on tropospheric composition.

Orbe, Clara↗

Approximation of Engine Casing Temperature Constraints for Casing Mounted Electronics

The performance of propulsion engine systems is sensitive to weight and volume considerations. This can severely constrain the configuration and complexity of the control system hardware. Distributed Engine Control technology is a response to these concerns by providing more flexibility in designing the control system, and by extension, more functionality leading to higher performing engine systems. Consequently, there can be a weight benefit to mounting modular electronic hardware on the engine core casing in a high temperature environment. This paper attempts to quantify the in-flight temperature constraints for engine casing mounted electronics. In addition, an attempt is made at studying heat soak back effects. The Commercial Modular Aero Propulsion System Simulation 40k (C-MAPSS40k) software is leveraged with real flight data as the inputs to the simulation. A two-dimensional (2-D) heat transfer model is integrated with the engine simulation to approximate the temperature along the length of the engine casing. This modification to the existing C-MAPSS40k software will provide tools and methodologies to develop a better understanding of the requirements for the embedded electronics hardware in future engine systems. Results of the simulations are presented and their implications on temperature constraints for engine casing mounted electronics is discussed.

Engine↗

Approximation of Engine Casing Temperature Constraints for Casing Mounted Electronics

The performance of propulsion engine systems is sensitive to weight and volume considerations. This can severely constrain the configuration and complexity of the control system hardware. Distributed Engine Control technology is a response to these concerns by providing more flexibility in designing the control system, and by extension, more functionality leading to higher performing engine systems. Consequently, there can be a weight benefit to mounting modular electronic hardware on the engine core casing in a high temperature environment. This paper attempts to quantify the in-flight temperature constraints for engine casing mounted electronics. In addition, an attempt is made at studying heat soak back effects. The Commercial Modular Aero Propulsion System Simulation 40k (C-MAPSS40k) software is leveraged with real flight data as the inputs to the simulation. A two-dimensional (2-D) heat transfer model is integrated with the engine simulation to approximate the temperature along the length of the engine casing. This modification to the existing C-MAPSS40k software will provide tools and methodologies to develop a better understanding of the requirements for the embedded electronics hardware in future engine systems. Results of the simulations are presented and their implications on temperature constraints for engine casing mounted electronics is discussed.

Engine↗

Orbit-averaged behavior of magnetic control laws for momentum unloading

Analytical formulas are derived for orbit-averaged behavior of magnetic control laws for unloading the excess angular momentum of a spacecraft reaction wheel control system in the presence of secular environmental torques. The specific example of an axially symmetric spacecraft with an inertially fixed attitude for which the dominant environmental torque is the gravity-gradient torque is treated in detail, but extensions of the general approach to other inertially fixed and earth-pointing spacecraft are discussed. The analytical formulas are compared to detailed simulations performed for the Solar Maximum Mission spacecraft, and agreement to within 10% is found. The analytical formulas can be used in place of detailed simulations for preliminary studies, and can be used to find selected cases giving the most stringent tests of momentum unloading capability for which detailed simulations may be performed.

Camillo, P. J.↗

Trajectory optimization and guidance law development for national aerospace plane applications

The work completed to date is comprised of the following: a simple vehicle model representative of the aerospace plane concept in the hypersonic flight regime, fuel-optimal climb profiles for the unconstrained and dynamic pressure constrained cases generated using a reduced order dynamic model, an analytic switching condition for transition to rocket powered flight as orbital velocity is approached, simple feedback guidance laws for both the unconstrained and dynamic pressure constrained cases derived via singular perturbation theory and a nonlinear transformation technique, and numerical simulation results for ascent to orbit in the dynamic pressure constrained case.

Calise, A. J.↗

Direct simulation of hypersonic rarefied flow about a delta wing

Three-dimensional hypersonic rarefied flow about a delta wing at angle of incidence has been studied using the direct simulation Monte Carlo technique. Results are obtained for a transitional flow case which has been investigated in a nitrogen wind-tunnel experiment. The simulations are performed using a nonreacting, single-species gas model that accounts for rotational and vibrational internal energies. The computations yield an attached leeside flow associated with supersonic expansion over the leeside of the wing. Results are presented for the computed flowfield and surface quantities and overall aerodynamic coefficients.

Celenligil, M. Cevdet↗

Numerical simulation of frontogenesis in a moist atmosphere

This paper describes the effects of condensation and evaporation on mesoscale frontal circulations in a two-dimensional numerical model. Utilizing an explicit scheme for the prediction of water vapor, cloud water and rainwater, the model is used to investigate the interactions between convection and the larger-scale environment. The model results are qualitatively compared with results of theoretical and observational studies, including those from the recent Severe Environmental Storms and Mesoscale Experiment-Atmospheric Variability Experiment (SESAME-AVE). Three major differences are observed in a comparison of the moist and dry simulations: (1) The speed of the upper- and lower-level jets was significantly higher in the moist case, (2) The intensity of the ageostrophic circulations in the moist simulation was much stronger, (3) The vertical velocity field in the moist case was characterized by a banded structure not present in the dry case.

Hsie, E.-Y.↗

Modeling the Impact of Drizzle and 3D Cloud Structure on Remote Sensing of Effective Radius

Remote sensing of cloud particle size with passive sensors like MODIS is an important tool for cloud microphysical studies. As a measure of the radiatively relevant droplet size, effective radius can be retrieved with different combinations of visible through shortwave infrared channels. MODIS observations sometimes show significantly larger effective radii in marine boundary layer cloud fields derived from the 1.6 and 2.1 pm channel observations than for 3.7 pm retrievals. Possible explanations range from 3D radiative transport effects and sub-pixel cloud inhomogeneity to the impact of drizzle formation on the droplet distribution. To investigate the potential influence of these factors, we use LES boundary layer cloud simulations in combination with 3D Monte Carlo simulations of MODIS observations. LES simulations of warm cloud spectral microphysics for cases of marine stratus and broken stratocumulus, each for two different values of cloud condensation nuclei density, produce cloud structures comprising droplet size distributions with and without drizzle size drops. In this study, synthetic MODIS observations generated from 3D radiative transport simulations that consider the full droplet size distribution will be generated for each scene. The operational MODIS effective radius retrievals will then be applied to the simulated reflectances and the results compared with the LES microphysics.

Platnick, Steven↗