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

The magnetotail boundary and energy transfer processes

A particle code is used to simulate the magnetopause region in the high latitude geomagnetic tail in which the magnetic field undergoes a significant increase in going from the magnetosheath to the magnetotail lobe. The simulation indicates that plasma can flow from the magnetosheath to the lobe, which is accompanied by a drop in pressure and density. In the earth's inertial frame, the particles do work against the convection electric field. Hence the magnetopause region serves as a dynamo. The simulation also shows that the width of the transition region increases with time. In the earth's inertial frame this is seen as an expansion of the magnetopause thickness in the antisunward direction.

Swift, D. W.↗

Cooperative control of two arms in the transport of an inertial load in zero gravity

In designing a robot control system for dual arm configurations, the control engineer is faced with two challenges: to derive the equations of motion for a given situation, and to meet certain desired control requirements (for instance, minimum energy). The former may involve closed kinematic chains, such as the case when the two arms are grasping a common object. The latter usually involves nonlinear optimization. These issues are considered in the context of transporting an inertial load using two planar three-link arms. A generalized 'reduction transformation' is applied to the dynamics to remove the singularity in the system equations. A suboptimal minimum energy method is presented to reduce a difficult 12-state, six-control nonlinear optimization to two independent, nonconflicting suboptimizations. A simulation example is provided to illustrate the degree of energy reduction possible using the optimal arm torque distribution that was developed.

Carignan, Craig R.↗

GCR Simulator Reference Field and a Spectral Approach for Laboratory Simulation

The galactic cosmic ray (GCR) simulator at the NASA Space Radiation Laboratory (NSRL) is intended to deliver the broad spectrum of particles and energies encountered in deep space to biological targets in a controlled laboratory setting. In this work, certain aspects of simulating the GCR environment in the laboratory are discussed. Reference field specification and beam selection strategies at NSRL are the main focus, but the analysis presented herein may be modified for other facilities. First, comparisons are made between direct simulation of the external, free space GCR field and simulation of the induced tissue field behind shielding. It is found that upper energy constraints at NSRL limit the ability to simulate the external, free space field directly (i.e. shielding placed in the beam line in front of a biological target and exposed to a free space spectrum). Second, variation in the induced tissue field associated with shielding configuration and solar activity is addressed. It is found that the observed variation is likely within the uncertainty associated with representing any GCR reference field with discrete ion beams in the laboratory, given current facility constraints. A single reference field for deep space missions is subsequently identified. Third, an approach for selecting beams at NSRL to simulate the designated reference field is presented. Drawbacks of the proposed methodology are discussed and weighed against alternative simulation strategies. The neutron component and track structure characteristics of the simulated field are discussed in this context.

Slaba, Tony C.↗

In Vivo Simulator for Microwave Treatment

Method and apparatus are provided for propagating microwave energy into heart tissues to produce a desired temperature profile therein at tissue depths sufficient for thermally ablating arrhythmogenic cardiac tissue to treat ventricular tachycardia and other arrhythmias while preventing excessive heating of surrounding tissues, organs, and blood. A wide bandwidth double-disk antenna is effective for this purpose over a bandwidth of about 6 GHz. A computer simulation provides initial screening capabilities for an antenna such as antenna. frequency, power level, and power application duration. The simulation also allows optimization of techniques for specific patients or conditions. In operation, microwave energy between about 1 GHz and 12 GHz is applied to monopole microwave radiator having a surface wave limiter. A test setup provides physical testing of microwave radiators to determine the temperature profile created in actual heart tissue or ersatz heart tissue. Saline solution pumped over the heart tissue with a peristaltic pump simulates blood flow. Optical temperature sensors disposed at various tissue depths within the heart tissue detect the temperature profile without creating any electromagnetic interference. The method may be used to produce a desired temperature profile in other body tissues reachable by catheter such as tumors and the like.

Arndt, G. Dickey↗

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗

Installation and Characterization of Charged Particle Sources for Space Environmental Effects Testing

Charged particle sources are integral devices used by Marshall Space Flight Center s Environmental Effects Branch (EM50) in order to simulate space environments for accurate testing of materials and systems. By using these sources inside custom vacuum systems, materials can be tested to determine charging and discharging properties as well as resistance to sputter damage. This knowledge can enable scientists and engineers to choose proper materials that will not fail in harsh space environments. This paper combines the steps utilized to build a low energy electron gun (The "Skevington 3000") as well as the methods used to characterize the output of both the Skevington 3000 and a manufactured Xenon ion source. Such characterizations include beam flux, beam uniformity, and beam energy. Both sources were deemed suitable for simulating environments in future testing.

Skevington, Jennifer L.↗

Phenolic Polymer Interactions with Water and Ethylene Glycol Solvents

Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final post-cured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through ab initio techniques and molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H2O. Three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H2O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H2O, and other similar solvents.

Haskins, Justin B.↗

Orion Landing Simulation Eight Soil Model Comparison

LS-DYNA finite element simulations of a rigid Orion Crew Module (CM) were used to investigate the CM impact behavior on eight different soil models. Ten different landing conditions, characterized by the combination of CM vertical and horizontal velocity, hang angle, and roll angle were simulated on the eight different soils. The CM center of gravity accelerations, pitch angle, kinetic energy, and soil contact forces were the outputs of interest. The simulation results are presented, with comparisons of the CM behavior on the different soils. The soils analyzed in this study can be roughly categorized as soft, medium, or hard, according to the CM accelerations that occur when landing on them. The soft group is comprised of the Carson Sink Wet soil and the Kennedy Space Center (KSC) Low Density Dry Sand. The medium group includes Carson Sink Dry, the KSC High Density In-Situ Moisture Sand and High Density Flooded Sand, and Cuddeback B. The hard soils are Cuddeback A and the Gantry Unwashed Sand. The softer soils were found to produce lower peak accelerations, have more stable pitch behavior, and to be less sensitive to the landing conditions. This investigation found that the Cuddeback A soil produced the highest peak accelerations and worst stability conditions, and that the best landing performance was achieved on the KSC Low Density Dry Sand.

Mark, Stephen D.↗

Real-valued composite filters for correlation-based optical pattern recognition

Advances in the technology of optical devices such as spatial light modulators (SLMs) have influenced the research and growth of optical pattern recognition. In the research leading to this report, the design of real-valued composite filters that can be implemented using currently available SLMs for optical pattern recognition and classification was investigated. The design of real-valued minimum average correlation energy (RMACE) filter was investigated. Proper selection of the phase of the output response was shown to reduce the correlation energy. The performance of the filter was evaluated using computer simulations and compared with the complex filters. It was found that the performance degraded only slightly. Continuing the above investigation, the design of a real filter that minimizes the output correlation energy and the output variance due to noise was developed. Simulation studies showed that this filter had better tolerance to distortion and noise compared to that of the RMACE filter. Finally, the space domain design of RMACE filter was developed and implemented on the computer. It was found that the sharpness of the correlation peak was slightly reduced but the filter design was more computationally efficient than the complex filter.

Rajan, P. K.↗

Energetics of a sudden stratospheric warming simulated with a three-dimensional, spectral, quasi-geostrophic model

The energetics of a three dimensional, quasi-geostrophic simulation of a sudden stratospheric warming which developed spontaneously during an annual cycle simulation are described. Daily values of the stratospheric temperatures, zonal winds, heating rates, energies, and energy conversions are discussed and compared with those for observed warmings. It is shown that, like observed warmings, the simulated warming was preceded by an increased vertical flux of eddy kinetic energy from the troposphere and the polar heating resulted because of the strong convergence of the horizontal, eddy heat flux which was only partially balanced by adiabatic and diabatic cooling. There is a significant similarity between the energetics of the simulated and observed warmings. In addition, the warming was spontaneous and the model did not develop a major warming in each winter of the simulation. These facts suggest that this model may be useful for studying not only the warming process but also the conditions that favor its development.

Haggard, K. V.↗

Multifractal spectra in homogeneous shear flow

Employing numerical simulations of 3-D homogeneous shear flow, the associated multifractal spectra of the energy dissipation, scalar dissipation and vorticity fields were calculated. The results for (128) cubed simulations of this flow, and those obtained in recent experiments that analyzed 1- and 2-D intersections of atmospheric and laboratory flows, are in some agreement. A two-scale Cantor set model of the energy cascade process which describes the experimental results from 1-D intersections quite well, describes the 3-D results only marginally.

Deane, A. E.↗

A Computational Methodology for Simulating Thermal Loss Testing of the Advanced Stirling Convertor

The U.S. Department of Energy (DOE) and Lockheed Martin Space Systems Company (LMSSC) have been developing the Advanced Stirling Radioisotope Generator (ASRG) for use as a power system for space science missions. This generator would use two highefficiency Advanced Stirling Convertors (ASCs), developed by Sunpower Inc. and NASA Glenn Research Center (GRC). The ASCs convert thermal energy from a radioisotope heat source into electricity. As part of ground testing of these ASCs, different operating conditions are used to simulate expected mission conditions. These conditions require achieving a particular operating frequency, hot end and cold end temperatures, and specified electrical power output for a given net heat input. In an effort to improve net heat input predictions, numerous tasks have been performed which provided a more accurate value for net heat input into the ASCs, including the use of multidimensional numerical models. Validation test hardware has also been used to provide a direct comparison of numerical results and validate the multi-dimensional numerical models used to predict convertor net heat input and efficiency. These validation tests were designed to simulate the temperature profile of an operating Stirling convertor and resulted in a measured net heat input of 244.4 W. The methodology was applied to the multi-dimensional numerical model which resulted in a net heat input of 240.3 W. The computational methodology resulted in a value of net heat input that was 1.7 percent less than that measured during laboratory testing. The resulting computational methodology and results are discussed.

Reid, Terry V.↗

Diffusion on Cu surfaces

Understanding surface diffusion is essential in understanding surface phenomena, such as crystal growth, thin film growth, corrosion, physisorption, and chemisorption. Because of its importance, various experimental and theoretical efforts have been directed to understand this phenomena. The Field Ion Microscope (FIM) has been the major experimental tool for studying surface diffusion. FIM have been employed by various research groups to study surface diffusion of adatoms. Because of limitations of the FIM, such studies are only limited to a few surfaces: nickel, platinum, aluminum, iridium, tungsten, and rhodium. From the theoretical standpoint, various atomistic simulations are performed to study surface diffusion. In most of these calculations the Embedded Atom Method (EAM) along with the molecular static (MS) simulation are utilized. The EAM is a semi-empirical approach for modeling the interatomic interactions. The MS simulation is a technique for minimizing the total energy of a system of particles with respect to the positions of its particles. One of the objectives of this work is to develop the EAM functions for Cu and use them in conjunction with the molecular static (MS) simulation to study diffusion of a Cu atom on a perfect as well as stepped Cu(100) surfaces. This will provide a test of the validity of the EAM functions on Cu(100) surface and near the stepped environments. In particular, we construct a terrace-ledge-kink (TLK) model and calculate the migration energies of an atom on a terrace, near a ledge site, near a kink site, and going over a descending step. We have also calculated formation energies of an atom on the bare surface, a vacancy in the surface, a stepped surface, and a stepped-kink surface. Our results are compared with the available experimental and theoretical results.

Karimi, Majid↗

On the Modeling of Shells in Multibody Dynamics

Energy preserving/decaying schemes are presented for the simulation of the nonlinear multibody systems involving shell components. The proposed schemes are designed to meet four specific requirements: unconditional nonlinear stability of the scheme, a rigorous treatment of both geometric and material nonlinearities, exact satisfaction of the constraints, and the presence of high frequency numerical dissipation. The kinematic nonlinearities associated with arbitrarily large displacements and rotations of shells are treated in a rigorous manner, and the material nonlinearities can be handled when the, constitutive laws stem from the existence of a strain energy density function. The efficiency and robustness of the proposed approach is illustrated with specific numerical examples that also demonstrate the need for integration schemes possessing high frequency numerical dissipation.

Bauchau, Olivier A.↗

Massive Stars and the Energy Balance of the Interstellar Medium: The Impact of an Isolated 60 M. Star - 1

We present results of numerical simulations carried out with a two-dimensional radiation hydrodynamics code in order to study the impact of massive stars on their surrounding interstellar medium. This first paper deals with the evolution of the circumstellar gas around an isolated 60 M. star. The interaction of the photo- ionized H II region with the stellar wind bubble forms a variety of interesting structures like shells, clouds, fingers, and spokes. These results demonstrate that complex structures found in H II regions are not necessarily relics from the time before the gas became ionized but may result from dynamical processes during the course of the H II region evolution. We have also analyzed the transfer and deposit of the stellar wind and radiation energy into the circumstellar medium until the star explodes as a supernova. Although the total mechanical wind energy supplied by the star is negligible compared to the accumulated energy of the Lyman continuum photons, the kinetic energy imparted to the circumstellar gas over the star s lifetime is 4 times higher than for a comparable windless simulation. Furthermore, the thermal energy of warm photoionized gas is lower by some 55%). Our results document the necessity to consider both ionizing radiation and stellar winds for an appropriate description of the interaction of OB stars with their circumstellar environment.

Freyer, Tim↗

Laboratory electron exposure of TSS-1 thermal control coating

RM400, a conductive thermal control coating, was developed for use on the exterior shell of the tethered satellite. Testing was performed by the Engineering Physics Division to quantify effects of the space environment on this coating and its conductive and optical properties. Included in this testing was exposure of RM400 to electrons with energies ranging from 0.1 to 1 keV, to simulate electrons accelerated from the ambient space plasma when the tethered satellite is fully deployed. During this testing, the coating was found to luminesce, and a prolonged exposure of the coating to high-energy electrons caused the coating to darken. This report describes the tests done to quantify the degradation of the thermal control properties caused by electron exposure and to measure the luminescence as a function of electron energy and current density to the satellite.

Vaughn, J. A.↗

Neutralizer Hollow Cathode Simulations and Comparisons with Ground Test Data

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially ionized gas in these devices. A summary of the comparisons between simulation results and Langmuir probe measurements is provided. The model has also been employed to provide insight into recent ground test observations of the neutralizer cathode in NEXT. It is found that a likely cause of the observed keeper voltage drop is cathode orifice erosion. However, due to the small magnitude of this change, is approx. 0.5 V (less than 5% of the beginning-of-life value) over 10 khrs, and in light of the large uncertainties of the cathode material sputtering yield at low ion energies, other causes cannot be excluded. Preliminary simulations to understand transition to plume mode suggest that in the range of 3-5 sccm the existing 2-D model reproduces fairly well the rise of the keeper voltage in the NEXT neutralizer as observed in the laboratory. At lower flow rates the simulation produces oscillations in the keeper current and voltage that require prohibitively small time-steps to resolve with the existing algorithms.

electric propulsion↗

Calculating Free Energies Using Average Force

A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.

Darve, Eric↗