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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Measurements of particle emission from discharge sites in Teflon irradiated by high energy electron beams

Anomalous behavior of synchronous orbit satellites manifested by overall degradation of system performance and reduced operating life is associated with electrical discharges resulting from differential charging of the spacecraft surface by fluxes of high energy electrons. During a laboratory simulation silver-backed Teflon samples have been irradiated by electron beams having energies in the range 16-26 keV. Charged particles emitted from the resultant electrical discharges have been measured with a biased Faraday cup and retarding potential analyser. Measurements indicate the presence of two distinct fluxes of particles, the first being an early pulse (0-600ns) of high energy (about 7keV) electrons, while the second is a late pulse (1-5 microseconds) of low energy electrons (less than 1eV) and ions (70eV) leaving the discharge site as a quasi plasma. Calculations indicate an electrostatic field as the dominant accelerating mechanism for charged particles.

Hazelton, R. C.↗
NASA NTRS2013

Solar Flare X-Ray Source Motion as a Response to Electron Spectral Hardening

Context. Solar flare hard X-rays (HXRs) are thought to be produced by nonthermal coronal electrons stopping in the chromosphere, or remaining trapped in the corona. The collisional thick target model (CTTM) predicts that more energetic electrons penetrate to greater column depths along the flare loop. This requires that sources produced by harder power-law injection spectra should appear further down the legs or footpoints of a flareloop. Therefore, the frequently observed hardening of the injected power-law electron spectrum during flare onset should be concurrent with a descending hard X-ray source. Aims. To test this implication of the CTTM by comparing its predicted HXR source locations with those derived from observations of a solar flare which exhibits a nonthermally-dominated spectrum before the peak in HXRs, known as an early impulsive event. Methods. HXR images and spectra of an early impulsive C-class flare were obtained using the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). Images were reconstructed to produce HXR source height evolutions for three energy bands. Spatially-integrated spectral analysis was performed to isolate nonthermal emission, and to determine the power-law index of the electron injection spectrum. The observed height-time evolutions were then fit with CTTM-based simulated heights for each energy, using the electron spectral indices derived from the RHESSI spectra. Results. The flare emission was found to be dominantly nonthermal above approx. 7 keV, with emission of thermal and nonthermal X-rays likely to be simultaneously observable below that energy. The density structure required for a good match between model and observed source heights agreed with previous studies of flare loop densities. Conclusions. The CTTM has been used to produce a descent of model HXR source heights that compares well with observations of this event. Based on this interpretation, downward motion of nonthermal sources should indeed occur in any flare where there is spectral hardening in the electron distribution during a flare. However, this would often be masked by thermal emission associated with flare plasma pre-heating. As yet, flare models that predict transfer of energy from the corona to the chromosphere by means other than a flux of nonthermal electrons do not predict this observed source descent. Therefore, flares such as this will be key in explaining this elusive energy transfer process.

X-rays↗
NASA NTRS2015

Modeling a Full Coronal Loop Observed with Hinode EIS and SDO AIA

Physical parameters measured from an observation of a coronal loop from Gupta et al. (2015) using Hinode/EIS and SDO/AIA were used as input for the hydrodynamic, impulsively heating NRLSOFM 1-­d loop model. The model was run at eight different energy inputs and used the measured quantities of temperature (0.73 MK), density (10(sup 8.5)cm(sup -3) and minimum loop lifetime to evaluate the success of the model at recreating the observations. The loop was measured by us to have an unprojected length of 236 Mm and was assumed to be almost perpendicular to the solar surface (tilt of 3.5 degrees) and have a dipolar geometry. Our results show that two of our simulation runs (with input energies of 0.01 and 0.02 ergs cm(sup -3)S(sup -1) closely match the temperature/density combination exhibited by the loop observation. However, our simulated loops only remain in the temperature sensitive region of the Mg 278.4 Angstrom filter for 500 and 800 seconds respectively which is less than the 1200 seconds that the loop is observed for with EIS in order to make the temperature/density measurements over the loop's entire length. This leads us to conclude that impulsive heating of a single loop is not complex enough to explain this observation. Additional steady heating or a collection of additional strands along the line-­of-­sight would help to align the simulation with the observation.

coronal loop↗
NASA NTRS2003

Computational Methods for Jet Noise Simulation

The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Goodrich, John W.↗
NASA NTRS2017

Simulation of Galactic Cosmic Rays and Dose-Rate Effects in RITRACKS

The NASA Space Radiation Laboratory (NSRL) facility has been used successfully for many years to generate ion beams for radiation research experiments by NASA investigators. Recently, modifications were made to the beam lines to allow rapid switching between different types of ions and energies, with the aim to simulate the Galactic Cosmic Rays (GCR) environment. As this will be a focus of space radiation research for upcoming years, the stochastic radiation track structure code RITRACKS (Relativistic Ion Tracks) was modified to simulate beams of various ion types and energies during time intervals specified by the user at the microscopic and nanoscopic scales. For example, particle distributions of a mixed 344.1-MeV protons (18.04 cGy) and 950-MeV/n iron (5.64 cGy) beam behind a 20 g/cm(exp 2) aluminum followed by a 10 g/cm(exp 2) polyethylene shield as calculated by the code GEANT4 were used as an input field in RITRACKS. Similarly, modifications were also made to simulate a realistic radiation environment in a spacecraft exposed to GCR by sampling the ion types and energies from particle spectra pre-calculated by the code HZETRN. The newly implemented features allows RITRACKS to generate time-dependent differential and cumulative 3D dose voxel maps. These new capabilities of RITRACKS will be used to investigate dose-rate effects and synergistic interactions of various types of radiations for many end points at the microscopic and nanoscopic scales such as DNA damage and chromosome aberrations.

Plante, Ianik↗
NASA NTRS1984

Nonthermal synchrotron radiation in a strong magnetic field

A synchrotron source model is presented for the origin of gamma-ray bursts from neutron stars. The field strength and length of time needed for electrons to lose sufficient energy in emitting the gamma rays are calculated, together with the efficiency of the process, limiting frequencies, and the energy output. A Monte Carlo simulation of the synchrotron radiation spectrum of an electron undergoing a self-consistent energy loss demonstrates that a process involving thermalized electrons can be sufficiently efficient to produce 30 keV gamma rays. Problems expected in using the current capabilities of detectors to scan for the projected energy spectra from fields of at least 3 x 10 to the 12th Gauss are discussed.

Bussard, R. W.↗
NASA NTRS2005

Microwave Medical Treatment Apparatus and Method

Methods, simulations, and apparatus are provided that may be utilized for medical treatments which are especially suitable for treatment of benign prostatic hyperplasia (BPH). In a preferred embodiment, a plurality of separate microwave antennas are utilized to heat prostatic tissue to promote necrosing of the prostatic tissue that relieves the pressure of the prostatic tissue against the urethra as the body reabsorbs the necrosed or dead tissue. By utilizing constructive and destructive interference of the microwave transmission, the energy can be deposited on the tissues to be necrosed while protecting other tissues such as the urethra. Saline injections to alter the conductivity of the tissues may also be used to further focus the energy deposits. A computer simulation is Provided that can be used to Predict the resulting temperature profile produced in the prostatic tissue. By changing the various control features of one or more catheters and the methods of applying microwave energy, a temperature profile can be predicted and produced that is similar to the temperature profile desired for the particular patient.

Arndt, G. Dickey↗
NASA NTRS2022

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗
NASA NTRS1986

Magnetic pumping by magnetosonic waves in the presence of noncompressive electromagnetic fluctuations

Numerical simulations of the damping of magnetosonic waves via magnetic pumping in the presence of electromagnetic fluctuations that can pitch-angle scatter the plasma particles are presented. From the first simulation it is found that the magnetosonic-wave energy is transferred to high-energy particles. In the second type of simulation, magnetosonic waves produce a hot surface layer on the plasma that is ablated by the wave energy. Solution of a Fokker-Planck equation for the magnetic-pumping process is found to adequately represent magnetic pumping by small-amplitude magnetosonic waves.

Borovsky, Joseph E.↗
NASA NTRS1992

Phase change energy storage for solar dynamic power systems

This paper presents the results of a transient computer simulation that was developed to study phase change energy storage techniques for Space Station Freedom (SSF) solar dynamic (SD) power systems. Such SD systems may be used in future growth SSF configurations. Two solar dynamic options are considered in this paper: Brayton and Rankine. Model elements consist of a single node receiver and concentrator, and takes into account overall heat engine efficiency and power distribution characteristics. The simulation not only computes the energy stored in the receiver phase change material (PCM), but also the amount of the PCM required for various combinations of load demands and power system mission constraints. For a solar dynamic power system in low earth orbit, the amount of stored PCM energy is calculated by balancing the solar energy input and the energy consumed by the loads corrected by an overall system efficiency. The model assumes an average 75 kW SD power system load profile which is connected to user loads via dedicated power distribution channels. The model then calculates the stored energy in the receiver and subsequently estimates the quantity of PCM necessary to meet peaking and contingency requirements. The model can also be used to conduct trade studies on the performance of SD power systems using different storage materials.

Chiaramonte, F. P.↗
NASA NTRS2014

Bidirectional Energy Cascades and the Origin of Kinetic Alfvenic and Whistler Turbulence in the Solar Wind

The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore how the free energy released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfvénic and whistler turbulence are excited that evolve through inverse and forward magnetic energy cascades.

plasmas↗
NASA NTRS1991

Performance characteristics of a thermal energy storage module - A transient PCM/forced convection conjugate analysis

The performance of a thermal energy storage module is simulated numerically. The change of phase of the phase-change material (PCM) and the transient forced convective heat transfer for the transfer fluid with low Prandtl numbers are solved simultaneously as a conjugate problem. A parametric study and a system optimization are conducted. The numerical results show that module geometry is crucial to the design of a space-based thermal energy storage system.

Cao, Y.↗
NASA NTRS2001

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗
NASA NTRS1976

Conceptual design of a parabolic dish solar collector using simulation techniques

The development of solar concentrators in recent years have produced a wide variety of collectors for the utilization of solar energy. This paper presents the simulation techniques used to predict the optical and thermal performance of a paraboloid of revolution type solar collector. Conceptual design of a dish concentrator with a fixed receiver size is obtained by parametrically examining the significant variables.

Gupta, B. P.↗
NASA NTRS1983

Simulating sunflower canopy temperatures to infer root-zone soil water potential

A soil-plant-atmosphere model for sunflower (Helianthus annuus L.), together with clear sky weather data for several days, is used to study the relationship between canopy temperature and root-zone soil water potential. Considering the empirical dependence of stomatal resistance on insolation, air temperature and leaf water potential, a continuity equation for water flux in the soil-plant-atmosphere system is solved for the leaf water potential. The transpirational flux is calculated using Monteith's combination equation, while the canopy temperature is calculated from the energy balance equation. The simulation shows that, at high soil water potentials, canopy temperature is determined primarily by air and dew point temperatures. These results agree with an empirically derived linear regression equation relating canopy-air temperature differential to air vapor pressure deficit. The model predictions of leaf water potential are also in agreement with observations, indicating that measurements of canopy temperature together with a knowledge of air and dew point temperatures can provide a reliable estimate of the root-zone soil water potential.

Choudhury, B. J.↗
NASA NTRS1987

Helicity fluctuations in incompressible turbulent flows

Results from direct numerical simulations of several homogeneous flows and fully developed turbulent channel flow indicate that the probability distribution function (pdf) of relative helicity density exhibits at most a 20 percent deviation from a flat distribution. Isotropic flows exhibit a slight helical nature but the presence of mean strain in homogeneous turbulence suppresses helical behavior. All the homogeneous turbulent flows studied show no correlation between relative helicity density and the dissipation of turbulent kinetic energy. The channel flow simulations indicate that, except for low-dissipation regions near the outer edge of the buffer layer, there is no tendency for the flow to be helical. The strong peaks in the relative helicity density pdf and the association of these peaks with regions of low dissipation found in previous simulations by Pelz et al.(1985) are not observed.

Rogers, Michael M.↗
NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗
NASA NTRS2018

Heasim and Skyback Simulation Tools and Their Application to the Hitomi Mission

We present an introduction to the heasim multimission observation and skyback background, high-energy pseudo Monte Carlo astrophysical simulation tools. Heasim may be used to accurately and efficiently construct flexible image transport system (FITS) event files for simple or composite sources with a wide range of standard and user-defined spatial, spectral, and temporal characteristics. Skyback is designed to enable users to assess the impact of background discrete and diffuse emission on prospective observations, and skyback output may be directly input into heasim. We present a brief overview of heasim and skyback input, algorithms, usage, and output. We also introduce the sxsbranch tool that computes Hitomi soft X-ray spectrometer resolution grade branching ratios, emphasizing its application to simulations. We include several examples of particular relevance to the Hitomi mission.

Loewenstein, Michael↗