Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗
NASA NTRS1974

Study of ceramic products and processing techniques in space

An analysis of the solidification kinetics of beta alumina in a zero-gravity environment was carried out, using computer-simulation techniques, in order to assess the feasibility of producing high-quality single crystals of this material in space. The two coupled transport processes included were movement of the solid-liquid interface and diffusion of sodium atoms in the melt. Results of the simulation indicate that appreciable crystal-growth rates can be attained in space. Considerations were also made of the advantages offered by high-quality single crystals of beta alumina for use as a solid electrolyte; these clearly indicate that space-grown materials are superior in many respects to analogous terrestrially-grown crystals. Likewise, economic considerations, based on the rapidly expanding technological applications for beta alumina and related fast ionic conductors, reveal that the many superior qualities of space-grown material justify the added expense and experimental detail associated with space processing.

Markworth, A. J.↗
NASA NTRS1977

The effect of ambient conditions on carbon monoxide emissions from an idling gas turbine combustor

A test program employing a gas turbine combustor is outlined; the results of which quantize the effects of changes in ambient temperature and humidity on carbon monoxide emissions at simulated idle operating conditions. A comparison of the experimental results with analytical results generated by a kinetic model of the combustion process, and reflecting changing ambient conditions, is given. It is demonstrated that for a complete range of possible ambient variations, significant changes do occur in the amount of carbon monoxide emitted by a gas turbine at idle, and that the analytical model is reasonably successful in predicting changes.

Subramanian, A. K.↗
NASA NTRS1981

Strong Langmuir turbulence in one and two dimensions

A comparison of one- and two-dimensional simulations of the strong turbulence evolution of Lagmuir waves such as would be generated by the kinetic beam-plasma instability is reported. It is shown that the initial formation of the localized fields is one dimensional.

Rowland, H. L.↗
NASA NTRS1995

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗
NASA NTRS2004

Montmorillonite, oligonucleotides, RNA and origin of life

Na-montmorillonite prepared from Volclay by the titration method facilitates the self-condensation of ImpA, the 5'-phosphorimidazolide derivative of adenosine. As was shown by AE-HPLC analysis and selective enzymatic hydrolysis of products, oligo(A)s formed in this reaction are 10 monomer units long and contain 67% 3',5'-phosphodiester bonds (Ferris and Ertem, 1992a). Under the same reaction conditions, 5'-phosphorimidazolide derivatives of cytidine, uridine and guanosine also undergo self-condensation producing oligomers containing up to 12-14 monomer units for oligo(C)s to 6 monomer units for oligo(G)s. In oligo(C)s and oligo(U)s, 75-80% of the monomers are linked by 2',5'-phosphodiester bonds. Hexamer and higher oligomers isolated from synthetic oligo(C)s formed by montmorillonite catalysis, which contain both 3',5'- and 2',5'-linkages, serve as catalysts for the non-enzymatic template directed synthesis of oligo(G)s from activated monomer 2-MeImpG, guanosine 5'-phospho-2-methylimidazolide (Ertem and Ferris, 1996). Pentamer and higher oligomers containing exclusively 2',5'-linkages, which were isolated from the synthetic oligo(C)s, also serve as templates and produce oligo(G)s with both 2',5'- and 3',5'-phosphodiester bonds. Kinetic studies on montmorillonite catalyzed elongation rates of oligomers using the computer program SIMFIT demonstrated that the rate constants for the formation of oligo(A)s increased in the order of 2-mer < 3-mer < 4-mer ... < 7-mer (Kawamura and Ferris, 1994). A decameric primer, dA(pdA)8pA bound to montmorillonite was elongated to contain up to 50 monomer units by daily addition of activated monomer ImpA to the reaction mixture (Ferris, Hill and Orgel, 1996). Analysis of dimer fractions formed in the montmorillonite catalyzed reaction of binary and quaternary mixtures of ImpA, ImpC, 2-MeImpG and ImpU suggested that only a limited number of oligomers could have formed on the primitive Earth rather than equal amounts of all possible isomers (Ertem and Ferris, 2000). Formation of phosphodiester bonds between mononucleotides by montmorillonite catalysis is a fascinating discovery, and a significant step forward in efforts to find out how the first RNA-like oligomers might have formed in the course of chemical evolution. However, as has been pointed out in several publications, these systems should be regarded as models rather than a literal representation of prebiotic chemistry (Orgel, 1998; Joyce and Orgel, 1999; Schwartz, 1999).

RNA/chemical synthesis↗
NASA NTRS2000

A New Eddy Dissipation Rate Formulation for the Terminal Area PBL Prediction System(TAPPS)

The TAPPS employs the MASS model to produce mesoscale atmospheric simulations in support of the Wake Vortex project at Dallas Fort-Worth International Airport (DFW). A post-processing scheme uses the simulated three-dimensional atmospheric characteristics in the planetary boundary layer (PBL) to calculate the turbulence quantities most important to the dissipation of vortices: turbulent kinetic energy and eddy dissipation rate. TAPPS will ultimately be employed to enhance terminal area productivity by providing weather forecasts for the Aircraft Vortex Spacing System (AVOSS). The post-processing scheme utilizes experimental data and similarity theory to determine the turbulence quantities from the simulated horizontal wind field and stability characteristics of the atmosphere. Characteristic PBL quantities important to these calculations are determined based on formulations from the Blackadar PBL parameterization, which is regularly employed in the MASS model to account for PBL processes in mesoscale simulations. The TAPPS forecasts are verified against high-resolution observations of the horizontal winds at DFW. Statistical assessments of the error in the wind forecasts suggest that TAPPS captures the essential features of the horizontal winds with considerable skill. Additionally, the turbulence quantities produced by the post-processor are shown to compare favorably with corresponding tower observations.

Charney, Joseph J.↗
NASA NTRS1988

Kinetic isolation tether experiment

Progress was made on the analysis of tether damping and on experimentation of the control system on the laboratory simulator. The damping analysis considers the dynamics of a long tether connecting two spacecraft in Earth orbit, one of the spacecraft having dominant mass. In particular, it considers the material damping of the tether. The results show that, with properly chosen tether material and braiding structure, longitudinal vibration of the tethered system is well damped. A particularly effective method of implementing attitude control for tethered satellites is to use the tether tension force to generate control torques by moving the tether attach point relative to the satellite center of mass. A scaled, one dimensional laboratory simulation of the KITE mission was built and preliminary experiments of the proposed attitude control system were performed. The simulator was built to verify theoretical predictions of attitude controllability, and to investigate the technological requirements in order to implement this concept. A detailed description of the laboratory apparatus is provided, and in addition, the results of the preliminary experiments are presented and discussed.

Powell, J. David↗
NASA NTRS2016

Optical and Near-Infrared Observations of SN 2013DX Associated with GRB 130702A

We present optical and near-infrared (NIR) light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be E(sub gamma, iso) = 6.4(+1.3/-1.0) x 10(exp 50) erg (1 keV to 10 MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed g'r'i'z' light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves approx. 20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined SNe Ic, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of approx. 21,000 km/s. We construct a quasi-bolometric (g'r'z'yJ) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond (Delta)t > 40 days. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a Ni-56 mass of M(sub Ni) = 0.37+/- 0.01 Stellar Mass, an ejecta mass of M(sub ej) = 3.1+/- 0.1 Stellar Mass, and a kinetic energy of E(sub K) = (8.2+/- 0.43) x 10(exp 51) erg (statistical uncertainties only), consistent with previous GRB-associated supernovae. When considering the ensemble population of GRB-associated supernovae, we find no correlation between the mass of synthesized Ni-56 and high-energy properties, despite clear predictions from numerical simulations that M(sub Ni) should correlate with the degree of asymmetry. On the other hand, M(sub Ni) clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

gamma-ray burst: individual (GRB 130702A) – supe↗
NASA NTRS2019

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗
NASA NTRS1996

Kinetic Space Weather: Toward a Global Hybrid Model of the Polar Ionosphere-Lower Magnetosphere Plasma Transport

During the indicated period of performance, we had a number of publications concerned with kinetic polar ionosphere-lower magnetosphere plasma transport. For the IUGG 1991-4 Quadrennial Report, we reviewed aspects of U.S. accomplishments concerned with polar plasma transport, among other issues. In another review, we examined the computer simulations of multiple-scale processes in space plasmas, including polar plasma outflow and transport. We also examined specifically multiscale processes in ionospheric outflows. We developed a Generalized Semi-Kinetic(GSK) model for the topside-lower magnetosphere which explored the synergistic action of wave heating and electric potentials in the formation of auroral Ion conics, in particular the "pressure cooker" mechanism. We extended the GSK model all the way down to 120 km and applied this code to illustrate the response of the ionosphere- magnetosphere to soft-electron precipitation and convection-driven frictional ion heating, respectively. Later, the convection-driven heating work was extended to a paper for the Journal of Geophysical Research. In addition to the above full published papers, we also presented the first developments of the coupled fluid-semikinetic model for polar plasma transport during this period. The results from a steady-state treatment were presented, with the second presentation being concerned with the effects of photo-electrons on the polar wind, and the first garnering an outstanding student paper award from the American Geophysical Union. We presented the first results from a time-dependent version of this coupled fluid-semikinetic model.

Horwitz, James L.↗
NASA NTRS1981

Simulation of the planetary boundary layer with the UCLA general circulation model

A planetary boundary layer (PBL) model is presented which employs a mixed layer entrainment formulation to describe the mass exchange between the mixed layer with the upper, laminar atmosphere. A modified coordinate system couples the mixed layer model with large scale and sub-grid scale processes of a general circulation model. The vertical coordinate is configured as a sigma coordinate with the lower boundary, the top of the PBL, and the prescribed pressure level near the tropopause expressed as coordinate surfaces. The entrainment mass flux is parameterized by assuming the dissipation rate of turbulent kinetic energy to be proportional to the positive part of the generation by convection or mechanical production. The results of a simulation of July are presented for the entire globe.

Suarez, M. J.↗
NASA NTRS1990

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗
NASA NTRS2024

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗
NASA NTRS2002

Simulation of Energy Response of the ATIC Calorimeter

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. For reconstruction of primary spectra from spectra of energy deposits measured in the experiment, correlations between kinetic energy of a primary particle E(sub kin) and energy deposit in the calorimeter E(sub d) should be known. For this purpose, simulations of energy response of the calorimeter on energy spectra of different nuclei were done. The simulations were performed by GEANT-3.21 code with QGSM generator for nucleus - nucleus interactions. The incident flux was taken as isotropic in the ATIC aperture. Primary spectra power-law by momentum were used as inputs according to standard models of cosmic ray acceleration. These spectra become power-law by kinetic energy at E(sub kin) higher than approx.20Mc(sup 2), where M is primary nucleus mass. It should be noted that energy deposit spectra measured by ATIC illustrate similar behavior. Distributions of ratio E(sub kin)/E(sub d) are presented for different energy deposits and for a set of primaries. For power-law regions of energy spectra at E(sub d)> or equal to 20Mc(sup 2) the obtained mean value of E(sub kin)/E(sub d) increases from approx.2.4 for protons to approx.3.1 for iron, while rms/ decreases from 50% for protons to about 15% for iron. These values were obtained for the spectral index gamma=1.6

Batkov, K. E.↗
NASA NTRS1990

Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

Goussis, D. A.↗
NASA NTRS1991

Direct simulation of compressible turbulence in a shear flow

Compressibility effects on the turbulence in homogeneous shear flow are investigated. The growth of the turbulent kinetic energy was found to decrease with increasing Mach number: a phenomenon which is similar to the reduction of turbulent velocity intensities observed in experiments on supersonic free shear layers. An examination of the turbulent energy budget shows that both the compressible dissipation and the pressure-dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation is predominantly negative in homogeneous shear flow, in contrast to its predominantly positive behavior in isotropic turbulence. The different signs of the pressure-dilatation are explained by theoretical consideration of the equations for the pressure variance and density variance. Previously, the following results were obtained for isotropic turbulence: (1) the normalized compressible dissipation is of O(M(sub t)(exp 2)); and (2) there is approximate equipartition between the kinetic and potential energies associated with the fluctuating compressible mode. Both of these results were substantiated in the case of homogeneous shear. The dilatation field is significantly more skewed and intermittent than the vorticity field. Strong compressions seem to be more likely than strong expansions.

Sarkar, S.↗
NASA NTRS2004

Computation Molecular Kinetics Model of HZE Induced Cell Cycle Arrest

Cell culture models play an important role in understanding the biological effectiveness of space radiation. High energy and charge (HZE) ions produce prolonged cell cycle arrests at the G1/S and G2/M transition points in the cell cycle. A detailed description of these phenomena is needed to integrate knowledge of the expression of DNA damage in surviving cells, including the determination of relative effectiveness factors between different types of radiation that produce differential types of DNA damage and arrest durations. We have developed a hierarchical kinetics model that tracks the distribution of cells in various cell phase compartments (early G1, late G1, S, G2, and M), however with transition rates that are controlled by rate-limiting steps in the kinetics of cyclin-cdk's interactions with their families of transcription factors and inhibitor molecules. The coupling of damaged DNA molecules to the downstream cyclin-cdk inhibitors is achieved through a description of the DNA-PK and ATM signaling pathways. For HZE irradiations we describe preliminary results, which introduce simulation of the stochastic nature of the number of direct particle traversals per cell in the modulation of cyclin-cdk and cell cycle population kinetics. Comparison of the model to data for fibroblast cells irradiated photons or HZE ions are described.

Cucinotta, Francis A.↗