Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 397 records · Page 22

Numerical Simulation of AFP Nip Point Temperature Prediction for Complex Geometries

Material placement at the ideal nip point temperature over complex surfaces with uniformity across the width of the compaction rollers results in optimized part properties for Automated Fiber Placement (AFP) processes. However, current AFP systems utilize heat control models and methodologies, based on multiple process parameters such as feed-rate and orientation, that are mostly open-loop. Here, infrared (IR) heater input is calibrated as a function of process parameters during machine qualification. This work presents a numerical simulation to predict arrayed-infrared (AIR) emitter radiation onto a substrate that includes view factor implementation, IR radiative heat flow calculation, energy rate balance, and a transient heat transfer model. The purpose of this numerical model is to predict nip point temperature on complex surfaces, serving as a baseline for a new arrayed-infrared (AIR) thermoset heater to improve AFP process control. It is anticipated that this simulation will accurately control the temperature for high-speed AFP layup of complex geometries. An anticipated result of an AIR heater system is that material calibration and testing will be reduced as temperature is instantaneously monitored and controlled. Therefore, temperature across the roller width will be uniform during placement of complex parts, independent of their geometry.

Xia, Kaishu↗

Experiment and numerical simulation of RF heating in the Tandem Mirror plasma propulsion device

Recent results of experimental and theoretical studies of ICRF heating of the plasma in the Tendem Mirror rocket are presented. The radial and axial profile of the magnetic field of the wave have been measured, and the data agree with the results from numerical simulation. One very important new finding is that the wave damped when it approached the resonance plane. This is a strong indication of beach heating effect and that RF power is absorbed by the ions in the plasma as expected. This power absorption phenomenon was also further confirmed by the refined analytical study of the wave propagation in a slab model. The electron density in the central cell has been measured by a microwave interferometer. The existence of two types of discharge similar to those of H-alpha emission have also been observed from the density measurements. A preliminary design of a flight system for a 100-day Mars transit has been undertaken. The specific weight for a 10 MW rocket is only 0.04 kg/kW.

Yang, T. F.↗

A mesoscale numerical simulation of an index for use in the forecasting of the potential severity of wet microbursts

This paper demonstrates the feasibility of using the wet microburst potential intensity index developed by Proctor (1989) for predicting the onset of a wet microburst. After describing the index and the principles involved in its derivation, the version 4.0 of the Mesoscale Atmospheric Simulation System (MASS) model described by Karyampudi et al. (1988) was used to illustrate how the index would perform as a forecast product. An analysis is presented of the use of the microburst potential intensity index in the numerical simulation of the DFR aircraft accident on August 2, 1985 in a downburst environment.

Zack, John W.↗

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

High Ice Water Content in Tropical Cyclones during NASA/FAA Radar Flight Campaigns with Comparison to Numerical Simulations

Results from two recent flight campaigns into High Ice Water Content (HIWC) are summarized. HIWC was detected with airborne radar and measured via microphysical probes mounted on NASA’s DC-8. Results from the campaign demonstrate detectability of HIWC using airborne radar. In the cases examined, tropical storms are more proficient at producing large areas of HIWC than hurricanes. Three-dimensional, numerical simulations of tropical cyclones encountered during the campaign are also presented and are compared to measurements in order to understand conditions associated with HIWC. Within Hurricane Lane (2018), the highest ice water contents were found in regenerating thermal plumes located in the eye wall and in feeder bands. Within Tropical Storm Danny (2015), large areas of ice water content greater than 1 g m-3 were fed by broad-areas of convection occurring downwind from the center of circulation.

Aviation hazards↗

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗

Gas Core Reactor Numerical Simulation Using a Coupled MHD-MCNP Model

Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.

Kazeminezhad, F.↗

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗

Numerical simulation of the nonlinear response of composite plates under combined thermal and acoustic loading

A time-domain study of the random response of a laminated plate subjected to combined acoustic and thermal loads is carried out. The features of this problem also include given uniform static inplane forces. The formulation takes into consideration a possible initial imperfection in the flatness of the plate. High decibel sound pressure levels along with high thermal gradients across thickness drive the plate response into nonlinear regimes. This calls for the analysis to use von Karman large deflection strain-displacement relationships. A finite element model that combines the von Karman strains with the first-order shear deformation plate theory is developed. The development of the analytical model can accommodate an anisotropic composite laminate built up of uniformly thick layers of orthotropic, linearly elastic laminae. The global system of finite element equations is then reduced to a modal system of equations. Numerical simulation using a single-step algorithm in the time-domain is then carried out to solve for the modal coordinates. Nonlinear algebraic equations within each time-step are solved by the Newton-Raphson method. The random gaussian filtered white noise load is generated using Monte Carlo simulation. The acoustic pressure distribution over the plate is capable of accounting for a grazing incidence wavefront. Numerical results are presented to study a variety of cases.

Mei, Chuh↗

Impact of dilution on microbial community structure and functional potential: comparison of numerical simulations and batch culture experiments

A series of microcosm experiments was performed using serial dilutions of a sewage microbial community to inoculate a set of batch cultures in sterile sewage. After inoculation, the dilution-defined communities were allowed to regrow for several days and a number of community attributes were measured in the regrown assemblages. Based upon a set of numerical simulations, community structure was expected to differ along the dilution gradient; the greatest differences in structure were anticipated between the undiluted-low-dilution communities and the communities regrown from the very dilute (more than 10(-4)) inocula. Furthermore, some differences were expected among the lower-dilution treatments (e.g., between undiluted and 10(-1)) depending upon the evenness of the original community. In general, each of the procedures used to examine the experimental community structures separated the communities into at least two, often three, distinct groups. The groupings were consistent with the simulated dilution of a mixture of organisms with a very uneven distribution. Significant differences in community structure were detected with genetic (amplified fragment length polymorphism and terminal restriction fragment length polymorphism), physiological (community level physiological profiling), and culture-based (colony morphology on R2A agar) measurements. Along with differences in community structure, differences in community size (acridine orange direct counting), composition (ratio of sewage medium counts to R2A counts, monitoring of each colony morphology across the treatments), and metabolic redundancy (i.e., generalist versus specialist) were also observed, suggesting that the differences in structure and diversity of communities maintained in the same environment can be manifested as differences in community organization and function.

Non-NASA Center↗

High Ice Water Content in Tropical Cyclones During NASA/FAA Radar Flight Campaigns with Comparison to Numerical Simulations

Results from two recent flight campaigns into High Ice Water Content (HIWC) are summarized. HIWC was detected with airborne radar and measured via microphysical probes mounted on NASA’s DC-8. Results from the campaign demonstrate detectability of HIWC using airborne radar. In the cases examined, tropical storms are more proficient at producing large areas of HIWC than hurricanes. Three-dimensional, numerical simulations of tropical cyclones encountered during the campaign are also presented and are compared to measurements in order to understand conditions associated with HIWC. Within Hurricane Lane (2018), the highest ice water contents were found in regenerating thermal plumes located in the eye wall and in feeder bands. Within Tropical Storm Danny (2015), large areas of ice water content greater than 1 g/cu. m were fed by broad-areas of convection occurring downshear from the center of circulation.

Fred H. Proctor↗

Numerical simulation of a catastrophe model for coronal mass ejections

Van Tend and Kuperus (1978) proposed a simple catastrophe model for magnetically driving coronal mass ejections, prominence eruptions, and two-ribbon flares. Their model, which is based on simple circuit concepts, suggests that a stable configuration containing a current filament will lose equilibrium when the filament current exceeds a critical value. Here, a two-dimensional numerical simulation is used to test how the Van Tend-Kuperus model works in an ideal MHD fluid. The simulation exhibits the expected loss of mechanical equilibrium near the predicted critical value, but the current filament moves only a short distance upward before coming to rest at a new equilibrium. However, this new equilibrium contains a current sheet which is resistively unstable to magnetic reconnection, and if magnetic reconnection occurs rapidly, the filament can continue to move upward at Alfvenic speeds.

Forbes, T. G.↗

Study of turbulence structure through numerical simulations: The perspective of drag and reduction

The first part of this note concerns organized structures in the near-wall region of a turbulent boundary layer. A particular emphasis is given to those organized motions believed to be responsible for the turbulence production process in the wall region, and hence, indirectly responsible for the viscous drag. Examples are selected to illustrate how the analyses of numerically generated databases have contributed to improve our understanding of the organized structures of turbulent boundary layers. In the second part, results from an exploratory study based on the direct numerical simulation of the concepts for active control of turbulent boundary layers are presented. A significant drag reduction is achieved when the surface boundary condition is modified such that it could suppress the organized motion in the wall region. This drag reduction is accompanied by a significant reduction in the intensity of the organized structures and in the magnitude of Reynolds stresses throughout the flow.

Kim, J. J.↗

Numerical simulations of positively-biased probes and dielectric-conductor disks in a plasma

Plasma densities in the low-earth orbit range may be sufficient to cause difficulties for spacecraft operating at high voltages in this environment. The present investigation is concerned with the results of a continuing effort to develop a particle-in-cell (PIC), cylindrically-symmetric, 2-1/2 dimensional, self-consistent numerical simulation code. The simulation has the objective to explore the interactions of an ambient plasma with a conducting disk, which may be partially covered by a dielectric material. The disk and the surrounding dielectric material represent a hole in an insulator covering a conductor. Attention is given to a review of the simulation model, plain disk calculations, the disk and dielectric configuration, and the 'pinhole' effect.

Brandon, S. T.↗

Numerical simulation of filling a magnetic flux tube with a cold plasma: Anomalous plasma effects

Large-scale models of plasmaspheric refilling have revealed that during the early stage of the refilling counterstreaming ion beams are a common feature. However, the instability of such ion beams and its effect on refilling remain unexplored. In order to learn the basic effects of ion beam instabilities on refilling, we have performed numerical simulations of the refilling of an artificial magnetic flux tube. (The shape and size of the tube are assumed so that the essential features of the refilling problem are kept in the simulation and at the same time the small scale processes driven by the ion beams are sufficiently resolved.) We have also studied the effect of commonly found equatorially trapped warm and/or hot plasma on the filling of a flux tube with a cold plasma. Three types of simulation runs have been performed.

Singh, Nagendra↗

Visualization of numerical simulation results: An interactive tool for mission oriented theory

With the advent of multispacecraft missions, it became obvious that theoretical models are needed, not only to interpret local spacecraft observations, but also to link single point measurements. Models that address space plasma phenomena on local and global scales, that involve nonlinear processes, and that use self-consistent approaches, were developed. An interactive system for visualizing and analyzing numerical simulation results was also developed. This system is based on visualization modules which use the application visualization system and the NCAR graphic packages. Examples from recent simulations illustrate how these modules can be used for displaying and manipulating simulation results.

Berchem, Jean↗

Importance of inlet boundary conditions for numerical simulation of combustor flows

Fluid dynamic computer codes for the mathematical simulation of problems in gas turbine engine combustion systems are required as design and diagnostic tools. To eventually achieve a performance standard with these codes of more than qualitative accuracy it is desirable to use benchmark experiments for validation studies. Typical of the fluid dynamic computer codes being developed for combustor simulations is the TEACH (Teaching Elliptic Axisymmetric Characteristics Heuristically) solution procedure. It is difficult to find suitable experiments which satisfy the present definition of benchmark quality. For the majority of the available experiments there is a lack of information concerning the boundary conditions. A standard TEACH-type numerical technique is applied to a number of test-case experiments. It is found that numerical simulations of gas turbine combustor-relevant flows can be sensitive to the plane at which the calculations start and the spatial distributions of inlet quantities for swirling flows.

Sturgess, G. J.↗