SEARCH · Engineering Papers
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.
Low dispersion IUE spectra of 3 cataclysmic variables (V442 Oph, V794 Aq1 and H2215-086) are discussed in terms of current disk models. The range of continuum fluxes, line emission and disk parameters of these three novalike systems are compared with past observations of dwarf novae at outburst and quiescence. Evidence of variability on orbital time scales is presented for V442 Oph and H2215-086.
Observations were obtained at wavelengths from 1.5 to 7.5 microns with beams varying in diameter from 4 to 28 arcsec of infrared hydrogen recombination lines toward the Mon R2 IRS1 H II region. It is found that the data cannot be fitted with the extinction law which characterizes the interstellar medium unless the obscuring matter is clumped on a small scale of not greater than 0.3 arcsec; in which case considerable fluctuations in the amount of extinction on scales smaller than 1 arcsec are expected. The data of Simon et al. (1983) suggest a dip in the extinction about 5 arcsec from the 2-micron and radio continuum peak, and rule out models with uniform dust and clump distributions.
Analytic solutions are derived for simple models of equatorial wave-CISK with conditional (positive-only) heating. Conditional heating favors an unstable mode having a single wet region of infinitesimal width. This result helps explain model simulations of tropical superclusters and intraseasonal oscillations. Mechanisms of scale selection and propagation are discussed, and the model is generalized to include evaporation-wind feedback, multiple vertical-mode CISK, and scale-dependent damping. The first two mechanisms cause zonal propagation, without scale selection. Scale-dependent damping breaks the inviscid continuum (wherein modes exist as a continuous function of wet region width) into discrete modes with unique horizontal structure and phase speed. Results are contrasted with the hierarchy of scales in observed intraseasonal oscillations.
We present a detailed spectral and temporal analysis of a long ROSAT PSPC pointing centered on the bright Seyfert 1 galaxy NGC 5548. The spectrum is particularly complex, with a two-component model favored, confirming the existence of the soft excess in this source. In addition, an absorption feature arising from highly ionized oxygen in the line of sight has been detected, consistent with an origin in the 'warm absorber' which has been found to be a common feature in Seyfert galaxies observed by Ginga. The energy and depth of the feature should allow us to constrain the density and geometry of the absorbing gas, which may lie close to the central continuum source. NGC 5548 was in a low X-ray state during the observation, brightening by a factor of about 2 on a time-scale of about 2 d. Time-resolved spectral analysis shows a change in the ratio of the hard and soft fluxes, consistent with reprocessing models.
Recent observations of the magnetotail plasma sheet have shown it to be a dynamic and turbulent region. Research has found strong turbulence in the plasma sheet at approximately 20 Earth's Radius tailward of Earth; the turbulence is observed at all activity levels. The existence of strong turbulence in the plasma sheet in the region associated with substorm onset might be thought difficult to reconcile with the coherence and repeatability of the substorm cycle. We review a variety of evidence that strongly suggests the magnetotail is driven, through magnetic flux transfer, into a state of "self-organized criticality" (SOC). It is an important property of physical systems that evolve into SOC that they self-organize into a unique, global dynamic state. This global state is inevitable, and repeatable. In this state, however, small-spatiotemporal-scale system fluctuations are unpredictable and can be only described statistically. This is the basis, we think, for the global coherence and repeatability of the substorm phenomenon in the turbulent plasma sheet. At, or near, substorm onset the plasma sheet can be described by a global SOC state containing significant small scale turbulence. In several recent studies, "sandpile" models were driven into SOC and then shown to reproduce various measures of substorm activity. We discuss the plasma physical foundation of these sandpile models. The evolution of simple continuum plasma sheet models into SOC-like states of many small reconnection events in the turbulent plasma sheet under the will be demonstrated. We view the substorm phenomenon as an avalanche assumption that the plasma sheet is in a SOC state.
Validity of the assumptions relating the applicability of continuum shell theories to the global mechanical behavior of carbon nanotubes is examined. The present study focuses on providing a basis that can be used to qualitatively assess the appropriateness of continuum-shell models for nanotubes. To address the effect of nanotube structure on their deformation, all nanotube geometries are divided into four major classes that require distinct models. Criteria for the applicability of continuum models are presented. The key parameters that control the buckling strains and deformation modes of these classes of nanotubes are determined. In an analogy with continuum mechanics, mechanical laws of geometric similitude are presented. A parametric map is constructed for a variety of nanotube geometries as a guide for the applicability of different models. The continuum assumptions made in representing a nanotube as a homogeneous thin shell are analyzed to identify possible limitations of applying shell theories and using their bifurcation-buckling equations at the nano-scale.
Researchers have shown that photoionization by the continuum of the hidden Seyfert I nucleus in NGC 1068 can have a significant effect on the ionization state and energetics of this disk's Interstellar Medium (ISM). Photoionization models with appropriate power law spectra can produce (NII) lambda lambda 6538, 6584/H alpha line ratios of 1.25 for ionization parameters Q approx. 10 (exp -12). However the data indicate large regions where the (NII)/H alpha ratio is 1 to 3. Since the abundances are known to be solar, there must be additional heating sources. Hardening of the incident radiation field by intervening absorption should be able to raise T sub e, thereby raising the (NII)/H alpha ratio. Heating with moderate efficiency by the intense starburst ring should also be a significant factor in raising the temperature of the ISM. The photoionization models with additional heating predict enhanced emission from other forbidden lines including (OII) lambda 3727 and (SII) lambda 6731.
In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.
A methodology for implementing parallel DSMC codes in a heterogeneous computing environment is described. The methodology involves the use of a common message-passing software library together with recently developed software that handles the actual interprocessor communications in a standard manner across a variety of computing platforms. Benchmark tests using a simple DSMC model problem were performed on an Intel iPSC/860, a Cray-YMP and a group of Sun workstations. The approach was found to give speedups that scaled linearly with problem size on all the computing platforms tested. This methodology was then incorporated into a production-type DSMC code to allow the simulation of problems that would not otherwise have been practical. The application of this production code to simulations of hypersonic shear flows and shock-lip interactions under near-continuum conditions is described. Synchronous and asynchronous models for implementing parallelism into DSMC simulations are also described and both models are shown to produce the same steady-state result.
Numerical simulation of the early stages of planet growth show that a few bodies nearly 1000 km in diameter may have formed within approx. 100,000 yr after solid material grew into km scale planetesimals by gravitational instability. Even after such large bodies formed, the bulk of the mass of the future terrestrial planet zone resided in small bodies. Subsequent evolution is difficult to model because it requires simultaneous consideration of continuum (multitudinous small bodies) and discrete (a few large bodies) evolution. Some relevant issues include definition of accretional feeding zones, evaluation of the range of gravitational influence, viscous transport and diffusion, orbital commensurabilities, role of gas, etc. The first large bodies may have been (1) the embryos of the final planets, which grew by accreting tiny planetesimals, or (2) merely the first of many 1000+ km bodies, which grew independently and later collided to form the planets. Models of late stage accretion that assume all bodies to be initially nearly Moon sized provide insight into relevant collisional and dynamical processes. The chief point in this research is that the correct size distribution during the later stages of planet growth remains unknown.
Two scaled models of the Aeroassist Flight Experiment (AFE) vehicle were tested in two air wind tunnels and one CF4 tunnel. The tests were to determine the static longitudinal aerodynamic characteristics, and shock shapes for the configuration in hypersonic continuum flow. The tests were conducted with a range of angle of attack to evaluate the effects of Mach number, Reynolds numbers, and normal shock density ratio.
Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.
Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.
High-speed U-band photometry of the peculiar emission object MWC 560 obtained with the ground-based instrumentation, and V-band photometry obtained with the International Ultraviolet Explorer-Fine Error Sensor indicates irregular brightness variations are quasi-periodic. Multiple peaks of relative brightness power indicate statistically significant quasi periods existing in a range of 3-35 minutes, that are superposed on slower hourly varying components. We present a preliminary model that explains the minute and hourly time-scale variations in MWC 560 in terms of a velocity-shear instability that arises because a white dwarf magnetosphere impinges on an accretion disk. We also find evidence for Fe II multiplet pseudocontinuum absorption opacity in far-UV spectra of CH Cygni which is also present in MWC 560. Both CH Cyg and MWC 560 may be in an evolutionary stage that is characterized by strong UV continuum opacity which changes significantly during outburst, occurring before they permanently enter the symbiotic nebular emission phase.
The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.
A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.
Twisted bilayer MoTe 2 (tMoTe 2 ) is an emergent platform for exploring exotic quantum phases driven by the interplay between nontrivial band topology and strong electron correlations. Direct experimental access to its momentum-resolved electronic structure is essential for uncovering the microscopic origins of the correlated topological phases therein. Here, we report angle-resolved photoemission spectroscopy measurements of tMoTe 2 , revealing pronounced twist-angle-dependent band reconstruction shaped by orbital character, interlayer coupling, and moiré potential modulation. Density functional theory captures the qualitative evolution, yet underestimates key energy scales across twist angles, highlighting the importance of electronic correlations. Notably, the hole effective mass at the 𝐾 point exhibits a nonmonotonic dependence on twist angle, peaking near 2°, consistent with band flattening at the magic angle predicted by continuum models. Via electrostatic gating and surface dosing, we further visualize the evolution of electronic structure versus doping, enabling direct observation of the conduction band minimum and confirm tMoTe 2 as a direct band gap semiconductor. These results establish a spectroscopic foundation for modeling and engineering emergent quantum phases in this moiré platform.
Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.