Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JuMP model”

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 145 records · Page 8

Flood basalt eruptions, comet showers, and mass extinction events

A chronology of initiation dates of the major continental flood basalt episodes has been established from compilation of published K-Ar and Ar-Ar ages of basaltic flows and related basic intrusions. The dating is therefore independent of the biostratigraphic and paleomagnetic time scales, and the estimated errors of the inititation dates are approximately + or - 4 pct. There are 11 distinct episodes of continental flood basalts known during the past 250 Myr. The data show that flood basalt episodes are generally relatively brief geologic events, with intermittent eruptions during peak output periods lasting ony 2 to 3 Myr or less. Statistical analyses suggest that these episodes may have occurred quasi-periodically with a mean cycle time of 32 + or - 1 Myr. The initiation dates of the flood basalts are close to the estimated dates of marine mass extinctions and impact-crater clusters. Although a purely internal forcing might be argued for the flood basalt volcanism, quasi-periodic comet impacts may be the trigger for both the flood basalts and the extinctions. Impact cratering models suggest that large-body impactors lead to deep initial cratering, and therefore may cause mantle disturbances and initiate mantle plume activity. The flood basalt episodes commonly mark the initiation or jump of a mantle hotspot, and are often followed by continental rifting and separation. Evidence from dynamical studies of impacts, occurrences of craters and hotspots, and the geochemistry of boundary layers is synthesized to provide a possible model of impact-generated volcanism. Flood basalt eruptions may themselves have severe effects on climate, and possibly on life. Impacts might, as a result, have led to mass extinctions through direct atmospheric disturbances, and/or indirectly through prolonged flood basalt volcanism.

Rampino, Michael R.↗

Short Wavelength Electromagnetic Perturbations Excited Near the Solar Probe Plus Spacecraft in the Inner Heliosphere: 2.5D Hybrid Modeling

A 2.5D numerical plasma model of the interaction of the solar wind (SW) with the Solar Probe Plus spacecraft (SPPSC) is presented. These results should be interpreted as a basic plasma model derived from the SW-interaction with the spacecraft (SC), which could have consequences for both plasma wave and electron plasma measurements on board the SC in the inner heliosphere. Compression waves and electric field jumps with amplitudes of about 1.5 V/m and (12-18) V/m were also observed. A strong polarization electric field was also observed in the wing of the plasma wake. However, 2.5D hybrid modeling did not show excitation of whistler/Alfven waves in the upstream connected with the bidirectional current closure that was observed in short-time 3D modeling SPPSC and near a tether in the ionosphere. The observed strong electromagnetic perturbations may be a crucial point in the electromagnetic measurements planned for the future Solar Probe Plus (SPP) mission. The results of modeling electromagnetic field perturbations in the SW due to shot noise in absence of SPPSC are also discussed.

Lipatov, Alexander S.↗

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗

A stochastic model for eye movements during fixation on a stationary target.

A stochastic model describing small eye movements occurring during steady fixation on a stationary target is presented. Based on eye movement data for steady gaze, the model has a hierarchical structure; the principal level represents the random motion of the image point within a local area of fixation, while the higher level mimics the jump processes involved in transitions from one local area to another. Target image motion within a local area is described by a Langevin-like stochastic differential equation taking into consideration the microsaccadic jumps pictured as being due to point processes and the high frequency muscle tremor, represented as a white noise. The transform of the probability density function for local area motion is obtained, leading to explicit expressions for their means and moments. Evaluation of these moments based on the model is comparable with experimental results.

Vasudevan, R.↗

MHD Simulations of Earth's Bow Shock at low Mach Numbers: Standoff Distances

Global, three-dimensional, ideal MHD simulations of Earth's bow shock are reported for low Alfven Mach numbers M(sub A) and quasi-perpendicular magnetic field orientations. The simulations use a hard, infinitely conducting magnetopauause obstacle, with axisymmetric three-dimensional location given by scaled standard model, to directly address previous gasdynamic (GD) and field-aligned MHD (FA-MHD) work. Tests of the simulated shocks' density jumps X for 1.4 approx. less than MA approx. less than 10 and the high M(sub A) shock location, and reproduction of the GD relation between magnetosheath thickness and X for quasi-gasdynamic MHD runs with M(sub A) much greater than M(sub s), confirm that the MHD code is working correctly. The MHD simulations show the standoff distance a(sub s), increasing monotonically with decreasing M(sub A). Significantly larger a(sub s), are found at low M(sub A) than predicted by GD and phenomenological MHD models and FA-MHD simulations, as required qualitatively by observations. The GD and FA-MHD predictions err qualitatively, predicting either constant or decreasing a(sub s), with decreasing M(sub A). This qualitative difference between quasi- perpendicular MHD and FA-MHD simulations is direct evidence for a(sub s), depending on the magnetic field orientation Theta. The enhancement factor over the phenomenological MHD predictions at MA approx. 2.4 agrees quantitatively with one observatiorial estimate. A linear relationship is found between the magnetosheath thickness and X, modified both quantitatively and intrinsically by MHD effects from the GD result. The MHD and GD results agree in the high M(sub A) limit. An MHD theory is developed for a(sub s), restricted to sufficiently perpendicular Theta and high sonic Mach numbers M(sub s). It explains the simulation results with excellent accuracy. Observational and further simulation testing of this MHD theory, and of its predicted M(sub A), Theta, and M(sub s) effects, is desirable.

Cairns, Iver H.↗

A program for identification of linear systems

A program has been written for the identification of parameters in certain linear systems. These systems appear in biomedical problems, particularly in compartmental models of pharmacokinetics. The method presented here assumes that some of the state variables are regularly modified by jump conditions. This simulates administration of drugs following some prescribed drug regime. Parameters are identified by a least-square fit of the linear differential system to a set of experimental observations. The method is especially suited when the interval of observation of the system is very long.

Buell, J.↗

On the structure of contact binaries. I - The contact discontinuity

The problem of the interior structure of contact binaries is reviewed, and a simple resolution of the difficulties which plague the theory is suggested. It is proposed that contact binaries contain a contact discontinuity between the lower surface of the common envelope and the Roche lobe of the cooler star. This discontinuity is maintained against thermal diffusion by fluid flow, and the transition layer is thin to the extent that the dynamical time scale is short in comparison with the thermal time scale. The idealization that the transition layer has infinitesimal thickness allows a simple formulation of the structure equations which are closed by appropriate jump conditions across the discontinuity. The further imposition of the standard boundary conditions suffices to define a unique model for the system once the chemical composition, the masses of the two stars, and the orbital separation are specified.

Shu, F. H.↗

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

Altitude transitions in energy climbs

The aircraft energy-climb trajectory for configurations with a sharp transonic drag rise is well known to possess two branches in the altitude/Mach-number plane. Transition in altitude between the two branches occurs instantaneously, a 'corner' in the minimum-time solution obtained with the energy-state model. If the initial and final values of altitude do not lie on the energy-climb trajectory, then additional jumps (crude approximations to dives and zooms) are required at the initial and terminal points. With a singular-perturbation approach, a 'boundary-layer' correction is obtained for each altitude jump, the transonic jump being a so-called 'internal' boundary layer, different in character from the initial and terminal layers. The determination of this internal boundary layer is examined and some computational results for an example presented.

Weston, A. R.↗

Effect of mid-altitude ion heating on ion outflow at polar latitudes

The effect of ion heating on polar ion outflow, when either the parallel or perpendicular (or both) ion temperatures at the exobase are elevated above values typical for the ionosphere, was investigated using a modified semikinetic steady-state model of Barakat and Schunk (1983) that allowed for anisotropic ion heating at the exobase and eliminated lower boundary potential jumps. In addition, the relative impact of the ion heating vs electron heating on the oxygen escape fluxes was investigated by examining the flux of O(+) ions for various combinations of electron and ion temperatures. It is demonstrated that the O(+) escape flux can be increased, to levels as high as were obtained by Barakat and Schunk (1983) with the electron temperatures of 10,000 K, by raising, instead, the ion temperatures (but to values considerably less than the 100,000 K observed by Moore et al., 1986).

Li, Peng↗

Contact-impact analysis of deployable space systems

This paper presents simulation results of a flexible, deployable space system with multiple contact-impact events, using three different modeling approaches. In the first two approaches, the equations of motion of the system are written in terms of generalized coordinates of each component, with modal coordinates used to account for flexibility in the components. Two different methods are used to model contact-impact events in the systems. In the momentum balance method, as contact-impact occurs, momentum equations of the systems are solved for the jump discontinuity in velocity of generalized coordinates. Integration of the system is restarted with the new initial conditions from solutions of momentum equations. In the logical function method, a virtual spring connecting the two contact points is assumed. The virtual spring is active only when impact occurs. No special treatment is needed in the integration algorithm using this method. In the third approach, instead of using modal coordinates, a finite element based formulation for the transient dynamics of mechanical system is used to model flexibility of the components. The logical function method is used to model contact-impact in the third approach.

Wu, Shih-Chin↗

Lunar and Planetary Science XXXV: Origin of Planetary Systems

The session "Origin of Planetary Systems" included the following reports: (12753) Povenmire - Standard Comparison Small Main Belt Asteroid?; Gravitational Frequencies of Extra-Solar Planets; 'Jumping Jupiters' in Binary Star Systems; Hermes, Asteroid 2002 SY50 and the Northern Cetids - No Link Found!; What Kind of Accretion Model is Required for the Solar System; and Use of an Orbital Phase Curve of Extrasolar Planet for Specification of its Mass.

Source record↗

Can TOPEX/Poseidon Data be Used for Monitoring the Indonesian Throughflow?

The effects of the Indonesian throughflow(ITF) are studied in a reduced gravity, primitive equation, sigma coordinate model. The surface heat fluxes are provided by coupling the ocean GCAF to an advective atmospheric mixed layer model. There is an ENSO related signal in the ITF but the correlation with the Southern Oscillation Index (SOI) is only -0.31. When the winds over the Indian Ocean are held to climatology, this correlation jumps to -0.65 indicating that the non-EASO signal in the ITF is caused by the downstream winds. On interannual time-scales the ITF can be explained in terms of sea level differences between the western Pacific and eastern Indian Oceans when appropriate representative locations are chosen as demonstrated in both model and TOPEX data. It is shown that the main climatological effect of the ITF is to warm the Indian Ocean and cool the Pacific. Mile the main ENSO indices, NIN03 and NIN04, are almost identical with and without the ITF, the total SSTs show significant ENSO dependence.

Murtugudde, Ragu↗

HPTAM, a two-dimensional Heat Pipe Transient Analysis Model, including the startup from a frozen state

A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.

Tournier, Jean-Michel↗

Application of the algebraic RNG model for transition simulation

The algebraic form of the RNG model of Yakhot and Orszag (1986) is investigated as a transition model for the Reynolds averaged boundary layer equations. It is found that the cubic equation for the eddy viscosity contains both a jump discontinuity and one spurious root. A yet unpublished transformation to a quartic equation is shown to remove the numerical difficulties associated with the discontinuity, but only at the expense of merging both the physical and spurious root of the cubic. Jumps between the branches of the resulting multiple-valued solution are found to lead to oscillations in flat plate transition calculations. Aside from the oscillations, the transition behavior is qualitatively correct.

Lund, Thomas S.↗

Astrophysical implications and observational prospects of X-ray polarimetry

X-ray polarimetry is a prime tool for investigating the physics of compact objects, which has not been adequately exploited thus far. However, current low-cost technology and modest launch requirements could provide a large number of positive observations with a sensitivity factor at least 100 greater than 10 years ago.The amount of astrophysical information potentially to be gained from this is enormous. The introduction of polarimetric information (direction and degree) would bring a quantum jump in the parameter space used to investigate compact objects, from the current two (spectra and time variability) to four independent parameters that models need to satisfy. This should greatly improve our ability to discriminate between various possible models. Such observations could lead to an elucidation of the rotation-powered and accretion-powered pulsar radiation mechanisms, could help clinch the identification of black hole canditates, and could decide between thermal and nonthermal AGN radiation models, as well as pin down the geometry of the accretion flows in both galactic and extragalactic sources.

Meszaros, P.↗

Computation of turbulent flow in a thin liquid layer of fluid involving a hydraulic jump

Numerically computed flow fields and free surface height distributions are presented for the flow of a thin layer of liquid adjacent to a solid horizontal surface that encounters a hydraulic jump. Two kinds of flow configurations are considered: two-dimensional plane flow and axisymmetric radial flow. The computations used a boundary-fitted moving grid method with a k-epsilon model for the closure of turbulence. The free surface height was determined by an optimization procedure which minimized the error in the pressure distribution on the free surface. It was also checked against an approximate procedure involving integration of the governing equations and use of the MacCormack predictor-corrector method. The computed film height also compared reasonably well with previous experiments. A region of recirculating flow was found to be present adjacent to the solid boundary near the location of the jump, which was caused by a rapid deceleration of the flow.

Rahman, M. M.↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗