Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spreading models”

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 793 records · Page 44

A non-PRE double-peaked burst with oscillations: burning front propagation and stalling

Non-photospheric-radius-expansion (non-PRE) double-peaked bursts may be explained in terms of spreading (and temporary stalling) of thermonuclear flames from a rotational pole on the neutron star surface, as we argued in a previous study. Here we analyze Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) data of such a burst from the low mass X-ray binary (LMXB) system 4U 1636-536, and show that our model (with ignition at high latitudes) can qualitatively explain the observed burst profile, and spectral evolution. Moreover, the evolution of the source radius inferred from the data shows a strong signature of temporary stalling of the burning front, which is an essential ingredient of our model. This implies that an understanding of thermonuclear flame spreading on neutron stars can be achieved by a simultaneous study of the evolution of intensity and spectrum of these bursts. We also report the discovery of millisecond period brightness oscillations from this burst, which is the first such observation from a non-PRE double-peaked burst. Our model can explain the corresponding oscillation amplitude during the first (weaker) peak, and the absence of oscillations during the second peak. We discuss how observations of oscillations during non-PRE double-peaked bursts provide an additional t 001 for understanding thermonuclear flame spreading successfully.

Bhattacharyya, Sudip↗

TLife-LSTM: Forecasting Future COVID-19 Progression with Topological Signatures of Atmospheric Conditions

Understanding the impact of atmospheric conditions on SARS-CoV2 is critical to model COVID-19 dynamics and sheds a light on the future spread around the world. Furthermore, geographic distri- butions of expected clinical severity of COVID-19 may be closely linked to prior history of respiratory diseases and changes in humidity, tem- perature, and air quality. In this context, we postulate that by tracking topological features of atmospheric conditions over time, we can provide a quanti?able structural distribution of atmospheric changes that are likely to be related to COVID-19 dynamics. As such, we apply the machinery of persistence homology on time series of graphs to extract topological signatures and to follow geographical changes in relative humidity and temperature. We develop an integrative machine learning framework named Topological Lifespan LSTM (TLife-LSTM) and test its predictive capabilities on forecasting the dynamics of SARS-CoV2 cases. We validate our framework using the number of con?rmed cases and hospitalization rates recorded in the states of Washington and California in the USA. Our results demonstrate the predictive potential of TLife-LSTM in forecasting the dynamics of COVID-19 and modeling its complex spatio-temporal spread dynamics.

Gel, Yulia R.↗

Reflight of the Solid Surface Combustion Experiment: Flame Radiation Near Extinction

In flame spreading in quiescent and low-velocity opposing flows, effects of surface reradiation and flame radiation are important in establishing the spread rate and whether the flame, once ignited, survives to steady spread or extinguishes after a time long compared to the ignition event. A reflight of the Solid Surface Combustion Experiment (SSCE), supported by modelling, demonstrates that for thick, flat fuels, the ultimate fate of the flame is extinction rather than steady spread. A mismatch between the thermal scale in the gas, driven by radiation, and the species diffusion scale, driven by mass diffusion, develops such that the high temperature regions of the flame are ultimately located in a region to which oxygen cannot be supplied at a sufficient rate to sustain reaction, and extinction occurs. Results of the experiment conducted on Space Shuttle mission STS 85 on 9 August 1997 are reviewed. For the flat surface geometry, while the hydrodynamic phenomena associated with opposed-flow flame spread may be treated two dimensionally, the radiative effects are three dimensional, and so modelling the radiative processes, with the mismatch in dimensionality, is difficult. The cylindrical geometry at least one long compared to the radius, provides a configuration in which the radiative processes for spread in the axial direction are two dimensional, thus simplifying the modelling. The cylindrical geometry allows for the development of more sophisticated radiative models without the complication of dimensionality concerns, e.g., discrete transfer, which is discussed in detail by Bundy (1998). Additionally, the cylindrical geometry results in a "focussing" of the heat transfer to the surface and may allow for steady spread for radii that for thick fuels of the same half-thickness there is no steady spread.

Altenkirch, R. A.↗

Spread Across Liquids: The World's First Microgravity Combustion Experiment on a Sounding Rocket

The Spread Across Liquids (SAL) experiment characterizes how flames spread over liquid pools in a low-gravity environment in comparison to test data at Earth's gravity and with numerical models. The modeling and experimental data provide a more complete understanding of flame spread, an area of textbook interest, and add to our knowledge about on-orbit and Earthbound fire behavior and fire hazards. The experiment was performed on a sounding rocket to obtain the necessary microgravity period. Such crewless sounding rockets provide a comparatively inexpensive means to fly very complex, and potentially hazardous, experiments and perform reflights at a very low additional cost. SAL was the first sounding-rocket-based, microgravity combustion experiment in the world. It was expected that gravity would affect ignition susceptibility and flame spread through buoyant convection in both the liquid pool and the gas above the pool. Prior to these sounding rocket tests, however, it was not clear whether the fuel would ignite readily and whether a flame would be sustained in microgravity. It also was not clear whether the flame spread rate would be faster or slower than in Earth's gravity.

Source record↗

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied here. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗

Opposed-Flow Flame Spreading in Reduced Gravity

Experimental results obtained in drop towers and in Space Shuttle based experiments coupled with modelling efforts are beginning to provide information that is allowing an understanding to be developed of the physics of opposed-flow flame spread at reduced gravity where the spread rate and flow velocity are comparable and of the role played by radiative and diffusive processes in flame spreading in microgravity. Here we describe one Space Shuttle based experiment on flame spreading in a quiescent environment, the Solid Surface Combustion Experiment, SSCE, one planned microgravity experiment on flame spreading in a radiatively-controlled, forced opposing flow environment, the Diffusive and Radiative Transport in Fires Experiment, DARTFire, modelling efforts to support these experiments, and some results obtained to date.

Altenkirch, Robert A.↗

Integrable Many-Body Quantum Floquet-Thouless Pumps

We construct an interacting integrable Floquet model featuring quasiparticle excitations with topologically nontrivial chiral dispersion. This model is a fully quantum generalization of an integrable classical cellular automaton. Here, we write down and solve the Bethe equations for the generalized quantum model and show that these take on a particularly simple form that allows for an exact solution: essentially, the quasiparticles behave like interacting hard rods. The generalized thermodynamics and hydrodynamics of this model follow directly, providing an exact description of interacting chiral particles in the thermodynamic limit. Although the model is interacting, its unusually simple structure allows us to construct operators that spread with no butterfly effect; this construction does not seem possible in other interacting integrable systems. This model exemplifies a new class of exactly solvable, interacting quantum systems specific to the Floquet setting.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Contributions of Different Cloud Types to Feedbacks and Rapid Adjustments in CMIP5*

When using five climate model simulations of the response to an abrupt quadrupling of CO 2 , the authors perform the first simultaneous model intercomparison of cloud feedbacks and rapid radiative adjustments with cloud masking effects removed, partitioned among changes in cloud types and gross cloud properties. After CO 2 quadrupling, clouds exhibit a rapid reduction in fractional coverage, cloud-top pressure, and optical depth, with each contributing equally to a 1.1 W m -2 net cloud radiative adjustment, primarily from shortwave radiation. Rapid reductions in midlevel clouds and optically thick clouds are important in reducing planetary albedo in every model. As the planet warms, clouds become fewer, higher, and thicker, and global mean net cloud feedback is positive in all but one model and results primarily from increased trapping of longwave radiation. As was true for earlier models, high cloud changes are the largest contributor to intermodel spread in longwave and shortwave cloud feedbacks, but low cloud changes are the largest contributor to the mean and spread in net cloud feedback. The importance of the negative optical depth feedback relative to the amount feedback at high latitudes is even more marked than in earlier models. Furthermore, the authors show that the negative longwave cloud adjustment inferred in previous studies is primarily caused by a 1.3 W m -2 cloud masking of CO 2 forcing. Properly accounting for cloud masking increases net cloud feedback by 0.3 W m -2 K -1 , whereas accounting for rapid adjustments reduces by 0.14 W m -2 K -1 the ensemble mean net cloud feedback through a combination of smaller positive cloud amount and altitude feedbacks and larger negative optical depth feedbacks.

58 GEOSCIENCES↗

Evaluation of the Simulation of Arctic and Antarctic Sea Ice Coverages by Eleven Major Global Climate Models

Comparison of polar sea ice results from 11 major global climate models and satellite-derived observations for 1979-2004 reveals that each of the models is simulating seasonal cycles that are phased at least approximately correctly in both hemispheres. Each is also simulating various key aspects of the observed ice cover distributions, such as winter ice not only throughout the central Arctic basin but also throughout Hudson Bay, despite its relatively low latitudes. However, some of the models simulate too much ice, others too little ice (in some cases varying depending on hemisphere and/or season), and some match the observations better in one season versus another. Several models do noticeably better in the Northern Hemisphere than in the Southern Hemisphere, and one does noticeably better in the Southern Hemisphere. In the Northern Hemisphere all simulate monthly average ice extents to within +/-5.1 x 10(exp 6)sq km of the observed ice extent throughout the year; and in the Southern Hemisphere all except one simulate the monthly averages to within +/-6.3 x 10(exp 6) sq km of the observed values. All the models properly simulate a lack of winter ice to the west of Norway; however, most do not obtain as much absence of ice immediately north of Norway as the observations show, suggesting an under simulation of the North Atlantic Current. The spread in monthly averaged ice extents amongst the 11 model simulations is greater in the Southern Hemisphere than in the Northern Hemisphere and greatest in the Southern Hemisphere winter and spring.

Parksinson, Claire↗

Vortex ring bubbles

Toroidal bubbles with circulation are studied numerically and by means of a physically motivated model equation. Two series of computations are performed by a boundary-integral method. One set shows the starting motion of an initially spherical bubble as a gravitationally driven liquid jet penetrates through the bubble from below causing a toroidal geometry to develop. The jet becomes broader as surface tension increases and fails to penetrate if surface tension is too large. The dimensionless circulation that develops is not very dependent on the surface tension. The second series of computations starts from a toroidal geometry, with circulation determined from the earlier series, and follows the motion of the rising and spreading vortex ring. Some modifications to the boundary-integral formulation were devised to handle the multiply connected geometry. The computations uncovered some unexpected rapid oscillations of the ring radius. These oscillations and the spreading of the ring are explained by the model equation which provides a more general description of vortex ring bubbles than previously available.

Lundgren, T. S.↗

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

Agent-Based Models for COVID-19 [Slides]

Goals: 1) Evaluate "what if” scenarios and mitigation strategies to combat COVID-19 spread within LANL and 2) perform retrospective analysis to assess accuracy of LANL COVID-19 agent-based model (ABM).

60 APPLIED LIFE SCIENCES↗

An Efficient End-To-End Model for the Kepler Photometer

The Kepler Mission is designed to characterize the frequency of Earth-sized planets in the habitable zones of solar-like stars in the solar galactic neighborhood by observing >100,000 main-sequence stars in a >100 square degree field of view (FOV) and seeking evidence of transiting planets. As part of the system engineering effort, we have developed an End-To-End Model (ETEM) of the photometer to better characterize the expected performance of the instrument and to guide us in making design trades. This model incorporates engineering information such as the point spread function, time histories of pointing offsets, operating temperature, quantization noise, the effects of shutterless readout, and read noise. Astrophysical parameters, such as a realistic distribution of stars vs. magnitude for the chosen FOV, zodiacal light, and cosmic ray events are also included. For a given set of design and operating parameters, ETEM generates pixel time series for all pixels of interest for a single CCD channel of the photometer. These time series are then processed to form light curves for the target stars and the impact of various noise sources on the combined differential photometric precision can be determined. This model is of particular value when investigating the effects of noise sources that cannot be easily subjected to direct analysis, such as residual pointing offsets, thermal drift or cosmic ray effects. This version of ETEM features extremely efficient computation times relative to the previous version while maintaining a high degree of fidelity with respect to the realism of the relevant phenomena.

performance modeling↗

Diffusive and Radiative Transport in Fires (DARTFire): Opposed-Flow Flame Spread in Low-Velocity Flows

For flames spreading into a low-velocity flow that can only be obtained in microgravity, we have observed behavior that is different from that which is obtained at higher velocities where radiative effects are unimportant and species transport is relatively fast. Unfortunately, lack of a large body of low-gravity flame spread experimental data inhibits progress in developing an understanding of the physics of low-velocity, opposed-flow flame spread phenomena. Recent DARTFire sounding rocket experimental studies though, coupled with developing theory and modelling, have allowed some strides in understanding to be made, on which we report here. Four launches to date have resulted in six experiments for opposed-flow flame spread over a thick PMMA sample. During the 6 min microgravity period, the PMMA samples were ignited, and steady flame spread was studied under varied flow velocity, oxidizer atmospheric conditions, and, because radiative heat transfer is so important in these slowly spreading flames, external radiant flux. These were the first attempts at such experimental control and measurement in microgravity. A recent reflight of the Solid Surface Combustion Experiment (SSCE) has demonstrated, as modelling predicts, that for the thick fuel of the DARTFire experiment, flame spread in a quiescent environment is a transient process evolving from ignition to extinction on the order of 600 s (Altenkirch et al., 1999). Further study then of the effects of radiation in the very low-velocity opposing flows is of interest in understanding the transition from steady, sustained spread to the unsteady evolution to extinction as the opposing flow is reduced further and eventually removed.

Altenkirch, R. A.↗

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material in the limit of long and moderately long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗

Novel Analysis of the Multiwavelength Structure of the Relativistic Jet in Quasar 3C 273

We present a detailed analysis of the best-quality multiwavelength data gathered for the large-scale jet in the core-dominated quasar 3C 273. We analyze all the archival observations of the target with the Chandra X-ray Observatory, the far-ultraviolet observations with the Hubble Space Telescope, and the 8.4 GHz map obtained with the Very Large Array. In our study, we focus on investigating the morphology of the outflow at different frequencies, and therefore we apply various techniques for the image deconvolution, paying particular attention to a precise modeling of the Chandra and Hubble point-spread functions. We find that the prominent brightness enhancements in the X-ray and far-ultraviolet jet of 3C 273—the “knots”—are not point-like, and can be resolved transversely as extended features with sizes of about ≃0.5 kpc. Also, the radio outflow is wider than the deconvolved X-ray/ultraviolet jet. We have also found circumstantial evidence that the intensity peaks of the X-ray knots are located systematically upstream of the corresponding radio intensity peaks, with the projected spatial offsets along the jet ranging from ≲0.2 kpc up to ≃1 kpc. We discuss our findings in the wider context of multi-component models for the emission and structure of large-scale quasar jets, and speculate on the physical processes enabling an efficient acceleration of the emitting ultrarelativistic electrons along the entire jet length that exceeds 100 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗

Northeast US Ecological Forecasting: Modeling Invasive Plant Habitat Suitability to Support Management Efforts in the American Northeast

Invasive plant species threaten environmental and economic interests when they spread into new areas, outcompete native species, and disrupt ecosystem services. If the spread is not controlled early, species can become well-established and increasingly difficult to manage. The National Park Service (NPS) Invasive Plant Management Teams (IPMTs) strive for an “early detection, rapid response” approach to reducing invasive species spread. Management teams can better prioritize their work with the help of species distribution models (SDMs), which map habitat suitability by combining species occurrences with environmental predictor variables. Scarce invaded range data for newly arrived invasive species presents a particular challenge for producing accurate models. To improve future modeling efforts, this project compared SDM methods using different spatial scales to model two plant species invasive to the Northeast US: the well-established Japanese stiltgrass (Microstegium vimineum) and newer invasive species wavyleaf basketgrass (Oplismenus undulatifolius). The team used NASA Earth observations and climate datasets to model occurrence data and predictor layers at a US-specific extent (90m2 spatial resolution) and global extent (1 km2 spatial resolution). Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), and Landsat 8 Operational Land Imager (OLI) provided data for US Normalized Difference Moisture Indices (NDMI), while global NDMI and topographic predictor layers were derived from Shuttle Radar Topography Mission (SRTM) and Terra Moderate Resolution Imaging Spectroradiometer (MODIS). The resulting models indicated important predictor variables for each species and explored the benefits and tradeoffs of using global data to model habitat suitability for new-arrival invasive species.

Rebecca Ohman↗

A nonlinear model for gas chromatograph systems

Fundamental engineering design techniques and concepts were studied for the optimization of a gas chromatograph-mass spectrometer chemical analysis system suitable for use on an unmanned, Martian roving vehicle. Previously developed mathematical models of the gas chromatograph are found to be inadequate for predicting peak heights and spreading for some experimental conditions and chemical systems. A modification to the existing equilibrium adsorption model is required; the Langmuir isotherm replaces the linear isotherm. The numerical technique of Crank-Nicolson was studied for use with the linear isotherm to determine the utility of the method. Modifications are made to the method eliminate unnecessary calculations which result in an overall reduction of the computation time of about 42 percent. The Langmuir isotherm is considered which takes into account the composition-dependent effects on the thermodynamic parameter, mRo.

Feinberg, M. P.↗