Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 271 records · Page 15

Ultraviolet variability of NGC 5548 - Dynamics of the continuum production region and geometry of the broad-line region

Data from the 1989-1990 IUE monitoring of the Seyfert galaxy NGC 5548 are used here to analyze the continuum variability properties of the galaxy and to derive the structure or its emission-line region. The mean shape of the UV continuum is well fit by an accretion disk model with a given black hole mass and an additional component required to reproduce the observed soft X-ray flux. The continuum fluctuation power spectrum is very steep, with most of the variance coming from about 1 yr time scales. The entire optical/UV continuum rises and falls almost simultaneously, so that the logarithmic slope of the power spectrum is nearly the same for all bands, but the flux at higher photon frequencies varies with larger amplitude. The emission-line material around the nucleus may best be described by a highly ionized inner zone of high and nearly constant pressure that stretches about 4-14 light-days from the center and an outer, more weakly ionized zone of considerably lower ionization at least 20-30 light-days out.

Krolik, J. H.↗

IUE observations of eight dwarf novae - A study of the outburst cycle from 0.12 to 3.5 microns

The dwarf novae are a class of close binary systems with a G-M star transferring mass to a white dwarf. The transferred material accumulates in a disk and accretes onto the white dwarf. Dwarf novae are characterized by outbursts on time scales of weeks to months possibly triggered by instabilities in the disk. To test the steady state models and the behavior of the disk at various stages of the outburst cycle, IUE observations were conducted for eight dwarf novae. Selected UV continuum points were joined to nonsimultaneous optical and IR spectrophotometry and broad-band photography to yield flux distributions over a large wavelength range. These flux distributions are discussed in terms of available accretion disk models.

Szkody, P.↗

A general mechanistic framework for cross-scale understanding of hot spots and hot moments in carbon and water fluxes

Semi-arid ecosystems, like those in the American Southwest, exert a massive impact on the interannual variability of carbon and water cycling. Unfortunately, these carbon and water fluxes are notoriously difficult to predict due to their high spatial and temporal variability, which is poorly captured by the current generation of vegetation models. Indeed, this region is exemplified by the ‘hot spots and hot moments’ concept, which states that small areas in space (‘hot spots’) and transient moments in time (‘hot moments’) exert an outsized influence on biogeochemical cycling. However, the factors that regulate these pulses in biogeochemical activity are unknown, as is their variability across space and time. These uncertainties severely limit efforts to better represent hot spots and hot moments in models. Here, we seek to develop a generalized method for detecting and quantifying the importance of hot spots and hot moments from individual plant to regional scales. Underpinning this method is our recently developed statistical approach for identifying hot spots and hot moments. By applying this method to semi-continuous measurements of plant water status, a depth profile of soil water potential, and ecosystem fluxes via eddy covariance, we will track the fate of water through the soil-plant-atmosphere continuum and identify the mechanistic drivers of these transient pulses in biogeochemical activity. Then, we will expand this approach across a broad network of Ameriflux towers, and apply a machine learning approach that will allow us to upscale measurements of hot spots and hot moments across the American Southwest and quantify their impact on carbon and water cycles. These products will allow us to identify hot spots and hot moments across spatio-temporal scales and will serve as crucial data sources for validating a new generation of models that can better capture highly dynamic carbon and water fluxes. The proposed method will be easily transferable across biomes and will serve as a framework for future research on hot spots and hot moments across the plant ecophysiology, biometeorology, and vegetation modeling communities.

54 ENVIRONMENTAL SCIENCES↗

A radiative transfer model for surface radiation budget studies

A high resolution spectral radiative transfer model of the troposphere for computing downwelling radiance and flux density at the surface in the 5-200 micron region has been developed. The model is fast and accurate and takes into consideration all major and minor constituents active in this region, namely, water vapor, carbon dioxide, ozone, nitrous oxide, and methane. It uses the quasi-random band model of absorption for computing atmospheric transmittances due to various bands. Continuum absorption by water vapor in the 8-14 micron region has also been taken into account. The model has been used to examine the sensitivity of downwelling flux to variations of surface water vapor pressure (or surface relative humidity), water vapor scale-height and concentrations of other gaseous constituents of the atmosphere. In addition, sensitivity to changes in a number of cloud parameters, namely, fractional cloud cover, cloud height, and high-cloud emissivity, has also been examined. Results obtained with the present model compare well with those from the empirical and semi-empirical correlations. The present work demonstrates, for the first time, the feasibility of using a high-resolution narrow-band model for computing this flux for a meteorological application.

Gupta, S. K.↗

Optical polarization of a complete sample of radio sources

The radio and optical properties of a complete 5 GHz sample of strong radio sources are studied. High polarization is strongly correlated with the fraction of the total 5 GHz flux density found in a milliarcsec core. The detection rate of radio galaxies containing optical cores with high polarization is consistent with the hypothesis that all radio sources contain optically polarized cores, with strength proportional to the core radio flux density. High optical polarization, optical power-law fraction, line-to-continuum ratio, emission lines of small equivalent widths, and large amplitude flux variability are all strongly correlated with the fraction of the 5 GHz radio flux density that is unresolved on VLBI scales. The distributions of these properties are well matched by a model where radio and optical radiation have the same beaming geometry and with Doppler boosting of a factor of 10-100 in the optical. A striking alignment of the position angle of the VLBI structure axis and the position angle of optical polarization in highly polarized sources is confirmed.

Impey, C. D.↗

Continuum, cyclotron line, and absorption variability in the high-mass X-ray binary Vela X-1

Because of its complex clumpy wind, prominent cyclotron resonant scattering features, intrinsic variability, and convenient physical parameters (close distance, high inclination, and small orbital separation), which facilitate the observation and analysis of the system, Vela X-1 is one of the key systems for understanding accretion processes in high-mass X-ray binaries on all scales. We revisit Vela X-1 with two new observations taken with NuSTAR at orbital phases ∼0.68–0.78 and ∼0.36–0.52, which show a plethora of variability and allow us to study the accretion geometry and stellar wind properties of the system. We follow the evolution of spectral parameters down to the pulse period timescale using a partially covered power law continuum with a Fermi-Dirac cutoff to model the continuum and local absorption. We are able to confirm anti-correlations between the photon index and the luminosity and, for low fluxes, between the folding energy and the luminosity, implying a change of properties in the Comptonising plasma. We were not able to confirm a previously seen correlation between the cyclotron line energy and the luminosity of the source in the overall observation, but we observed a drop in the cyclotron line energy following a strong flare. We see strong variability in absorption between the two observations and within one observation (for the ∼0.36–0.52 orbital phases) that can be explained by the presence of a large-scale structure, such as accretion and photoionisation wakes in the system, and our variable line of sight through this structure.

C. M. Diez↗

Static adhesion hysteresis in elastic structures

Adhesive interactions between elastic structures such as graphene sheets, carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to irrecoverable energy loss associated with bond breakage, even in static (rate-independent) experiments. To understand this phenomenon, we start with a minimal theory for the peeling of a thin sheet from a substrate, coupling the local event of bond breaking to the nonlocal elastic relaxation of the sheet and show that this can drive static adhesion hysteresis over a bonding/debonding cycle. Using this model we quantify hysteresis in terms of the adhesion and elasticity parameters of the system. We find that this allows us to derive a scaling relation that preserves hysteresis at different levels of granularity while resolving a seeming paradox of lattice trapping in the continuum limit of a discrete fracture process. Finally, to verify our theory, we use new experiments to demonstrate and measure adhesion hysteresis in bundled microtubules.

36 MATERIALS SCIENCE↗

The central engine of quasars and active galactic nuclei Hadronic interactions of shock-accelerated relativistic protons

The detailed dynamics are presented of a model for AGNs and QSOs based on spherically symmetric accretion onto a black hole and dissipation of the accretion energy in a collisionless shock. The basic assumptions and equations relevant for the model are presented, and the relations between the physical parameters are given in terms of the shock parameters. The latter are calculated for varying Mach numbers, and the resulting scaling laws are discussed and compared to the data, in particular the recently discovered correlations between luminosity and mass. The model can naturally produce ultrarelativistic electrons with flat spectra, a feature which, when coupled with e(+)-e(-) pair production, can account for the observed continuum radiation from these objects over essentially the whole electromagnetic spectrum. The model also predicts the emission of high-energy neutrinos with a luminosity roughly equal to that of the photons.

Kazanas, D.↗

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)↗

Intermittent behavior of galactic dynamo activities

Recent observations by Beck and Golla of far-infrared and radio continuum emission from nearby spiral galaxies suggest that the galactic magnetic field strength is connected to the current star formation rate. The role of star formation on the generation of large-scale galactic magnetic field is studied in this paper. Using a simple galactic model, it is shown how the galactic dynamo depends strongly on the turbulent velocity of the interstellar medium. When the star formation efficiency is high, the ISM is churned which in turn amplifies the galactic magnetic field. Between active star formation epochs, the magnetic field is in dormant state and decays at a negligible rate. If density waves trigger star formation, then they also turn on the otherwise dormant dynamo.

Ko, C. M.↗

Extended JT supergravity and random matrix models: The power of the string equation

A number of supersymmetric Jackiw-Teitelboim (JT) gravity theories are known to be described (in the Euclidean path integral formulation) by double-scaled random matrix models. Such matrix models can be characterized using a certain “string equation”. It was shown recently that in extended supergravity, when the number of BPS states scales as e$^{S_0}$, where $S_0$ is the extremal entropy, a special ansatz for the leading order solution of the string equation yields the supergravity spectrum. Somewhat miraculously, the construction showed that the functional form of the non-BPS (continuum) sector predicts the precise form of the BPS sector, showing the robustness of the supergravity/matrix-model correspondence. In this paper, we refine the analysis and show that the string equation, combined with some simple requirements on solutions, are powerful tools for constraining the spectrum of extended JT supergravity theories. We re-explore the cases of $\mathcal{N} = 2$ and (small) $\mathcal{N} = 4$ JT supergravity, and then explore the new cases of spectra from $\mathcal{N} = 3$ and $\mathcal{N} = 4$ large JT supergravity (recently derived by Heydeman, Shi, and Turiaci) showing that our approach also works naturally for (nearly) all the models. Based on this success, we conjecture that these new supergravity models also have matrix model descriptions.

Extended Supersymmetry↗

COWALKER:EFFECTIVE TRANSPORT PROPERTIES OF COMPOSITE MATERIALS

SF-23-026 This software computes effective transport properties of composite materials involving fibers and nanoparticles using a random-walk algorithm that efficiently scales to an arbitrary number of processes and cores. Effective transport properties (thermal, electrical) are key to bridge the microstructure of complex materials with its macroscopic behavior. Traditional approaches either use effective medium approximations (closed mathematical expressions that are approximation for certain conditions) or continuum simulation models such as finite element or finite volume, which require the generation of a mesh for each configuration explored. cowalker leverages the equivalence between laplacian or heat equation-based models and random walks to compute the asymptotic transport properties from an ensemble of first sojourn times of a random walker moving through the composite material. This allows us to directly define a composite material as a collection of particles and use algorithms developed for molecular dynamics to quickly compute the intersection of the walker with the different interfaces in the material. cowalker is developed in C++, and it relies on the GNU Scientific Library for random generation. cowalker is currently delivered as source code, so the GSL library is not included in cowalker's distribution. A more userfriendly version, cowalker.jl is currently in development and will be released as part of cowalker.

YANGUAS-GIL, ANGEL↗

Cloud-particle galactic gas dynamics and star formation

Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.

Roberts, W. W., Jr.↗

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Interpretation of the microlensing event in QSO 2237 + 0305

A model of microlensing for image A of the gravitationally lensed QSO 2237 + 0305 for which Irwin et al. reported in 1989 an increase of the apparent luminosity by about 0.5 mag on a time scale of a few months is presented. The model, with the Salpeter mass function over the mass range of 0.1-1.0 solar mass and the transverse velocity of the lens (or observer) of 600 km/s, can reproduce the reported luminosity variation if the source of the optical continuum has a radius smaller than about 2 x 10 to the 15th cm. This size is compatible with the accretion disk interpretation of the big ultraviolet bump in quasar spectra. The model demonstrates a very large diversity of light curves while the source crosses individual microcaustics or clusters of microcaustics. It will take more than 100 yr before the full variety of light curves will be sampled by the observations.

Wambsganss, J.↗

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗