Search NASA⌕ Search

SEARCH · Search NASA

Results for “point kinetics equations”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Integrated Power and Attitude Control for a Spacecraft with Flywheels and Control Moment Gyroscopes

A law is designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as all integrated set of actuators for attitude control. General. nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as it means of compensating for damping exerted by rotor bearings. Two flywheel 'steering laws' are developed such that torque commanded by all attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude and illustrate the benefits of kinetic energy error feedback.

Roithmayr, Carlos M.↗
NASA NTRS2003

Dynamics and Control of Attitude, Power, and Momentum for a Spacecraft Using Flywheels and Control Moment Gyroscopes

Several laws are designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as an integrated set of actuators for attitude control. General, nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as a means of compensating for damping exerted by rotor bearings. Two flywheel steering laws are developed such that torque commanded by an attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude, and illustrate the benefits of kinetic energy error feedback. Control laws for attitude hold are also developed, and used to show the amount of propellant that can be saved when flywheels assist the CMGs. Nonlinear control laws for large-angle slew maneuvers perform well, but excessive momentum is required to reorient a vehicle like the International Space Station.

Roithmayr, Carlos M.↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
NASA NTRS1974

A study of the solar wind angular momentum including proton thermal anisotropy

The solution to the steady state magnetohydrodynamic equations governing the supersonic expansion of the solar corona into interplanetary space is obtained for various assumptions regarding the form in which proton thermal energy is carried away from the sun. The one-fluid, inviscid, formulation of the MHD equations is considered assuming that thermal energy is carried away by conduction from a heat source located at the base of the corona. Angular motion of the solar wind led to the existence of three critical points through which the numerical solutions must pass to extend from the sun's surface to large heliocentric distances. The results show that the amount of magnetic field energy converted into kinetic energy in the solar wind is only a small fraction of the total expansion energy flux and has little effect upon the final radial expansion velocity.

Acuna, M. H.↗
NASA NTRS1995

Test particle study of minor ions in solar wind turbulence

We perform a parameter study of the temporal evolution of a test particle distribution function in MHD turbulence. The turbulent fields are calculated using a pseudo-spectral method and periodic boundary conditions on a regular grid of 180(exp 3) points, appropriate for incompressible, homogeneous and isotropic turbulence. Initially, the kinetic and the magnetic energy are equal on the average. Both, deterministic and random initial conditions are used, in the former case with zeros of the magnetic field located at grid points, in the latter case located by interpolation between grid points. The evolution of the minor ion distribution function is studied in detail as these turbulent fields evolve, developing strong current and vorticity sheets. Using the full collisionless equation of motion for the test particles, the efficiency of nonlinear interactions can be studied. The results are compared to theoretical predictions and are then discussed in connection with the observations of the dynamical properties of solar wind minor ions derived from in situ observations.

Zurbuchen, Th.↗
DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2013

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗
NASA NTRS2010

An Attempt to Derive the epsilon Equation from a Two-Point Closure

The goal of this paper is to derive the equation for the turbulence dissipation rate epsilon for a shear-driven flow. In 1961, Davydov used a one-point closure model to derive the epsilon equation from first principles but the final result contained undetermined terms and thus lacked predictive power. Both in 1987 and in 2001, attempts were made to derive the epsilon equation from first principles using a two-point closure, but their methods relied on a phenomenological assumption. The standard practice has thus been to employ a heuristic form of the equation that contains three empirical ingredients: two constants, c(sub 1 epsilon), and c(sub 2 epsilon), and a diffusion term D(sub epsilon) In this work, a two-point closure is employed, yielding the following results: 1) the empirical constants get replaced by c(sub 1), c(sub 2), which are now functions of Kappa and epsilon; 2) c(sub 1) and c(sub 2) are not independent because a general relation between the two that are valid for any Kappa and epsilon are derived; 3) c(sub 1), c(sub 2) become constant with values close to the empirical values c(sub 1 epsilon), c(sub epsilon 2), (i.e., homogenous flows); and 4) the empirical form of the diffusion term D(sub epsilon) is no longer needed because it gets substituted by the Kappa-epsilon dependence of c(sub 1), c(sub 2), which plays the role of the diffusion, together with the diffusion of the turbulent kinetic energy D(sub Kappa), which now enters the new equation (i.e., inhomogeneous flows). Thus, the three empirical ingredients c(sub 1 epsilon), c(sub epsilon 2), D (sub epsilon)are replaced by a single function c(sub 1)(Kappa, epsilon ) or c(sub 2)(Kappa, epsilon ), plus a D(sub Kappa)term. Three tests of the new equation for epsilon are presented: one concerning channel flow and two concerning the shear-driven planetary boundary layer (PBL).

Canuto, V. M.↗
NASA NTRS1976

On the physical interpretation of the cosmic-ray transport equations

It is pointed out that an earlier paper (Jokipii and Parker, 1967) had discussed in detail the question of the interpretation of the energy-change terms in the general equation for cosmic-ray transport in the solar wind, so that no need exists to reassess the physical interpretation, as urged by Gleeson and Webb (1974). The question concerns the expression for the adiabatic deceleration, which is 'replaced' by an expression for the mean rate of change of momentum of fast particles. It is pointed out that these two expressions represent two essentially different quantities. The adiabatic deceleration refers to the mean rate loss of kinetic energy of an individual particle as a result of the expansion of the gas and magnetic field of the wind. This parameter can in no way depend on the gradient of the fast particle density, as does the mean rate of change of momentum of the fast particles.

Jokipii, J. R.↗
NASA NTRS1992

Navier-Stokes solutions with surface catalysis for Martian atmospheric entry

In this study numerical solutions have been obtained for two-dimensional axisymmetric hypersonic nonequilibrium CO2 flow over a high angle blunt cone with appropriate surface boundary conditions to account for energy and mass conservation at the body surface. The flowfield is described by the Navier-Stokes equations and multicomponent conservation laws which account for both translational and internal vibrational nonequilibrium effects. Complete forebody solutions have been obtained for the peak heating point of the Mars entry trajectory specified in the proposed NASA MESUR (Mars Environmental Survey) project. In these solutions, radiative equilibrium wall temperature and surface heating distributions are determined over the MESUR aeroshell forebody for entry velocity equal to 7 km/sec with varying degrees of surface catalysis. The effects of gas kinetics, surface catalysis, transport properties, and vibrational relaxation times on the surface heating are examined. The results identify some important issues in the prediction of surface heating for flows in thermochemical nonequilibrium and show that the Navier-Stokes code used herein is effective for thermal protection system design and materials selection.

Chen, Y.-K.↗
NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, r d a U production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembe1 (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and platelike), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and in the mid-latitude continent with different concentrations of CCN: a low "c1ean"concentration and a high "dirty" concentration. In addition, differences and similarities between bulk microphysics and spectral-bin microphysical schemes will be examined and discussed.

Tao, Wei-Kuo↗
NASA NTRS2003

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitable affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distribution parameterized as spectral bin microphysics are needed to explicitly study the effect of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembel (GCE) model. The formulation for the explicit spectral-bim microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), groupel and frozen drops/hall] Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region and in the mid-latitude using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Besides the initial differences in aerosol concentration, preliminary results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CCN case but has less cloud water mass aloft. Because the spectral-bim model explicitly calculates and allows for the examination of both the mass and number concentration of cpecies in each size category, a detailed analysis of the instantaneous size spectrum can be obtained for the two cases. It is shown that since the low CCN case produces fever droplets, larger size develop due to greater condencational and collectional growth, leading to a broader size spectrum in comparison to the high CCN case.

Tao, Wei-Kuo↗
NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembel (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and in the mid-latitude continent with different concentrations of CCN: a low "c1ean"concentration and a high "dirty" concentration. In addition, differences and similarities between bulk microphysics and spectral-bin microphysical schemes will be examined and discussed.

Tao, Wei-Kuo↗
NASA NTRS1988

Correlations of velocity and temperature fluctuations in the stagnation-point flow of circular cylinder in turbulent flow

The present analyses of boundary layer flow and turbulence transport attempt to characterize the influence of freestream turbulence on the surface heat-transfer rate and stagnation point region skin friction of a circular cross-section cylinder in turbulent flow. The Reynolds stress-transport equations and k-epsilon two-equation turbulence modeling are used, yielding time-averaged turbulence double-correlations, mean-flow properties, surface heat-transfer rate, and skin-friction with freestream isotropic turbulence. A comparison of analytical results with experimental data indicates that large Reynolds normal stresses are induced at the boundary layer edge by the kinetic energy of the turbulence.

Wang, C. R.↗
NASA NTRS1990

A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

Brush, L. N.↗
NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles (i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail). Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region, in the sub-tropics (Florida) and in the mid-latitude using identical thermodynamic conditions but with different concentrations of CCN: a low 'clean' concentration and a high 'dirty' concentration.

Tao, Wei-Kuo↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e.,pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.A spectral-bin microphysical model is very expensive from a from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Besides the initial differences in aerosol concentration, preliminary results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CCN case but has less cloud water mass aloft. Because the spectral-bin model explicitly calculates and allows for the examination of both the mass and number concentration of species in each size categor, a detailed analysis of the instantaneous size spectrum can be obtained for the two cases. It is shown that since the low CCN case produces fewer droplets, larger sized develop due to the greater condensational and collectional growth, leading to a broader size spectrum in comparison to the high CCN case.

Tao, Wei-Kuo↗