Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singularities in general relativity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calculation of the Heat Transfer Coefficient in the Outer Body for a Rotational Detonation

Unsteady heat transfer characterization on the combustion surfaces of Rotational Detonation Engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion anulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. The reacting flow fields that are characteristic of RDEs are inherently compressible, three dimensional, unsteady, and turbulent, having properties that change by orders of magnitude throughout the flow field. They will therefore contain multiple length scales and time scales operating everywhere in the flow during all times. Because of this it is not likely that the RDE flow fields will lend themselves to explanation using simple dimensionless parameters. The dimensionless groups have meaning only in situations where length scales and time scales are singular and well defined. In spite of this it may be possible to get a relatively good idea of what the convection heat transfer coefficient is. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the start-up flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE was calculated and compared with experimental data. The closet match can then be used to directly calculate the convection heat transfer coefficient on the inside surface of the RDE.

VanOsdol, John↗
DOE OSTI2023

Constraints on sequential discontinuities from the geometry of on-shell spaces

We present several classes of constraints on the discontinuities of Feynman integrals that go beyond the Steinmann relations. These constraints follow from a geometric formulation of the Landau equations that was advocated by Pham, in which the singularities of Feynman integrals correspond to critical points of maps between on-shell spaces. To establish our results, we review elements of Picard-Lefschetz theory, which connect the homotopy properties of the space of complexified external momenta to the homology of the combined space of on-shell internal and external momenta. An important concept that emerges from this analysis is the question of whether or not a pair of Landau singularities is compatible — namely, whether or not the Landau equations for the two singularities can be satisfied simultaneously. Under conditions we describe, sequential discontinuities with respect to non-compatible Landau singularities must vanish. Although we only rigorously prove results for Feynman integrals with generic masses in this paper, we expect the geometric and algebraic insights that we gain will also assist in the analysis of more general Feynman integrals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Open string amplitudes: singularities, asymptotics and new representations

Open string amplitudes at tree level have been studied for over fifty years. However, there is no known analytic form for general n-point amplitudes, and their conventional representation in terms of worldsheet integrals does not make many of their most basic physical properties manifest. Recently, a formulation of these amplitudes exposing the underlying “binary geometry” via the use of “u” variables has given us many insights into their basic features. In this paper, we initiate a systematic exploration of fundamental aspects of open string amplitudes from this new point of view. We begin by finding explicit expressions for the factorization of amplitudes at general massive levels, which are seen to be determined by products of lower-point massless amplitudes with shifted kinematics. We then study the asymptotic behavior when subsets of kinematic variables become large, delineating regimes with exponential (generalized hard scattering) and power-law (generalized Regge) behavior. We also give precise expressions for the asymptotics, which reveal another example of the recently observed property of factorization away from poles. We derive new recursion relations for the amplitude, which when repeatedly applied reduce to infinite series representations with a wider domain of convergence than the usual integral representations. For the five-point case, we present a new closed-form expression for the amplitude that for the first time gives its analytic continuation to all of kinematic space. We also discuss novel relations between amplitudes at different kinematic points following from the recently observed “split” factorizations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Community Resilience Indicator Analysis: Commonly Used Indicators from Peer-Reviewed Research (Updated for Research Published 2003-2021)

In 2017, FEMA’s National Integration Center (NIC) Technical Assistance (TA) Branch identified a need to establish a data-driven basis for prioritizing locations for TA investment and guiding local emergency management planning. To achieve this goal, FEMA tasked Argonne National Laboratory (Argonne) with identifying commonly used indicators of community resilience across the landscape of published peer-reviewed research. FEMA and Argonne completed the first Community Resilience Indicator Analysis (CRIA) in 2018 and repeated the process in 2022. The CRIA process begins with a literature review and cataloguing of published peer-reviewed assessment methodologies on social vulnerability and community resilience. The literature review findings are then filtered by inclusion criteria established by the CRIA research team to ensure the methodologies are: (1) Quantitative, (2) Data and methodology are publicly available, (3) Calculated at the county level or lower, (4) Examine generalized hazard risk (rather than a singular hazard), and (5) Focused on pre-disaster community conditions. After this, the research team identifies the commonly used indicators across these methodologies and selects the best data source for each indicator. Finally, the research team bins the data for visual display, conducts a correlation analysis and creates a composite index, the FEMA Community Resilience Index (FEMA CRI). In 2018, the CRIA identified eight resilience and vulnerability assessment methodologies and 20 commonly used indicators (indicators used in three or more of the eight methodologies). The FEMA CRI in 2018 was created from these 20 indicators and was produced for at the county level. The 2022 CRIA updated the literature review to expand the list of methodologies examined and followed the same process, resulting in an analysis of 14 methodologies published between 2003 and 2021 and 22 indicators identified as commonly used (indicators used in five or more of the 14 methodologies). In 2022, the research team produced the FEMA CRI at the county and the census tract levels. To make the CRIA data more accessible and more actionable, each individual indicator and the FEMA CRI is binned and included in FEMA’s Resilience Analysis and Planning Tool (RAPT). RAPT enables emergency managers and community partners to quickly visualize relative differences in potential resilience by county, tribe and census tract. By reviewing the data for each of these 22 indicators individually, emergency managers can gain insights for targeted outreach strategies, planning, mitigation investments and response and recovery operations. Communities, regional governments and others can use this data to better understand potential challenges to resilience. As the social science field of examining and validating indicators of resilience evolves, FEMA will update RAPT to provide emergency managers and community partners with additional data and tools to inform planning, mitigation, response and recovery. It is important to understand that the role of the emergency manager is not to change or to “improve” the data, but to plan appropriately for the community characteristics reflected in the data. These datasets are community characteristics that researchers have identified as important considerations for resilience. For example, people with disabilities may have greater challenges to be resilient to disasters. If a community has a high population of people with disabilities, the emergency manager(s) may need to create tailored preparedness outreach programs and strategies to ensure those residents have support if evacuation is necessary. Rather than label these indicators as an absolute measure of resilience, FEMA considers “potential challenges to resilience” a better frame to understand these indicators. Everyone is vulnerable to disasters. While scholars theorize that certain characteristics may make an individual or a household more socially vulnerable, the data does not reflect measures that individuals and/or communities have taken to address potential challenges, such as emergency management planning and outreach or household preparedness measures. To aid emergency managers in understanding how to use these indicators, calling them potential challenges to resilience supports a more positive and strategic application of the data in all phases of emergency management.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Closing the Loops on Solar Photovoltaics Modules: An Agent-Based Modeling Approach for the Study of Circular Economy Strategies

Solar photovoltaics (PV) installed capacity have grown exponentially since the early 2000s (average annual growth rate of 50%). With 4,700 GW projected installed capacity by 2050, the volume of PV panels waste is also expected to become substantial. Though renewables are a sine qua non to the establishment of a truly circular economy (CE), the issue arising from their end-of-life management needs to be resolved. The management of PV end-of-life also represent a singular opportunity to create value, from recovered valuable, rare or critical materials (e.g. silver, tellurium, indium). Moreover, the PV case illustrates some of the current barriers to CE; for instance, the need for a common definition of waste (PV are defined as e-waste in the European Union but as general waste in the United States) and potential loss of innovations' advantages (PV average efficiency has continuously grown the past 10 years). In this study, an agent-based modeling (ABM) approach is proposed to simulate PV end-of-life management in the United States. The model explores how the decisions of the various actors involved in handling PV waste affect the quantities of PV that are reused, recycled, or land-filled. In the model, manufacturers, residential, and nonresidential PV owners, installers, and recyclers are represented by agents while governmental policies and regulations are treated as exogenous variables. PV-market data are used to define agents' characteristics (e.g., installed PV capacities or producing costs). Agents' decision rules draw on the literature related to industrial symbiosis (IS) and peoples' waste behaviors as they appear as the most widely used model to implement CE principles at a meso level. Specifically, the concepts of mutual trust between IS actors and knowledge about the IS philosophy are yielded to define the agents' decision process. The primary outputs of the ABM are the costs and volume of waste associated with each end-of-life pathway. Preliminaries results indicate recycling of residential PV is affected by the cost, perceived difficulty of recycling behaviors, and social norms. Moreover, the type of network connecting agents and the initial recycling rate strongly influence results. The model also highlights the crucial role of recycling behavior adoption: in case of early failure of PV, recycled volumes increase by about 40%. This represents an additional 5.7 billion USD from potential recovered silver. Finally, although results for the PV case are presented here, the developed ABM aims at being a general tool for the study of CE strategies and the CE transition. Furthermore, steps of this research include extending the model to encompass circular economy strategies (e.g., design for recycling or lifetime extension) and validating its general architecture from various case studies.

14 SOLAR ENERGY↗
DOE OSTI2019

Kinematic focus point method for particle mass measurements in missing energy events

We investigate the solvability of the event kinematics in missing energy events at hadron colliders, as a function of the particle mass ansatz. To be specific, we reconstruct the neutrino momenta in dilepton tt¯-like events, without assuming any prior knowledge of the mass spectrum. We identify a class of events, which we call extreme events, with the property that the kinematic boundary of their allowed region in mass parameter space passes through the true mass point. We develop techniques for recognizing extreme events in the data and demonstrate that they are abundant in a realistic data sample, due to expected singularities in phase space. We propose a new method for mass measurement whereby we obtain the true values of the mass parameters as the focus point of the kinematic boundaries for all events in the data sample. Since the masses are determined from a relatively sharp peak structure (the density of kinematic boundary curves), the method avoids some of the systematic errors associated with other techniques. We show that this new approach is complementary to previously considered methods in the literature where one studies the solvability of the kinematic constraints throughout the mass parameter space. In particular, we identify a problematic direction in mass space of nearly 100% solvability, and then show that the focus point method is effective in lifting the degeneracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Regularization prescription for the mixing between nonlocal gluon and quark operators

It is well known that in the study of mixing between nonlocal gluon and quark bilinear operators there exists an ambiguity when relating coordinate space and momentum space results. In this work, we show that this ambiguity is due to the lack of a proper regularization prescription of the singularity that arises when the separation between the gluon/quark fields approaches zero. We then demonstrate that dimensional regularization resolves this issue and yields consistent results in both coordinate and momentum space. This prescription is also compatible with lattice extractions of parton distributions from nonlocal operators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Celestial operator product expansions and w 1+∞ symmetry for all spins

The operator product expansion of massless celestial primary operators of arbitrary spin is investigated. Poincaré symmetry is found to imply a set of recursion relations on the operator product expansion coefficients of the leading singular terms at tree-level in a holomorphic limit. The symmetry constraints are solved by an Euler beta function with arguments that depend simply on the right-moving conformal weights of the operators in the product. These symmetry-derived coefficients are shown not only to match precisely those arising from momentum-space tree-level collinear limits, but also to obey an infinite number of additional symmetry transformations that respect the algebra of w 1+∞ . In tree-level minimally-coupled gravitational theories, celestial currents are constructed from light transforms of conformally soft gravitons and found to generate the action of w 1+∞ on arbitrary massless celestial primaries. Results include operator product expansion coefficients for fermions as well as those arising from higher-derivative non-minimal couplings of gluons and gravitons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Cosmological Bootstrap: Spinning correlators from symmetries and factorization

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ~1.3×10 9 parameters to obtain accurate dynamics.

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

Unidirectional amplification in the frozen mode regime enabled by a nonlinear defect

A stationary inflection point (SIP) is a spectral singularity of the Bloch dispersion relation ω(k) of a periodic structure where the first and the second derivatives of ω with respect to k vanish. An SIP is associated with a third-order exceptional point degeneracy in the spectrum of the unit-cell transfer matrix, where there is a collapse of one propagating and two evanescent Bloch modes. At the SIP frequency, the incident wave can be efficiently converted into the frozen mode with greatly enhanced amplitude and vanishing group velocity. This can be very attractive for applications, including light amplification. Due to its non-resonant nature, the frozen mode regime (FMR) has fundamental advantages over common cavity resonances. Furthermore, we propose, a novel, to the best of our knowledge, scheme for FMR-based unidirectional amplifiers by leveraging a tailored amplification/attenuation mechanism and a single nonlinear defect. The defect breaks the directional symmetry of the periodic structure and enables nonlinearity-related unidirectional amplification/attenuation in the vicinity of the SIP frequency. We demonstrate the robustness of the amplification mechanism to local impurities and parasitic nonlinearity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

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

Past-completeness of inflationary spacetimes

Here, we discuss the question of whether or not inflationary spacetimes can be geodesically complete in the infinite past. Geodesic completeness is a necessary condition for averting an initial singularity during eternal inflation. It is frequently argued that cosmological models which are expanding sufficiently fast (having average Hubble expansion rate H avg > 0) must be incomplete in null and timelike past directions. This well-known conjecture relies on specific bounds on the integral of the Hubble parameter over a past- directed timelike or null geodesic. As stated, we show this claim is an open issue. We show that the calculation of H avg yields a continuum of results for a given spacetime predicated upon the underlying topological assumptions. We present an improved definition for H avg and introduce an uncountably infinite cohort of cosmological solutions which are geodesically complete despite having H avg > 0. We discuss a standardized definition for inflationary spacetimes as well as quantum (semiclassical) cosmological concerns over physically reasonable scale factors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Elastocaloric signatures of symmetric and antisymmetric strain-tuning of quadrupolar and magnetic phases in DyB 2 C 2

The adiabatic elastocaloric effect measures the temperature change of a given system with strain and provides a thermodynamic probe of the entropic landscape in the temperature-strain space. Here, we demonstrate that the DC bias strain-dependence of AC elastocaloric effect allows decomposition of the latter into symmetric (rotation-symmetry-preserving) and antisymmetric (rotation-symmetry-breaking) strain channels, using a tetragonal f -electron intermetallic DyB 2 C 2 —whose antiferroquadrupolar order breaks local fourfold rotational symmetries while globally remaining tetragonal—as a showcase example. We capture the strain evolution of its quadrupolar and magnetic phase transitions using both singularities in the elastocaloric coefficient and its jumps at the transitions, and the latter we show follows a modified Ehrenfest relation. We find that antisymmetric strain couples to the underlying order parameter in a biquadratic (linear-quadratic) manner in the antiferroquadrupolar (canted antiferromagnetic) phase, which are attributed to a preserved (broken) global tetragonal symmetry, respectively. The broken tetragonal symmetry in the magnetic phase is further evidenced by elastocaloric strain-hysteresis and optical birefringence. Additionally, within the staggered quadrupolar order, the observed elastocaloric response reflects a quadratic increase of entropy with antisymmetric strain, analogous to the role magnetic field plays for Ising antiferromagnetic orders by promoting pseudospin flips. Our results demonstrate AC elastocaloric effect as a compact and incisive thermodynamic probe into the coupling between electronic degrees of freedom and strain in free energy, which holds the potential for investigating and understanding the symmetry of a wide variety of ordered phases in broader classes of quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Space–time reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems

A classical reduced order model for dynamical problems involves spatial reduction of the problem size. However, temporal reduction accompanied by the spatial reduction can further reduce the problem size without losing much accuracy, which results in a considerably more speed-up than the spatial reduction only. Recently, a novel space–time reduced order model for dynamical problems has been developed [17], where the space–time reduced order model shows an order of a hundred speed-up with a relative error of 10 –4 for small academic problems. However, in order for the method to be applicable to a large-scale problem, an efficient space–time reduced basis construction algorithm needs to be developed. Here we present the incremental space–time reduced basis construction algorithm. The incremental algorithm is fully parallel and scalable. Additionally, the block structure in the space–time reduced basis is exploited, which enables the avoidance of constructing the reduced space–time basis. These novel techniques are applied to a large-scale particle transport simulation with million and billion degrees of freedom. The numerical example shows that the algorithm is scalable and practical. Also, it achieves a tremendous speed-up, maintaining a good accuracy. Finally, error bounds for space-only and space–time reduced order models are derived.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗