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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nonequilibrium information entropy approach to ternary fission of actinides

Ternary fission of actinides probes the state of the nucleus at scission. Light clusters are produced in space and time very close to the scission point. Within the nonequilibrium statistical operator method, a generalized Gibbs distribution is constructed from the information given by the observed yields of isotopes. Using this relevant statistical operator, yields are calculated taking excited states and continuum correlations into account, in accordance with the virial expansion of the equation of state. Furthermore, clusters with mass number A≤10 are well described using the nonequilibrium generalizations of temperature and chemical potentials. Improving the virial expansion, in-medium effects may become of importance in determining the contribution of weakly bound states and continuum correlations to the intrinsic partition function. Yields of larger clusters, which fail to reach this quasiequilibrium form of the relevant distribution, are described by nucleation kinetics, and a saddle-to-scission relaxation time of about 7000 fm/c is inferred. Light-charged particle emission, described by reaction kinetics and virial expansions, may therefore be regarded as a very important tool to probe the nonequilibrium time evolution of actinide nuclei during fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Optimal aeroassisted coplanar orbital transfer using an energy model

The atmospheric portion of the trajectories for the aeroassisted coplanar orbit transfer was investigated. The equations of motion for the problem are expressed using reduced order model and total vehicle energy, kinetic plus potential, as the independent variable rather than time. The order reduction is achieved analytically without an approximation of the vehicle dynamics. In this model, the problem of coplanar orbit transfer is seen as one in which a given amount of energy must be transferred from the vehicle to the atmosphere during the trajectory without overheating the vehicle. An optimal control problem is posed where a linear combination of the integrated square of the heating rate and the vehicle drag is the cost function to be minimized. The necessary conditions for optimality are obtained. These result in a 4th order two-point-boundary-value problem. A parametric study of the optimal guidance trajectory in which the proportion of the heating rate term versus the drag varies is made. Simulations of the guidance trajectories are presented.

Halyo, Nesim↗
NASA NTRS2003

Termolecular ion-molecule reactions in Titan's atmosphere. IV. A search made at up to 1 micron in pure hydrocarbons

The results of a study of ion-molecule reactions occurring in pure methane, acetylene, ethylene, ethane, propyne, propene, propane, and diacetylene at pressures up to 40 microns of pressure are reported. A variety of experimental methods are used: The standard double resonance in an ICR, for determination of the precursor ions and the modulated double resonance ejection in an ICR, for the determination of the daughter ions. The FA-SIFT technique was used for validation and examination of termolecular reactions with rate coefficients that are less than 10(-26) cm(6) s(-1). An extensive database of reaction kinetics already exists for many of these reactions. The main point of this study was the determination of the accuracy of this database and to search for any missing reactions and reaction channels that may have been omitted from earlier investigations. A specific objective of this work was to extend the study to the highest pressures possible to find out if there were any important termolecular reaction channels occurring. A new approach was used here. In the pure hydrocarbon gases the mass spectra were followed as a function of the pressure changes of the gas. An initial guess was first made using the current literature as a source of the reaction kinetics that were expected. A model of the ion abundances was produced from the solution of the partial differential equations in terms of reaction rate coefficients and initial abundances. The experimental data was fitted to the model for all of the pressures by a least squares minimization to the reaction rate coefficients and initial abundances. The reaction rate coefficients obtained from the model were then compared to the literature values. Several new channels and reactions were discovered when the modeled fits were compared to the actual data. This is all explained in the text and the implications of these results are discussed for the Titan atmosphere.

Hydrocarbons, Acyclic/analysis↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Proposed Common Model of Multi-phase Strength and Equation of State for a Tri-laboratory Collaboration (Working Draft 1.1)

The proposed simple common model for multiphase strength and EoS (CMMP) is meant to be sufficiently simple that each of the collaborating labs can share in a common starting point. Another objective is to start with relatively simple assumptions, which will not necessarily capture details of the physical processes, and incrementally add complexity in order to identify the minimal-needed technical detail. Through this co-evolution of model and experiment, we will better learn the importance of various theoretical approximations and where to invest future resources in experiment and model development. This simple framework is based on pressure and temperature equilibrium of all co-existing phases combined with deviatoric stress averaging through a volume fraction weighted ow stress and a volume fraction weighted shear modulus. Implementations of the framework based on equilibrium phase fractions (i.e. instantaneous kinetic rate) and for finite rate transition kinetics are proposed.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Flight test results for the Daedalus and Light Eagle human powered aircraft

The results of the flight test program of the Daedalus and Light Eagle human powered aircraft in the winter of 1987/88 are given. The results from experiments exploring the Light Eagle's rigid body and structural dynamics are presented. The interactions of these dynamics with the autopilot design are investigated. Estimates of the power required to fly the Daedalus aircraft are detailed. The system of sensors, signal conditioning boards, and data acquisition equipment used to record the flight data is also described. In order to investigate the dynamics of the aircraft, flight test maneuvers were developed to yield maximum data quality from the point of view of estimating lateral and longitudinal stability derivatives. From this data, structural flexibility and unsteady aerodynamics have been modeled in an ad hoc manner and are used to augment the equations of motion with flexibility effects. Results of maneuvers that were flown are compared with the predictions from the flexibility model. To extend the ad hoc flexibility model, a fully flexible aeroelastic model has been developed. The model is unusual in the approximate equality of many structural natural frequencies and the importance of unsteady aerodynamic effects. the Gossamer Albatross. It is hypothesized that this inverse ground effect is caused by turbulence in the Earth's boundary layer. The diameters of the largest boundary layer eddies (which represent most of the turbulent kinetic energy) are proportional to altitude; thus, closer to the ground, the energy in the boundary layer becomes concentrated in eddies of smaller and smaller diameter. Eventually the eddies become sufficiently small (approximately 0.5 cm) that they trip the laminar boundary layer on the wing. As a result, a greater percentage of the wing area is covered with turbulent flow. Consequently the aircraft's drag and the pow er required both increase as the aircraft flies closer to the ground. The results of the flight test program are examined critically, especially with respect to future applications. Maneuvers that allow the observation of stability derivatives for flexible aircraft are detailed. Considerations for the design of autopilots for future human power aircraft and high-altitude RPV's are discussed. Techniques useful in estimating the power required to fly aircraft of very high lift to drag ratio are described.

Sullivan, R. Bryan↗
NASA NTRS1988

Turbulence modeling in hypersonic inlets

A study is conducted to analyze the performance of different turbulence models when applied to flow through a Mach 7.4 hypersonic inlet. The analysis, which is two-dimensional, is done by comparing computational results from a Parabolized Navier-Stokes code and a full Navier-Stokes code, with experimental data. The McDonald-Camarata (MC) and Baldwin-Lomax (BL) models were the two zero-equation models used in the study. The Turbulent Kinetic Energy (TKE) model was chosen as a representative higher order model. The MC model, when run with transition of flow, provides a solution which compares excellently with the data. Transition has a first order effect on the overall solution provided by the code. The BL model predicts separation of flow in the inlet, which contradicts experimental findings. The TKE model does not perform any better than the MC and BL models, despite the fact that it is a higher order turbulence model. The BL and TKE models predict transition in the inlet at a location which is much earlier than observed in the experiment. This may be attributed to the empirical constants used to determine the point of transition.

Ng, W. F.↗
NASA NTRS1986

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗
NASA NTRS2015

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗
NASA NTRS1985

Geocoronal structure - The effects of solar radiation pressure and the plasmasphere interaction

The theory of planetary exospheres is extended to incorporate solar radiation pressure in a rigorous manner, and an evaporative geocoronal prototype (classical, motionless exobase) is constructed using Liouville's theorem. Model calculations for density and kinetic temperature at points along the earth-sun axis (solar and antisolar directions) reveal an extensive satellite component, comprising approximately 2/3 of the total hydrogen density near 10 earth radii, and a temperature profile suggestive of an isotropic quasi-Maxwellian velocity distribution for the bound component. A geotail is also evident as an enhancement of the density at local midnight compared to local noon that increases outward (from approximately 25 percent at 10 earth radii to over 60 percent at 20 earth radii). Additional mechanisms acting upon the geocorona alter the basic evaporative case in notable ways. Solar ionization has been included in a simple fashion; the effect is to partially deplete the density without otherwise altering the structure. Interaction with a simple plasmasphere via the Boltzmann equation results in 'heating' the geocorona and enhancing the escape flux at the expense of the density of the bound component, an effect not appreciated in earlier studies; the geotail survives this interaction.

Bishop, J.↗
NASA NTRS1991

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗
DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis under Boosted Spark-ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel's performance irrespective of its chemical composition, was numerically investigated In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel's knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (= 98) and MON (= 90), were first generated using a non-linear regression model based on artificial neural network (ANN). Three unconventional bioblendstocks were included in the analysis: Di-isobutylene (DIB), Isobutanol and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed was tabulated for each blend as a function of pressure, temperature and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON and engine operating conditions, the TPRF-Anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (DI) expectation was shown to be caused by differences in laminar flame speed (LFS). However; it was found that relatively large fuel-specific differences in LFS (> 20%) would have to be present to cause any appreciable deviation from the 01 framework. Otherwise, RON and MON would still be robust enough to predict a fuel's knock-limited performance.

Pal, Pinaki↗
NASA NTRS2001

Thermal/Pyrolysis Gas Flow Analysis of Carbon Phenolic Material

Provided in this study are predicted in-depth temperature and pyrolysis gas pressure distributions for carbon phenolic materials that are externally heated with a laser source. Governing equations, numerical techniques and comparisons to measured temperature data are also presented. Surface thermochemical conditions were determined using the Aerotherm Chemical Equilibrium (ACE) program. Surface heating simulation used facility calibrated radiative and convective flux levels. Temperatures and pyrolysis gas pressures are predicted using an upgraded form of the SINDA/CMA program that was developed by NASA during the Solid Propulsion Integrity Program (SPIP). Multispecie mass balance, tracking of condensable vapors, high heat rate kinetics, real gas compressibility and reduced mixture viscosity's have been added to the algorithm. In general, surface and in-depth temperature comparisons are very good. Specie partial pressures calculations show that a saturated water-vapor mixture is the main contributor to peak in-depth total pressure. Further, for most of the cases studied, the water-vapor mixture is driven near the critical point and is believed to significantly increase the local heat capacity of the composite material. This phenomenon if not accounted for in analysis models may lead to an over prediction in temperature response in charring regions of the material.

Clayton, J. Louie↗
NASA NTRS2004

Update on the Status of the FLUKA Monte Carlo Transport Code

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. Here we review the progresses achieved in the last year on the physics models. From the point of view of hadronic physics, most of the effort is still in the field of nucleus--nucleus interactions. The currently available version of FLUKA already includes the internal capability to simulate inelastic nuclear interactions beginning with lab kinetic energies of 100 MeV/A up the the highest accessible energies by means of the DPMJET-II.5 event generator to handle the interactions for greater than 5 GeV/A and rQMD for energies below that. The new developments concern, at high energy, the embedding of the DPMJET-III generator, which represent a major change with respect to the DPMJET-II structure. This will also allow to achieve a better consistency between the nucleus-nucleus section with the original FLUKA model for hadron-nucleus collisions. Work is also in progress to implement a third event generator model based on the Master Boltzmann Equation approach, in order to extend the energy capability from 100 MeV/A down to the threshold for these reactions. In addition to these extended physics capabilities, structural changes to the programs input and scoring capabilities are continually being upgraded. In particular we want to mention the upgrades in the geometry packages, now capable of reaching higher levels of abstraction. Work is also proceeding to provide direct import into ROOT of the FLUKA output files for analysis and to deploy a user-friendly GUI input interface.

Pinsky, L.↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2002

Solutal Convection Around Growing Protein Crystal and Diffusional Purification in Space

This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities. Onset of convection was considered analytically and numerically. Crystal growth was characterized by slow surface incorporation kinetics, i.e. growth kinetic coefficient beta (cm/s) small as compared to the typical bulk diffusion rate, D(sub 1)/h, where D(sub 1) is diffusivity of major crystallizing protein and h is the crystal size. Scaling type analysis predicted two laws on how the convection rate, v, essentially the Peclet number, Pe exactly equal to vh/D(sub 1), depends on dimensionless kinetic coefficient a exactly equal to beta h/D(sub 1). Namely: Pe = C(sub 2/5)(aRa(sup 2/5)) and Pe = C(sub 1) aRa. Here, Reynolds number Ra = rho(sub 1)(sup 0)gh(sup 3)(rho(sub p) - rho(sub w))/rho(sup p)rho(sub 1)vD(sub 1), v being solution viscosity. The constants C(sub 2/5), exactly equal to 0.28 and C(sub 1) exactly equal to 10(exp -2) found from the full scale computer simulation for a cylindrical crystal inside big cylindrical vessel. The linear boundary conditions connecting protein and impurity concentration at the interface with the flux to/from the interface was applied. No-slip condition for Navier-Shocker equations was employed. With these conditions, flow and concentration distributions were calculated. Validity of the Pe(Ra) dependencies follows for wide range of parameters for which numerical calculations have been accomplished and presented by various points.

Chernov, A. A.↗