Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data intensive computing”

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 613 records · Page 34

Use of Apollo 17 Epoch Neutron Spectrum as a Benchmark in Testing LEND Collimated Sensor

The Apollo 17 neutron experiment LPNE provided a unique set of data on production of neutrons in the Lunar soil bombarded by Galactic Cosmic Rays (GCR). It serves as valuable "ground-truth" in the age of orbital remote sensing. We used the neutron data attributed to Apollo 17 epoch as a benchmark for testing the LEND's collimated sensor, as introduced by the geometry of collimator and efficiency of He3 counters. The latter is defined by the size of gas counter and pressure inside it. The intensity and energy spectrum of neutrons escaping the lunar surface are dependent on incident flux of Galactic Cosmic Rays (GCR) whose variability is associated with Solar Cycle and its peculiarities. We obtain first the share of neutrons entering through the field of view of collimator as a fraction of the total neutron flux by using the angular distribution of neutron exiting the Moon described by our Monte Carlo code. We computed next the count rate of the 3He sensor by using the neutron energy spectrum from McKinney et al. [JGR, 2006] and by consider geometry and gas pressure of the LEND sensor. Finally the neutron count rate obtained for the Apollo 17 epoch characterized by intermediate solar activity was adjusted to the LRO epoch characterized by low solar activity. It has been done by taking into account solar modulation potential, which affects the GCR flux, and in turn changes the neutron albedo flux.

Chin, Gordon↗

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 Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗

Velocity and Temperature Measurement in Supersonic Free Jets Using Spectrally Resolved Rayleigh Scattering

The flow fields of unheated, supersonic free jets from convergent and convergent-divergent nozzles operating at M = 0.99, 1.4, and 1.6 were measured using spectrally resolved Rayleigh scattering technique. The axial component of velocity and temperature data as well as density data obtained from a previous experiment are presented in a systematic way with the goal of producing a database useful for validating computational fluid dynamics codes. The Rayleigh scattering process from air molecules provides a fundamental means of measuring flow properties in a non-intrusive, particle free manner. In the spectrally resolved application, laser light scattered by the air molecules is collected and analyzed using a Fabry-Perot interferometer (FPI). The difference between the incident laser frequency and the peak of the Rayleigh spectrum provides a measure of gas velocity. The temperature is measured from the spectral broadening caused by the random thermal motion and density is measured from the total light intensity. The present point measurement technique uses a CW laser, a scanning FPI and photon counting electronics. The 1 mm long probe volume is moved from point to point to survey the flow fields. Additional arrangements were made to remove particles from the main as well as the entrained flow and to isolate FPI from the high sound and vibration levels produced by the supersonic jets. In general, velocity is measured within +/- 10 m/s accuracy and temperature within +/- 10 K accuracy.

Panda, J.↗

A Distributive, Non-Destructive, Real-Time Approach to Snowpack Monitoring

This invention is designed to ascertain the snow water equivalence (SWE) of snowpacks with better spatial and temporal resolutions than present techniques. The approach is ground-based, as opposed to some techniques that are air-based. In addition, the approach is compact, non-destructive, and can be communicated with remotely, and thus can be deployed in areas not possible with current methods. Presently there are two principal ground-based techniques for obtaining SWE measurements. The first is manual snow core measurements of the snowpack. This approach is labor-intensive, destructive, and has poor temporal resolution. The second approach is to deploy a large (e.g., 3x3 m) snowpillow, which requires significant infrastructure, is potentially hazardous [uses a approximately equal to 200-gallon (approximately equal to 760-L) antifreeze-filled bladder], and requires deployment in a large, flat area. High deployment costs necessitate few installations, thus yielding poor spatial resolution of data. Both approaches have limited usefulness in complex and/or avalanche-prone terrains. This approach is compact, non-destructive to the snowpack, provides high temporal resolution data, and due to potential low cost, can be deployed with high spatial resolution. The invention consists of three primary components: a robust wireless network and computing platform designed for harsh climates, new SWE sensing strategies, and algorithms for smart sampling, data logging, and SWE computation.

Frolik, Jeff↗

Direct simulation of high-speed mixing layers

A computational study of a nonreacting high-speed mixing layer is performed. A higher order algorithm with sufficient grid points is used to resolve all relevant scales. In all cases, a temporal free-stream disturbance is introduced. The resulting flow is time-sampled to generate a statistical cross section of the flow properties. The studies are conducted at two convective Mach numbers, three free-stream turbulence intensities, three Reynolds numbers, and two types of initial profiles-hyperbolic tangent (tanh) and boundary layer. The boundary-layer profile leads to more realistic predictions of the transition processes. The predicted transition Reynolds number of 0.18 x 10(exp 6) compares well with experimental data. Normalized vortex spacings for the boundary-layer case are about 3.5 and compare favorably with the 1.5 to 2.5 found in experimental measurements. The tanh profile produces spacings of about 10. The growth rate of the layer is shown to be moderately affected by the initial disturbance field, but comparison with experimental data shows moderate agreement. For the boundary-layer case, it is shown that noise at the Strouhal number of 0.007 is selectively amplified and shows little Reynolds number dependence.

Mukunda, H. S.↗

Computational Modeling of Space Physiology for Informing Spaceflight Countermeasure Design and Predictions of Efficacy

MOTIVATION: Spaceflight countermeasures mitigate the harmful effects of the space environment on astronaut health and performance. Exercise has historically been used as a countermeasure to physical deconditioning, and additional countermeasures including lower body negative pressure, blood flow occlusion and artificial gravity are being researched as countermeasures to spaceflight-induced fluid shifts. The NASA Digital Astronaut Project uses computational models of physiological systems to inform countermeasure design and to predict countermeasure efficacy.OVERVIEW: Computational modeling supports the development of the exercise devices that will be flown on NASAs new exploration crew vehicles. Biomechanical modeling is used to inform design requirements to ensure that exercises can be properly performed within the volume allocated for exercise and to determine whether the limited mass, volume and power requirements of the devices will affect biomechanical outcomes. Models of muscle atrophy and bone remodeling can predict device efficacy for protecting musculoskeletal health during long-duration missions. A lumped-parameter whole-body model of the fluids within the body, which includes the blood within the cardiovascular system, the cerebral spinal fluid, interstitial fluid and lymphatic system fluid, estimates compartmental changes in pressure and volume due to gravitational changes. These models simulate fluid shift countermeasure effects and predict the associated changes in tissue strain in areas of physiological interest to aid in predicting countermeasure effectiveness. SIGNIFICANCE: Development and testing of spaceflight countermeasure prototypes are resource-intensive efforts. Computational modeling can supplement this process by performing simulations that reduce the amount of necessary experimental testing. Outcomes of the simulations are often important for the definition of design requirements and the identification of factors essential in ensuring countermeasure efficacy.

biomedical data↗

Predictions of axisymmetric free turbulent shear flows using a generalized eddy-viscosity approach

The generalized eddy viscosity approach is described and results are presented of test cases which show that predictions obtained by this approach are adequate for most engineering applications. Because of the importance of starting computations from the injection station where experimentally determined mean and turbulence parameters are rarely available, a very simple core model applicable to simple step-type (slug) profiles was developed. Agreement between predicted and experimental mean profiles was generally almost as good for calculations made by using this model throughout the core region and the transition model for all subsequent regions as predictions made by starting from experimental profiles in the transition region. The generalized eddy-viscosity model, which was developed in part through correlation of turbulence parameters, successfully predicted turbulent shear stress, turbulent intensity, and mean velocity profiles for a 0.040-inch-diameter microjet. Therefore, successful scaling by the model was demonstrated since data used in its development was for jet areas up to 90,000 times as large as the microjet and velocities only 1/20th as high.

Morgenthaler, J. H.↗

Constructing an advanced software tool for planetary atmospheric modeling

Scientific model building can be an intensive and painstaking process, often involving the development of large and complex computer programs. Despite the effort involved, scientific models cannot be easily distributed and shared with other scientists. In general, implemented scientific models are complex, idiosyncratic, and difficult for anyone but the original scientist/programmer to understand. We believe that advanced software techniques can facilitate both the model building and model sharing process. In this paper, we describe a prototype for a scientific modeling software tool that serves as an aid to the scientist in developing and using models. This tool includes an interactive intelligent graphical interface, a high level domain specific modeling language, a library of physics equations and experimental datasets, and a suite of data display facilities. Our prototype has been developed in the domain of planetary atmospheric modeling, and is being used to construct models of Titan's atmosphere.

Keller, Richard M.↗

Vibrational Spectroscopy of Ions and Radicals Present in the Interstellar Medium and in Planetary Atmospheres: A Theoretical Study

Anharmonic vibrational frequencies and intensities are calculated for OH(H2O)n and H(H2O)n radicals (that form on icy particles of the interstellar medium), HCO radical (the main intermediate in the synthesis of organic molecules in space), NH2(-) and C2H(-) anions, H5(+) cation, and other systems relevant to interstellar chemistry. In addition to pure ions and radicals, their complexes with water are studied to assess the effects of water environment on infrared spectra. The calculations are performed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. Fundamental, overtone, and combination excitations are computed. The results are in good agreement with available experimental data and provide reliable predictions for vibrational excitations not yet measured in laboratory experiments. The data should be useful for interpretation of astronomically observed spectra and identification of ions and radicals present in the interstellar medium and in planetary atmospheres.

Chaban, Galina M.↗

Invasion of the Giant X-Ray Bubbles

More bizarre news from the bizarre world of neutron stars: They got dozens of giant X-ray bubbles taller than the Empire State Building forming and popping every second. Yet finding evidence of this bubble blowing was no childish endeavor. It took a team of Berkeley scientists just about all they could muster: three days' worth of observational data and 1,200 hours of number-crunching on the slickest of supercomputers at the Lawrence Livermore National Laboratory (LLNL). Mind you, these are the same powerful computers used to model nuclear reactions and calculate Bill Gates' salary. Neutron stars are already known as the bad boys of the cosmic neighborhood. They are the core remains of exploded stars, with a mass greater than our sun condensed into a sphere no bigger than Manhattan, A teaspoon of its densely-packed surface matter, the analogy goes, would weigh three billion tons back on earth. Its 10 trillion gauss magnetic field could wipe your credit card clean at a distance half-way to the moon. Its intense gravity lures gas from nearby stars, spinning the neutron star to speeds of up to a thousand times a second. But giant X-ray bubbles dancing on the surface? That even freaked out the discoverers. "Early on, I tended to disbelieve what the calculations were telling me," said Richard Klein, who splits his time between Berkeley and LLNL. "The bubbles 'appeared' and I wasn't convinced they were real."

Wanjek, Christopher↗

The Transition from Thick to Thin Plate Wake Physics: Whither Vortex Shedding?

The near and very near wake of a flat plate with a circular trailing edge is investigated with data from direct numerical simulations. Computations were performed for six different combinations of the Reynolds numbers based on plate thickness (D) and boundary layer momentum thickness upstream of the trailing edge (theta). Unlike the case of the cylinder, these Reynolds numbers are independent parameters for the flat plate. The separating boundary layers are turbulent in all the cases investigated. One objective of the study is to understand the changes in the wake vortex shedding process as the plate thickness is reduced (increasing theta/D). The value of D varies by a factor of 16 and that of theta by approximately 5 in the computations. Vortex shedding is vigorous in the low theta/D cases with a substantial decrease in shedding intensity in the large theta/D cases. Other shedding characteristics are also significantly altered with increasing theta/D. A visualization of the shedding process in the different cases is provided and discussed. The basic shedding mechanism is explored in depth. The effect of changing theta/D on the time-averaged, near-wake velocity statistics is also discussed. A functional relationship between the shedding frequency and the Reynolds numbers mentioned above is obtained.

CFD↗

Definition of the unsteady vortex flow over a wing/body configuration

A problem of current interest in computational aerodynamics is the prediction of unsteady vortex flows over aircraft at high angles of attack. A six-month experimental effort was conducted at the John H. Harper Wind Tunnel to acquire qualitative and quantitative information on the unsteady vortex flow over a generic wing-body configuration at high angles of attack. A double-delta flat-plate wing with beveled edges was combined with a slender sharp-nosed body-of-revolution fuselage to form the generic configuration. This configuration produces a strong attached leading edge vortex on the wing, as well as sharply-peaked flow velocity spectra above the wing. While it thus produces flows with several well-defined features of current interest, the model was designed for efficiency of representation in computational codes. A moderate number of surface pressure ports and two unsteady pressure sensors were used to study the pressure distribution over the wing and body surface at high angles of attack; the unsteady pressure sensing did not succeed because of inadequate signal-to-noise ratio. A pulsed copper vapor laser sheet was used to visualize the vortex flow over the model, and vortex trajectories, burst locations, mutual induction of vortex systems from the forebody, strake, and wing, were quantified. Laser Doppler velocimetry was used to quantify all 3 components of the time-average velocity in 3 data planes perpendicular to the freestream direction. Statistics of the instantaneous velocity were used to study intermittency and fluctuation intensity. Hot-film anemometry was used to study the fluctuation energy content in the velocity field, and the spectra of these fluctuations. In addition, a successful attempt was made to measure velocity spectra, component by component, using laser velocimetry, and these were compared with spectra measured by hot-film anemometry at several locations.

Liou, S. G.↗

Artificial intelligence support for scientific model-building

Scientific model-building can be a time-intensive and painstaking process, often involving the development of large and complex computer programs. Despite the effort involved, scientific models cannot easily be distributed and shared with other scientists. In general, implemented scientific models are complex, idiosyncratic, and difficult for anyone but the original scientific development team to understand. We believe that artificial intelligence techniques can facilitate both the model-building and model-sharing process. In this paper, we overview our effort to build a scientific modeling software tool that aids the scientist in developing and using models. This tool includes an interactive intelligent graphical interface, a high-level domain specific modeling language, a library of physics equations and experimental datasets, and a suite of data display facilities.

Keller, Richard M.↗

Sound Diffraction Modeling of Rotorcraft Noise Around Terrain

A new computational technique, Wave Confinement (WC), is extended here to account for sound diffraction around arbitrary terrain. While diffraction around elementary scattering objects, such as a knife edge, single slit, disc, sphere, etc. has been studied for several decades, realistic environments still pose significant problems. This new technique is first validated against Sommerfeld's classical problem of diffraction due to a knife edge. This is followed by comparisons with diffraction over three-dimensional smooth obstacles, such as a disc and Gaussian hill. Finally, comparisons with flight test acoustics data measured behind a hill are also shown. Comparison between experiment and Wave Confinement prediction demonstrates that a Poisson spot occurred behind the isolated hill, resulting in significantly increased sound intensity near the center of the shadowed region.

Stephenson, James H.↗

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

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

Qiuyu Xu↗

Isotopologue Consistency of Semi-Empirically Computed InfraRed Line Lists and Further Improvement for Rare Isotopologues: CO2 and SO2 Case Studies.

The semi-empirical molecular rovibrational IR line lists, such as ExoMol, TheoReTs, and Ames, combine the experimental accuracy and theoretical power to reach better than 0.1 cm-1 accuracy for line positions and better than 80-90% agreement for line intensities. The quality of these existing semi-empirical IR lists allows further improvements of intensity and line positions for those unobserved minor isotopologues. This paper presents our new BTRHE (Best Theory + Reliable High-resolution Experiment) strategy implementation. For line intensity, the isotopologue consistency and the patterns of mass dependence in the Ames-296K SO2 and CO2 IR lists are quantitatively presented along the mass-inverse coordinates. The consistency and patterns are better than those in existing experimental data. The methodology proposed here can be used to identify inconsistencies, outliers, and mistakes in intensities, and help improve Effective Dipole Model (EDM) and molecular IR databases. We call for an experimental study on the 50006 and 60007 bands of CO2 628. For line position predictions, a simple approach combining the variational IR line lists with Effective Hamiltonian (EH) model may refine the effective rotational constants A0/B0/C0 and quartic centrifugal distortion constants of minor isotopologues. The prediction accuracy may be improved by two orders of magnitude, i.e. reaching 0-5 MHz prediction accuracy in the range of J<20-30, Ka<10-20, and 0.01-0.02 MHz accuracy for A0/B0/C0. Several important factors have been systematically investigated and discussed, e.g. convergence, uncertainties, higher order terms, fixing EH parameters, mass coordinates, etc. Amicrowave (MW) line set consisting of 644,636 strong transitions for all 30 isotopologues and corresponding refined EH(Ames) parameters are reported in the supplementary material. This approach may be easily extended to rovibrational bands, hot bands, and other molecular systems.

Xinchuan Huang↗