Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy simulation”

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 505 records · Page 28

Subgrid-scale backscatter in transitional and turbulent flows

Most subgrid-scale (SGS) models for large-eddy simulations are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Direct numerical simulations of turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomena. In all flows considered roughly 50 percent of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type and Reynolds number. The amount of backscatter showed an increasing trend with Reynolds numbers. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fluctuation Mach number.

Piomelli, Ugo↗

Influence of 2D Finite Element Modeling Assumptions on Debonding Prediction for Composite Skin-stiffener Specimens Subjected to Tension and Bending

The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane strain elements as well as three different generalized plane strain type approaches were performed. The computed deflections, skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with lamination length. For more accurate predictions, however, a three-dimensional analysis is required.

Krueger, Ronald↗

Comparison of 2D Finite Element Modeling Assumptions with Results From 3D Analysis for Composite Skin-Stiffener Debonding

The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane-strain elements as well as three different generalized plane strain type approaches were performed. The computed skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with delamination length. For more accurate predictions, however, a three-dimensional analysis is required.

Krueger, Ronald↗

Assimilation of GOES Land Surface Data into a Mesoscale Models

A technique has been developed for assimilating Geostationary Operational Environmental Satellite (GOES)-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The assimilation technique has been applied to the Oklahoma-Kansas region during the spring-summer 2000 time period when dynamic changes in vegetation cover occur. In April, central Oklahoma is characterized by large NDVI associated with winter wheat while surrounding areas are primarily rangeland with lower NDVI. In July the vegetation pattern reverses as the central wheat area changes to low NDVI due to harvesting and the surrounding rangeland is greener than it was in April. The goal of this study is to determine if assimilating satellite land surface data can improve simulation of the complex spatial distribution of surface energy and water fluxes across this region. The PSU/NCAR NM5 V3 system is used in this study. The grid configuration consists of a 36-km CONUS domain and a 12-km nest over the area of interest. Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint indicates that assimilation of the satellite data results reduces both the bias and RMSE for both state variables. In addition, comparison of model data with ARM/CART EBBR flux observations reveals that the assimilation technique adjusts the bowen ratio in a realistic fashion.

Lapenta, William M.↗

The simulated features of heliospheric cosmic-ray modulation with a time-dependent drift model. III - General energy dependence

The time-dependent cosmic-ray transport equation is solved numerically in an axially symmetric heliosphere. Gradient and curvature drifts are incorporated, together with an emulated wavy neutral sheet. This model is used to simulate heliospheric cosmic-ray modulation for the period 1985-1989 during which drifts are considered to be important. The general energy dependence of the modulation of Galactic protons is studied as predicted by the model for the energy range 1 MeV to 10 GeV. The corresponding instantaneous radial and latitudinal gradients are calculated, and it is found that, whereas the latitudinal gradients follow the trends in the waviness of the neutral sheet to a large extent for all energies, the radial gradients below about 200 MeV deviate from this general pattern. In particular, these gradients increase when the waviness decreases for the simulated period 1985-1987.3, after which they again follow the neutral sheet by increasing rapidly.

Potgieter, M. S.↗

The ideal Kolmogorov inertial range and constant

The energy transfer statistics measured in numerically simulated flows are found to be nearly self-similar for wavenumbers in the inertial range. Using the measured self-similar form, an 'ideal' energy transfer function and the corresponding energy flux rate were deduced. From this flux rate, the Kolmogorov constant was calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗

Preliminary results from a four-working space, double-acting piston, Stirling engine controls model

A four working space, double acting piston, Stirling engine simulation is being developed for controls studies. The development method is to construct two simulations, one for detailed fluid behavior, and a second model with simple fluid behaviour but containing the four working space aspects and engine inertias, validate these models separately, then upgrade the four working space model by incorporating the detailed fluid behaviour model for all four working spaces. The single working space (SWS) model contains the detailed fluid dynamics. It has seven control volumes in which continuity, energy, and pressure loss effects are simulated. Comparison of the SWS model with experimental data shows reasonable agreement in net power versus speed characteristics for various mean pressure levels in the working space. The four working space (FWS) model was built to observe the behaviour of the whole engine. The drive dynamics and vehicle inertia effects are simulated. To reduce calculation time, only three volumes are used in each working space and the gas temperature are fixed (no energy equation). Comparison of the FWS model predicted power with experimental data shows reasonable agreement. Since all four working spaces are simulated, the unique capabilities of the model are exercised to look at working fluid supply transients, short circuit transients, and piston ring leakage effects.

Daniele, C. J.↗

Crops Models for Varying Environmental Conditions

New variable environment Modified Energy Cascade (MEC) crop models were developed for all the Advanced Life Support (ALS) candidate crops and implemented in SIMULINK. The MEC models are based on the Volk, Bugbee, and Wheeler Energy Cascade (EC) model and are derived from more recent Top-Level Energy Cascade (TLEC) models. The MEC models simulate crop plant responses to day-to-day changes in photosynthetic photon flux, photoperiod, carbon dioxide level, temperature, and relative humidity. The original EC model allows changes in light energy but uses a less accurate linear approximation. The simulation outputs of the new MEC models for constant nominal environmental conditions are very similar to those of earlier EC models that use parameters produced by the TLEC models. There are a few differences. The new MEC models allow setting the time for seed emergence, have realistic exponential canopy growth, and have corrected harvest dates for potato and tomato. The new MEC models indicate that the maximum edible biomass per meter squared per day is produced at the maximum allowed carbon dioxide level, the nominal temperatures, and the maximum light input. Reducing the carbon dioxide level from the maximum to the minimum allowed in the model reduces crop production significantly. Increasing temperature decreases production more than it decreases the time to harvest, so productivity in edible biomass per meter squared per day is greater at nominal than maximum temperatures, The productivity in edible biomass per meter squared per day is greatest at the maximum light energy input allowed in the model, but the edible biomass produced per light energy input unit is lower than at nominal light levels. Reducing light levels increases light and power use efficiency. The MEC models suggest we can adjust the light energy day-to- day to accommodate power shortages or Lise excess power while monitoring and controlling edible biomass production.

Jones, Harry↗

Computation of large-scale statistics in decaying isotropic turbulence

We have performed large-eddy simulations of decaying isotropic turbulence to test the prediction of self-similar decay of the energy spectrum and to compute the decay exponents of the kinetic energy. In general, good agreement between the simulation results and the assumption of self-similarity were obtained. However, the statistics of the simulations were insufficient to compute the value of gamma which corrects the decay exponent when the spectrum follows a k(exp 4) wave number behavior near k = 0. To obtain good statistics, it was found necessary to average over a large ensemble of turbulent flows.

Chasnov, Jeffrey R.↗

Design of the Experimental Exposure Conditions to Simulate Ionizing Radiation Effects on Candidate Replacement Materials for the Hubble Space Telescope

In this effort, experimental exposure times for monoenergetic electrons and protons were determined to simulate the space radiation environment effects on Teflon components of the Hubble Space Telescope. Although the energy range of the available laboratory particle accelerators was limited, optimal exposure times for 50 keV, 220 keV, 350 keV, and 500 KeV electrons were calculated that produced a dose-versus-depth profile that approximated the full spectrum profile, and were realizable with existing equipment. For the case of proton exposure, the limited energy range of the laboratory accelerator restricted simulation of the dose to a depth of .5 mil. Also, while optimal exposure times were found for 200 keV, 500 keV and 700 keV protons that simulated the full spectrum dose-versus-depth profile to this depth, they were of such short duration that the existing laboratory could not be controlled to within the required accuracy. In addition to the obvious experimental issues, other areas exist in which the analytical work could be advanced. Improved computer codes for the dose prediction- along with improved methodology for data input and output- would accelerate and make more accurate the calculational aspects. This is particularly true in the case of proton fluxes where a paucity of available predictive software appears to exist. The dated nature of many of the existing Monte Carlo particle/radiation transport codes raises the issue as to whether existing codes are sufficient for this type of analysis. Other areas that would result in greater fidelity of laboratory exposure effects to the space environment is the use of a larger number of monoenergetic particle fluxes and improved optimization algorithms to determine the weighting values.

Smith, L. Montgomery↗

Study of Kapton Degradation under Simulated Shuttle Environment

Weight loss and severe degradation of the surface of Kapton that occurs in low Earth orbit is studied. Atomic oxygen, the major ambient species at low Earth altitude and incident with approximately 5 eV energy in ram conditions, is the primary suspect, but a thorough study of oxygen-Kapton interactions has not yet been carried out. A low-energy ion source is used to simulate the shuttle low Earth orbit environment. This source, together with diagnostic tools including surface analysis and mass spectroscopic capability, is being used to carry out experiments from which quantum yields may be obtained.

Eck, T. G.↗

Battery simulation program

Computer program calculates spacecraft battery energy at specific times dictated by input sequence of recharge, discharge, and no activity phases. Recharge phases are assumed during times not specified, unless batteries are fully charged. Warnings are printed when energy falls below specified level. Program assumes two identical batteries.

Abreu, A.↗

Altitude Effects on Thermal Ice Protection System Performance; A Study of an Alternative Simulation Approach

The quest for more energy-efficient green aircraft, dictates that all systems, including the ice protection system (IPS), be closely examined for ways to reduce energy consumption and to increase efficiency. A thermal ice protection systems must protect the aircraft from the hazardous effects of icing, and yet it needs to do so as efficiently as possible. The system can no longer be afforded the degree of over-design in power usage they once were. To achieve these more exacting designs, a better understanding of the heat and mass transport phenomena involved during an icing encounter is needed.

thermal ice protection system↗

Probing the High-energy Gamma-Ray Emission Mechanism in the Vela Pulsar via Phase-resolved Spectral and Energy-dependent Light-curve Modeling

Recent kinetic simulations sparked a debate regarding the emission mechanism responsible for pulsed GeV γ-ray emission from pulsars. Some models invoke curvature radiation, while other models assume synchrotron radiation in the current sheet. We interpret the curved spectrum of the Vela pulsar as seen by H.E.S.S. II (up to ∼100 GeV) and the Fermi Large Area Telescope to be the result of curvature radiation due to primary particles in the pulsar magnetosphere and current sheet. We present phase-resolved spectra and energy-dependent light curves using an extended slot gap and current-sheet model, invoking a step function for the accelerating electric field as motivated by kinetic simulations. We include a refined calculation of the curvature radius of particle trajectories in the lab frame, impacting the particle transport, predicted light curves, and spectra. Our model reproduces the decrease of the flux of the first peak relative to the second one, evolution of the bridge emission, near-constant phase positions of peaks, and narrowing of pulses with increasing energy. We can explain the first of these trends because we find that the curvature radii of the particle trajectories in regions where the second γ-ray light-curve peak originates are systematically larger than those associated with the first peak, implying that the spectral cutoff of the second peak is correspondingly larger. However, an unknown azimuthal dependence of the E field, as well as uncertainty in the precise spatial origin of the GeV emission, precludes a simplistic discrimination of emission mechanisms.

Monica Barnard↗

A numerical study of fundamental shock noise mechanisms

The results of this thesis demonstrate that direct numerical simulation can predict sound generation in unsteady aerodynamic flows containing shock waves. Shock waves can be significant sources of sound in high speed jet flows, on helicopter blades, and in supersonic combustion inlets. Direct computation of sound permits the prediction of noise levels in the preliminary design stage and can be used as a tool to focus experimental studies, thereby reducing cost and increasing the probability of a successfully quiet product in less time. This thesis reveals and investigates two mechanisms fundamental to sound generation by shocked flows: shock motion and shock deformation. Shock motion is modeled by the interaction of a sound wave with a shock. During the interaction, the shock wave begins to move and the sound pressure is amplified as the wave passes through the shock. The numerical approach presented in this thesis is validated by the comparison of results obtained in a quasi-one dimensional simulation with linear theory. Analysis of the perturbation energy demonstrated for the first time that acoustic energy is generated by the interaction. Shock deformation is investigated by the numerical simulation of a ring vortex interacting with a shock. This interaction models the passage of turbulent structures through the shock wave. The simulation demonstrates that both acoustic waves and contact surfaces are generated downstream during the interaction. Analysis demonstrates that the acoustic wave spreads cylindrically, that the sound intensity is highly directional, and that the sound pressure level increases significantly with increasing shock strength. The effect of shock strength on sound pressure level is consistent with experimental observations of shock noise, indicating that the interaction of a ring vortex with a shock wave correctly models a dominant mechanism of shock noise generation.

Meadows, Kristine R.↗

Sensitivity analysis of cool-down strategies for a transonic cryogenic tunnel

Guidelines and suggestions substantiated by real-time simulation data to ensure optimum time and energy use of injected liquid nitrogen for cooling the Langley 0.3-Meter Transonic Cryogenic Tunnel (TCT) are presented. It is directed toward enabling operators and researchers to become cognizant of criteria for using the 0.3-m TCT in an energy- or time-efficient manner. Operational recommendations were developed based on information collected from a validated simulator of the 0.3-m TCT and experimental data from the tunnel. Results and trends, however, can be extrapolated to other similarly constructed cryogenic wind tunnels.

Thibodeaux, J. J.↗

A potential enstrophy and energy conserving scheme for the shallow water equations

To improve the simulation of nonlinear aspects of the flow over steep topography, a potential enstrophy and energy conserving scheme for the shallow water equations is derived. It is pointed out that a family of schemes can conserve total energy for general flow and potential enstrophy for flow with no mass flux divergence. The newly derived scheme is a unique member of this family, that conserves both potential enstrophy and energy for general flow. Comparison by means of numerical experiment with a scheme that conserves (potential) enstrophy for purely horizontal nondivergent flow demonstrated the considerable superiority of the newly derived potential enstrophy and energy conserving scheme, not only in suppressing a spurious energy cascade but also in determining the overall flow regime. The potential enstrophy and energy conserving scheme for a spherical grid is also presented.

Arakawa, A.↗