Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Fast time-resolved electrostatic force microscopy: Achieving sub-cycle time resolution

The ability to measure microsecond- and nanosecond-scale local dynamics below the diffraction limit with widely available atomic force microscopy hardware would enable new scientific studies in fields ranging from biology to semiconductor physics. However, commercially available scanning-probe instruments typically offer the ability to measure dynamics only on time scales of milliseconds to seconds. Here, we describe in detail the implementation of fast time-resolved electrostatic force microscopy using an oscillating cantilever as a means to measure fast local dynamics following a perturbation to a sample. We show how the phase of the oscillating cantilever relative to the perturbation event is critical to achieving reliable sub-cycle time resolution. We explore how noise affects the achievable time resolution and present empirical guidelines for reducing noise and optimizing experimental parameters. Specifically, we show that reducing the noise on the cantilever by using photothermal excitation instead of piezoacoustic excitation further improves time resolution. We demonstrate the discrimination of signal rise times with time constants as fast as 10 ns, and simultaneous data acquisition and analysis for dramatically improved image acquisition times.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

The magnetic recoil spectrometer (MRSt) for time-resolved measurements of the neutron spectrum at the National Ignition Facility (NIF)

The next-generation magnetic recoil spectrometer for time-resolved measurements of the neutron spectrum has been conceptually designed for the National Ignition Facility. This spectrometer, called MRSt, represents a paradigm shift in our thinking about neutron spectrometry for inertial confinement fusion applications, as it will provide simultaneously information about the burn history and time evolution of areal density (ρR), apparent ion temperature (T{sub i}), yield (Y{sub n}), and macroscopic flows during burn. From this type of data, an assessment of the evolution of the fuel assembly, hotspot, and alpha heating can be made. According to simulations, the MRSt will provide accurate data with a time resolution of ∼20 ps and energy resolution of ∼100 keV for total neutron yields above ∼10{sup 16}. At lower yields, the diagnostic will be operated at a higher-efficiency, lower-energy-resolution mode to provide a time resolution of ∼20 ps.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Real time time-dependent density functional theory using higher order finite-element methods

We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real time using higher order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time step. We employ spectral finite elements along with special reduced order quadrature to render the overlap matrix diagonal, thereby simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator, wherein the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher order finite elements over linear finite elements, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, we also demonstrate the competence of higher order finite elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Real-Time High-Accuracy Digital Wireless Time, Frequency, and Phase Calibration for Coherent Distributed Antenna Arrays

his work presents a fully-digital high-accuracy real-time calibration procedure for frequency and time alignment of open-loop wirelessly coordinated coherent distributed antenna array (CDA) modems, enabling radio frequency (RF) phase coherence of spatially separated commercial off-the-shelf (COTS) software-defined radios (SDRs) without cables or external references such as the global navigation satellite system (GNSS). Building on previous work using high-accuracy spectrally-sparse time of arrival (ToA) waveforms and a multistep ToA refinement process, a high-accuracy two-way time transfer (TWTT)-based time–frequency coordination approach is demonstrated. Due to the two-way nature of the high-accuracy TWTT approach, the time and frequency estimates are Doppler and multipath tolerant, so long as the channel is reciprocal over the synchronization epoch. This technique is experimentally verified using COTS SDRs in a lab environment in static and dynamic scenarios and with significant multipath scatterers. Time, frequency, and phase stability were evaluated by beamforming over coaxial cables to an oscilloscope which achieved time and phase precisions of ~60– 70 ps , with median coherent gains above 99% using optimized coordination parameters, and a beamforming frequency root-mean-square error (RMSE) of 3.73 ppb in a dynamic scenario. Finally, experiments were conducted to compare the performance of this technique with previous works using an analog continuous-wave two-tone (CWTT) frequency reference technique in both static and dynamic settings.

Clock synchronization↗
DOE OSTI2023

Assessment and Commissioning of Electrical Substation Grid Testbed with a Real-Time Simulator and Protective Relays/Power Meters in the Loop

Electrical utility substations are wired with intelligent electronic devices (IEDs), such as protective relays, power meters, and communication switches. Substation engineers commission these IEDs to assess the appropriate measurements for monitoring, control, power system protection, and communication applications. Like real electrical utility substations, complex electrical substation grid testbeds (ESGTs) need to be assessed for measuring current and voltage signals in monitoring, power system protection, control (synchro check), and communication applications that are limited by small measurement percentage errors. In the process of setting an ESGT with real-time simulators and IEDs in the loop, protective relays, power meters, and communication devices must be commissioned before running experiments. In this study, an ESGT with IEDs and distributed ledger technology was developed. The ESGT with a real-time simulator and IEDs in the loop was satisfactorily assessed and commissioned. The commissioning and problem-solving tasks of the testbed are described to define a method with flowcharts to assess possible trouble-shooting in ESGTs. This method was based on comparing the simulations versus IED measurements for the phase current and voltage magnitudes, three-phase phasor diagrams, breaker states, protective relay times with selectivity coordination at electrical faults, communication data points, and time-stamp sources.

Piesciorovsky, Emilio↗
DOE OSTI2024

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI↗
DOE OSTI2019

Utility of Atmospheric Transport Runs Done Backwards in Time for Source Term Estimation

In a general sense, the equations describing transport of an atmospheric contaminant are reversible in time. However, the time reversibility is not perfect. For a source term estimation problem, one atmospheric transport model (ATM) run done backwards in time to represent the influence of a single sample value can yield as much information as hundreds or thousands of forwards ATM runs. Thus, backwards ATM are often used to limit the computational burden of an analysis. A large number of paired ATM runs at 14 release locations were used to explore the agreement between dilution factors for the forward and backwards time directions. Four of the release locations have correlation coefficients between forward and backwards dilution factors greater than 0.9. None of the four have large elevation changes within 100 km of the release location. The two release locations with the lowest correlation coefficients are in mountainous terrain with large elevation differences. Elevation changes have a greater negative influence on the correlations than distance to the sampler or plume transit times. The locations of the releases were estimated using dilution factors from backwards ATM runs. Many of the location estimates are close to the release point, vindicating the general approach of using ATM runs backwards in time. Aggregated over the 14 release points, 26% of the location estimates were within 10 km of the release point, 61% were within 25 km, 80% were within 50 km, and 87% were within 75 km. Most of the remaining 13% (poor matches) occur for two release locations in mountainous terrain. Good time-reversibility cannot be guaranteed for a new analysis, so we recommend that any source term solution using backwards ATM runs should include comparisons based on forward ATM runs. As an example, calculations show good agreement between dilutions in different time-directions for samples taken after the announced nuclear test by the Democratic People’s Republic of Korea in 2013.

Eslinger, Paul W.↗
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2018

Fast Real-Time Time-Dependent Density Functional Theory Calculations with the Parallel Transport Gauge

Real-time time-dependent density functional theory (RT-TDDFT) is known to be hindered by the very small time step (attosecond or smaller) needed in the numerical simulation, because of the fast oscillation of electron wave functions, which significantly limits its range of applicability for the study of ultrafast dynamics. In this paper, we demonstrate that such oscillation can be considerably reduced by optimizing the gauge choice using the parallel transport formalism. RT-TDDFT calculations can thus be significantly accelerated using a combination of the parallel transport gauge and implicit integrators, and the resulting scheme can be used to accelerate any electronic structure software that uses a Schrödinger representation. Here, using absorption spectrum, ultrashort laser pulse, and Ehrenfest dynamics calculations for example, we show that the new method can utilize a time step that is on the order of 10–100 attoseconds using a planewave basis set. Thanks to the significant increase of the size of the time step, we also demonstrate that the new method is more than 10 times faster, in terms of the wall clock time, when compared to the standard explicit fourth-order Runge–Kutta time integrator for silicon systems ranging from 32 to 1024 atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

DETECTION OF BROAD Hα EMISSION LINES IN THE LATE-TIME SPECTRA OF A HYDROGEN-POOR SUPERLUMINOUS SUPERNOVA

iPTF13ehe is a hydrogen-poor superluminous supernova (SLSN) at z = 0.3434, with a slow-evolving light curve and spectral features similar to SN2007bi. It rises in 83–148 days to reach a peak bolometric luminosity of ∼1.3 × 10{sup 44} erg s{sup −1}, then decays slowly at 0.015 mag day{sup −1}. The measured ejecta velocity is ∼ 13,000 km s{sup −1}. The inferred explosion characteristics, such as the ejecta mass (70–220 M{sub ⊙}), and the total radiative and kinetic energy (E{sub rad} ∼ 10{sup 51} erg, E{sub kin} ∼ 2 × 10{sup 53} erg), are typical of slow-evolving H-poor SLSN events. However, the late-time spectrum taken at +251 days (rest, post-peak) reveals a Balmer Hα emission feature with broad and narrow components, which has never been detected before among other H-poor SLSNe. The broad component has a velocity width of ∼4500 km s{sup −1} and a ∼300 km s{sup −1} blueward shift relative to the narrow component. We interpret this broad Hα emission with a luminosity of ∼2 × 10{sup 41} erg s{sup −1} as resulting from the interaction between the supernova ejecta and a discrete H-rich shell, located at a distance of ∼4 × 10{sup 16} cm from the explosion site. This interaction causes the rest-frame r-band LC to brighten at late times. The fact that the late-time spectra are not completely absorbed by the shock-ionized H-shell implies that its Thomson scattering optical depth is likely ≤1, thus setting upper limits on the shell mass ≤30 M{sub ⊙}. Of the existing models, a Pulsational Pair Instability supernova model can naturally explain the observed 30 M{sub ⊙} H-shell, ejected from a progenitor star with an initial mass of (95–150) M{sub ⊙} about 40 years ago. We estimate that at least ∼15% of all SLSNe-I may have late-time Balmer emission lines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Time-Local Equation for the Exact Optimized Effective Potential in Time-Dependent Density Functional Theory

A long-standing challenge in the time-dependent density functional theory is to efficiently solve the exact time-dependent optimized effective potential (TDOEP) integral equation derived from orbital-dependent functionals, especially for the study of nonadiabatic dynamics in time-dependent external fields. In this Letter, we formulate a completely equivalent time-local TDOEP equation that admits a unique real-time solution in terms of time-dependent Kohn-Sham and effective memory orbitals. The time-local formulation is numerically implemented, with the incorporation of exponential memory loss to address the unaccounted for correlation component in the exact-exchange-only functional, to enable the study of the many-electron dynamics of a one-dimensional hydrogen chain. It is shown that the long time behavior of the electric dipole converges correctly and the zero-force theorem is fulfilled in the current implementation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Temperature sensitivity of soil organic carbon decomposition increased with mean carbon residence time: Field incubation and data assimilation

Abstract Temperature sensitivity of soil organic carbon ( SOC ) decomposition is one of the major uncertainties in predicting climate‐carbon (C) cycle feedback. Results from previous studies are highly contradictory with old soil C decomposition being more, similarly, or less sensitive to temperature than decomposition of young fractions. The contradictory results are partly from difficulties in distinguishing old from young SOC and their changes over time in the experiments with or without isotopic techniques. In this study, we have conducted a long‐term field incubation experiment with deep soil collars (0–70 cm in depth, 10 cm in diameter of PVC tubes) for excluding root C input to examine apparent temperature sensitivity of SOC decomposition under ambient and warming treatments from 2002 to 2008. The data from the experiment were infused into a multi‐pool soil C model to estimate intrinsic temperature sensitivity of SOC decomposition and C residence times of three SOC fractions (i.e., active, slow, and passive) using a data assimilation ( DA ) technique. As active SOC with the short C residence time was progressively depleted in the deep soil collars under both ambient and warming treatments, the residences times of the whole SOC became longer over time. Concomitantly, the estimated apparent and intrinsic temperature sensitivity of SOC decomposition also became gradually higher over time as more than 50% of active SOC was depleted. Thus, the temperature sensitivity of soil C decomposition in deep soil collars was positively correlated with the mean C residence times. However, the regression slope of the temperature sensitivity against the residence time was lower under the warming treatment than under ambient temperature, indicating that other processes also regulated temperature sensitivity of SOC decomposition. These results indicate that old SOC decomposition is more sensitive to temperature than young components, making the old C more vulnerable to future warmer climate.

Zhou, Xuhui↗
DOE OSTI2014

Cardiac C-arm computed tomography using a 3D + time ROI reconstruction method with spatial and temporal regularization

Purpose: Reconstruction of the beating heart in 3D + time in the catheter laboratory using only the available C-arm system would improve diagnosis, guidance, device sizing, and outcome control for intracardiac interventions, e.g., electrophysiology, valvular disease treatment, structural or congenital heart disease. To obtain such a reconstruction, the patient's electrocardiogram (ECG) must be recorded during the acquisition and used in the reconstruction. In this paper, the authors present a 4D reconstruction method aiming to reconstruct the heart from a single sweep 10 s acquisition. Methods: The authors introduce the 4D RecOnstructiOn using Spatial and TEmporal Regularization (short 4D ROOSTER) method, which reconstructs all cardiac phases at once, as a 3D + time volume. The algorithm alternates between a reconstruction step based on conjugate gradient and four regularization steps: enforcing positivity, averaging along time outside a motion mask that contains the heart and vessels, 3D spatial total variation minimization, and 1D temporal total variation minimization. Results: 4D ROOSTER recovers the different temporal representations of a moving Shepp and Logan phantom, and outperforms both ECG-gated simultaneous algebraic reconstruction technique and prior image constrained compressed sensing on a clinical case. It generates 3D + time reconstructions with sharp edges which can be used, for example, to estimate the patient's left ventricular ejection fraction. Conclusions: 4D ROOSTER can be applied for human cardiac C-arm CT, and potentially in other dynamic tomography areas. It can easily be adapted to other problems as regularization is decoupled from projection and back projection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Stabilization of Switched Systems on Non Uniform Time Domain with Dwell Time

In this paper, we will present new stability conditions for a special class of linear switched systems, that evolves on non-uniform time domain. The considered systems switches between continuous-time subsystems on intervals with variable lengths and discrete-time subsystems with variable step sizes. Time scale theory is introduced to derive conditions for exponential stability of this special class of switched systems by using the dwell time approach. The conditions are based on the existence of a Lyapunov function which is non-increasing at the switching instants. This shown that this class of switched systems can be stabilized if the dwell time of each continuous-time subsystem is greater than some bound, and if the gap of the discrete-time subsystem is bounded by some specific values. Numerical examples are presented to show the effectiveness of the proposed scheme.

Taousser, Fatima Zohra↗
DOE OSTI2017

Late-time Cooling of Neutron Star Transients and the Physics of the Inner Crust

An accretion outburst onto a neutron star transient heats the neutron star’s crust out of thermal equilibrium with the core. After the outburst, the crust thermally relaxes toward equilibrium with the neutron star core, and the surface thermal emission powers the quiescent X-ray light curve. Crust cooling models predict that thermal equilibrium of the crust will be established ≈1000 days into quiescence. Recent observations of the cooling neutron star transient MXB 1659-29, however, suggest that the crust did not reach thermal equilibrium with the core on the predicted timescale and continued to cool after ≈2500 days into quiescence. Because the quiescent light curve reveals successively deeper layers of the crust, the observed late-time cooling of MXB 1659-29 depends on the thermal transport in the inner crust. In particular, the observed late-time cooling is consistent with a low thermal conductivity layer near the depth predicted for nuclear pasta that maintains a temperature gradient between the neutron star’s inner crust and core for thousands of days into quiescence. As a result, the temperature near the crust–core boundary remains above the critical temperature for neutron superfluidity, and a layer of normal neutrons forms in the inner crust. We find that the late-time cooling of MXB 1659-29 is consistent with heat release from a normal neutron layer near the crust–core boundary with a long thermal time. We also investigate the effect of inner crust physics on the predicted cooling curves of the accreting transient KS 1731-260 and the magnetar SGR 1627-41.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Physical Origin of Long Gas Depletion Times in Galaxies

We present a model that explains why galaxies form stars on a timescale significantly longer than the timescales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short timescale, and thus the rate of this evolution does not limit the star formation rate (SFR). Instead, the SFR is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of a star-forming state multiple times, which results in a long timescale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the SFR in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated L {sub *}-sized galaxy simulation that reproduces the observed Kennicutt–Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density.

79 ASTRONOMY AND ASTROPHYSICS↗