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 271 records · Page 15

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

DyFK Simulation of Field-Aligned Ion Flows Observed by POLAR within Convecting Flux Over the Polar Ionosphere

Ion (O+ and H+) parallel flows along antisunward convecting flux tubes across the polar ionosphere from day to night side are simulated by an extended Dynamic Fluid semiKinetic (DyFK) model. The collision dominated portion of the flux tubes is treated with a moment-based fluid model for altitudes from 120 ?1100 km, while the generalized semikinetic model is used for the topside through 3 RE region. The effects of cleft/auroral soft electron precipitation and wave-driven transverse ion heating are incorporated into the generalized semi-kinetic treatment of topside ionosphere. The simulated evolution of field-aligned ion flow parameters is compared with observations made by the Thermal Ion Dynamics Experiment (TIDE) on board the POLAR satellite near 5000 km altitude over the southern hemisphere polar ionosphere.

Tu, J.-N.↗

A model for the kinetics of a solar-pumped long path laser experiment

A kinetic model for a solar-simulator pumped iodine laser system is developed and compared to an experiment in which the solar simulator output is dispersed over a large active volume (150 cu cm) with low simulator light intensity (approx. 200 solar constants). A trace foreign gas which quenches the upper level is introduced into the model. Furthermore, a constant representing optical absorption of the stimulated emission is introduced, in addition to a constant representing the scattering at each of the mirrors, via the optical cavity time constant. The non-uniform heating of the gas is treated as well as the pressure change as a function of time within the cavity. With these new phenomena introduced into the kinetic model, a best reasonable fit to the experimental data is found by adjusting the reaction rate coefficients within the range of known uncertainty by numerical methods giving a new bound within this range of uncertainty. The experimental parameters modeled are the lasing time, laser pulse energy, and time to laser threshold.

Stock, L. V.↗

Equilibrium and kinetic aspects of brittle fracture

The lattice statics method for simulating a brittle crack will be employed to determine: (1) the energy-based criterion for crack stability, (2) at a given stress intensity in the stable region, the kinetic barrier for crack motion, and (3) the assumptions that break down in the Griffith theory when generalized to the lattice model. A simple model is used in the analysis, but this model is general enough to provide a conceptual framework for more realistic models.

Esterling, D. M.↗

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↗

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

shock wave/boundary-layer interaction↗

A study of the kinetic energy generation with general circulation models

The history data of winter simulation by the GLAS climate model and the NCAR community climate model are used to examine the generation of atmospheric kinetic energy. The contrast between the geographic distributions of the generation of kinetic energy and divergence of kinetic energy flux shows that kinetic energy is generated in the upstream side of jets, transported to the downstream side and destroyed there. The contributions from the time-mean and transient modes to the counterbalance between generation of kinetic energy and divergence of kinetic energy flux are also investigated. It is observed that the kinetic energy generated by the time-mean mode is essentially redistributed by the time-mean flow, while that generated by the transient flow is mainly responsible for the maintenance of the kinetic energy of the entire atmospheric flow.

Chen, T.-C.↗

Zoledronate Prevents Simulated Weightlessness-Induced Bone Loss in the Cancellous Compartment While Blunting the Efficacy of a Mechanical Loading Countermeasure

Astronauts using high-force resistance training while weightless show a high-turnover remodeling state within the skeletal system, with resorption and formation biomarkers being elevated. One countermeasure for the skeletal health of astronauts includes an antiresorptive of the bisphosphonate (BP) drug class. We asked, does the combination of an anti-resorptive and high-force exercise during weightlessness have negative effects on bone remodeling and strength? In this study, we developed an integrated model to mimic the mechanical strain of exercise via cyclical loading (CL) in mice treated with the BP Zoledronate (ZOL) combined with hind limb unloading (HU) to simulate weightlessness. We hypothesized that ZOL prevents structural degradation from simulated weightlessness and that CL and ZOL interact to render CL less effective. Thirty-two C57BL/6 mice (male, 16 weeks old, n=8/group) were exposed to 3 weeks of either HU or normal ambulation (NA). Cohorts of mice received one subcutaneous injection of ZOL (45μg/kg), or saline vehicle (VEH), prior to the start of HU. The right tibia was axially compressed in vivo 60x/day to 9N (+1200μstrain on the periosteal surface) and repeated 3x/week during HU. Left tibiae served as a within subject, non-compressed control. Ex vivo μCT was performed on all subjects to determine cancellous and cortical architectural parameters. Static and dynamic histomorphometry were carried out for the left and right tibiae to determine osteoclast- and osteoblast relevant surfaces. Further, micro damage was assessed in select groups by basic-fuchsin staining to test whether CL had an effect. For all assays, a multivariate (2x2x2) ANCOVA model was used to account for body weight changes. Additionally, for the tibiae, we incorporated a random effect for the subject (hence, a mixed model) to account for observations of both left and right tibiae within each subject. P < 0.05 was considered significant. In the cancellous compartment of the proximal tibial metaphysis, we observed a main effect from each independent variable, as determined by structural and histomorphometric assays. Specifically, as expected, ZOL showed an increase in the cancellous bone volume to total volume fraction (BV/TV, +32%) and trabecular number (+18%) compared to the VEH. As expected, ZOL decreased osteoclast surface (OC/BS) by -45% compared to VEH. Surprisingly, ZOL reduced mineralizing surface (MS/BS) and bone formation rate (BFR), indicators of osteoblast activity, by -40% and -54%, respectively, compared to VEH. Altogether, ZOL-treated mice displayed a low turnover state of remodeling in the metaphysis. In the context of skeletal aging, we speculate that ZOL prevented age-related cancellous strut loss during the experiment. As a main effect, as expected, HU decreased BV/TV by - 31% via reductions in both trabecular thickness (-11%) and number (-22%) compared to NA controls. Additionally, HU decreased MS/BS by -38% and bone formation rate (BFR) by -50% compared to NA controls. Altogether, these data are consistent with structural degradation resulting from imbalanced remodeling that favors resorption. As a main effect, CL increased BV/TV by +15% via increased trabecular thickness (+12%) compared to the noncompressed limb. As expected, CL increased MS/BS (+20%) and BFR (+24%), indicating osteoblast mineralization contributed to bone gains. These data show that CL provided an anabolic stimulus to the cancellous tissue. We observed unique interactions in ZOL*CL and HU*CL. First, ZOL prevented CL-induced increases in BV/TV and trabecular number, as compared to VEH. In the context of skeletal aging, these data suggest no added benefit from CL in the ZOL-treated mice. Interestingly, no microdamage was observed in mice that were unloaded and treated with ZOL (independent of CL). Secondly, HU prevented CL-induced increases in BFR, as compared to NA controls. These data suggest that either exercise is less effective or the kinetics of formation are slower during simulated weightlessness. Osteoclast surface was unchanged by either treatment. Thus, in contrast to exercising astronauts, these data do not suggest a high-turnover state in the metaphysis. To assess mechanical properties as a function of HU or ZOL, we tested the left femur in three-point bending ex vivo. As expected, HU decreased stiffness (-30%) compared to NA, and ZOL increased stiffness compared to VEH (+28%). Interestingly, HU increased the post-yield displacement, related to collagenous tensile loading, compared to NA (+20%). ZOL increased yield force (+11%) and ultimate force (+17%), which seems to explain the significant effect of ZOL increasing total energy (work-to-fracture, +15%), while not affecting the post yield displacement. Taken together, ZOL did not have detrimental affect on mechanical properties. Our integrated model simulates the combination of weightlessness, exercise-induced mechanical strain, and anti-resorptive treatment that astronauts experience during space missions. We conclude that Zoledronate was an effective countermeasure against weightlessness-induced bone loss, though zoledronate, as well as weightlessness, rendered exercise-related mechanical loading less effective.

bone↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 2: The two-dimensional model

The two-dimensional model of stratospheric constituents is presented in detail. The derivation of pertinent transport parameters and the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations that represent the chemical kinetic terms, and appropriate methods for simulating the diurnal variations of the solar zenith angle and species concentrations are discussed. Predicted distributions of tracer constituents (ozone, carbon 14, nitric acid) are compared with observed distributions.

Whitten, R. C.↗

Sensitivity of bone cell populations to weightlessness and simulated weightlessness

A rat suspension model for simulating certain aspects of weightlessness is discussed. Perturbations in physiological systems induced by this head down suspension model are verified by flight data. Findings of a suppression of osteoblast differentiation help explain the inhibition of bone formation inflight and during Earth-bound simulations. Since the anatomical site for these studies was in the maxilla, which is gravity loaded but non weightbearing in ground-based simulations, the similarity of bone cell kinetic changes, both inflight and in the ground-based model, suggest that fluid shifts rather than unloading may play an important role in bone alterations, at least at this sampling site.

Roberts, W. E.↗

A computer simulation of thin film nucleation and growth: The Volmer-Weber case

The computer simulation of thin film nucleation and growth, which was previously performed for the case of single monolayer, was modified to include multilayer growth via the Volmer-Weber mechanism. The simulation results show that: (1) the kinetics of multilayer film growth is nearly identical to that of monolayer growth; (2) when no reevaporation takes place, the cluster density resulting from multilayer growth is higher at high coverage than that resulting from monolayer growth; (3) when reevaporation does take place, the cluster density resulting from multilayer growth is nearly identical to that resulting from monolayer growth. This is not due, however, to similarity in microstructure.

Salik, J.↗

Reynolds-stress and dissipation rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗

Reynolds-stress and dissipation-rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗

Generation of strong MHD Alfvenic turbulence

Strong Alfvenic turbulence containing a number of solitonlike structures propagating at super-Alfvenic speeds is generated self-consistently and studied by means of computer simulation. A one-dimensional hybrid (kinetic ions, fluid electrons) code is used to investigate the nonlinear evolution of an electromagnetic ion-beam instability that generates low-frequency Alfven-like waves. As the instability develops, the field-aligned hydromagnetic waves steepen, forming a soliton that bifurcates several times, leading to a fully turbulent state.

Akimoto, K.↗

Extinction of a Free Methanol Droplet in Microgravity

The combustion of isolated, unsupported methanol droplets under microgravity conditions has been studied in the NASA-Lewis 2.2-sec drop tower. The substitution of O2/He mixtures for O2/N2 mixtures results in a substantial increase in the oxygen index required for ignition, the rate of gasification, and the diameter at extinction. As a result, droplets with initial size larger than 1 mm can be produced, deployed, ignited, and burned to extinction, all under microgravity conditions in this facility. A comprehensive time dependent numerical model incorporating detailed multicomponent molecular transport mechanisms and complex elementary chemical kinetics was developed and applied to simulate the entire transient droplet combustion processes (ignition/quasi-steady burning/extinction). Model outputs show good agreement with the experimental measurements, without parameter adjustment. Analysis of modeling results shows that the thermal conductivity of the gas phase under combustion environments increased by a factor of two when the 50-percent O2/50 percent He mixture is used instead of air. The increased thermal conductivity accelerates the heat transfer from the flame front to the droplet surface, resulting in a factor of two increase in the burning rate.

Cho, S. Y.↗

Multiple Satellite Observations of High-Latitude Ionospheric Outflows

We will present reported observations of up- and down-flows of topside ionospheric thermal plasmas from multiple near-simultaneous tracks through the high-latitude topside ionosphere. From several Southern polar passes, it has been possible to construct plots of field-aligned flows of O+ observed by the Thermal Ion Dynamics Experiment(TIDE) on the POLAR spacecraft near 5000 km altitude together with vertical ion flow observations from one or more DMSP spacecraft near 800 km altitude. We will also involve simulations from our Dynamic Fluid-Kinetic(DyFK) modeling of polar plasma transport and display resulting altitude profiles of ion parallel velocities and densities, and overlay those profiles with the "conjunction" measurements by DMSP(800 km) and POLAR(approx. 5000 km altitude). We also will present simultaneous observations of POLAR auroral UVI images with field-aligned flows.

Horwitz, J.↗