Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocation”

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 199 records · Page 11

Impact of parametric uncertainty on estimation of the energy deposition into an irradiated brain tumor

Predictions of the total energy deposited into a brain tumor through X-ray irradiation are notoriously error-prone. We investigate how this predictive uncertainty is affected by uncertainty in both the location of the region occupied by a dose-enhancing iodinated contrast agent and the agent's concentration. This is done within the probabilistic framework in which these uncertain parameters are modeled as random variables. We employ the stochastic collocation (SC) method to estimate statistical moments of the deposited energy in terms of statistical moments of the random inputs, and the global sensitivity analysis (GSA) to quantify the relative importance of uncertainty in these parameters on the overall predictive uncertainty. A nonlinear radiation–diffusion equation dramatically magnifies the coefficient of variation of the uncertain parameters, yielding a large coefficient of variation for the predicted energy deposition. This demonstrates that accurate prediction of the energy deposition requires a proper treatment of even small parametric uncertainty. Our analysis also reveals that SC outperforms standard Monte Carlo, but its relative efficiency decreases as the number of uncertain parameters increases from one to three. A robust GSA ameliorates this problem by reducing this number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stellar surface as low-rank modification in iterative methods for binary neutron stars

We present a new multidomain spectral method for the treatment of non-spherical stellar surfaces in iterative methods for binary neutron stars. A stellar surface changes throughout the course of an iterative solution, potentially stalling the convergence. Our method affords low-complexity updates of the relevant subdomain preconditioners, thereby avoiding such stalling. Unlike current collocation (or nodal) approaches for treating surfaces (which rely on coordinate transformations to ensure that stellar surfaces arise at subdomain boundaries), our approach requires no regridding or nontrivial Jacobians. For polytropes with an equation of state specified by an integer polytropic index, our method delivers exponential accuracy with increased truncation, although for “stiff” equations of state (e.g. fractional) it suffers from the same accuracy loss as current methods. We have presented an outline of our approach before, but here present details with numerical tests.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An energy-conserving method for stochastic Maxwell equations with multiplicative noise

In this paper, it is shown that three-dimensional stochastic Maxwell equations with multiplicative noise are stochastic Hamiltonian partial differential equations possessing a geometric structure (i.e. stochastic multi-symplectic conservation law), and the energy of system is a conservative quantity almost surely. We propose a stochastic multi-symplectic energy-conserving method for the equations by using the wavelet collocation method in space and stochastic symplectic method in time. Numerical experiments are performed to verify the excellent abilities of the proposed method in providing accurate solution and preserving energy. The mean square convergence result of the method in temporal direction is tested numerically, and numerical comparisons with finite difference method are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov–Poisson system

The Vlasov–Poisson system, modeling the evolution of non-collisional plasmas in the electrostatic limit, is approximated by a semi-Lagrangian technique. Spectral methods of periodic type are implemented through a collocation approach. Groups of particles are represented by the Fourier Lagrangian basis and evolve, for a single timestep, along a high-order accurate representation of the local characteristic lines. The time-advancing technique is based on truncated Taylor series that can be, in principle, of any order of accuracy. A variant is obtained by coupling the phase space discretization with high-order accurate Backward Differentiation Formulas (BDF). At each timestep, particle displacements are reinterpolated and expressed in the original basis to guarantee the order of accuracy in all the variables at relatively low costs. Thus, these techniques combine the excellent features of spectral approximations with high-order time integration. The resulting method has excellent conservation properties. Indeed, it can be proven that the total number of particles, proportional to the total mass and charge, is conserved up to the machine precision. Series of numerical experiments are performed in order to assess the real performance. In particular, comparisons with standard benchmarks are examined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

What Is the Fractional Laplacian? A Comparative Review with New Results

The fractional Laplacian in $\mathbb{R}^{d}$, which we write as (–Δ) α/2 with, α ϵ (0, 2) has multiple equivalent characterizations. Furthermore, in bounded domains, boundary conditions must be incorporated in these characterizations in mathematically distinct ways, and there is currently no consensus in the literature as to which definition of the fractional Laplacian in bounded domains is most appropriate for a given application. The Riesz (or integral) definition, for example, admits a nonlocal boundary condition, where the value of a function must be prescribed on the entire exterior of the domain in order to compute its fractional Laplacian. Yet, the spectral definition requires only the standard local boundary condition. These differences, among others, lead us to ask the question: “What is the fractional Laplacian?” Beginning from first principles, we compare several commonly used definitions of the fractional Laplacian theoretically, through their stochastic interpretations as well as their analytical properties. Next, we present quantitative comparisons using a sample of state-of-the-art methods. We discuss recent advances on nonzero boundary conditions and present new methods to discretize such boundary value problems: radial basis function collocation (for the Riesz fractional Laplacian) and nonharmonic lifting (for the spectral fractional Laplacian). In our numerical studies, we aim to compare different definitions on bounded domains using a collection of benchmark problems. We consider the fractional Poisson equation with both zero and nonzero boundary conditions, where the fractional Laplacian is defined according to the Riesz definition, the spectral definition, the directional definition, and the horizon-based nonlocal definition. We verify the accuracy of the numerical methods used in the approximations for each operator, and we focus on identifying differences in the boundary behaviors of solutions to equations posed with these different definitions. Through our efforts, we aim to further engage the research community in open problems and assist practitioners in identifying the most appropriate definition and computational approach to use for their mathematical models in addressing anomalous transport in diverse applications.

97 MATHEMATICS AND COMPUTING↗

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↗

Pseudospectral convex optimization for on-ramp merging control of connected vehicles

It can be a daunting task for human drivers to merge into highways because of the intricate vehicle negotiations and potential risk within limited time and space. Connected vehicle (CV) technologies could be a solution to this problem and offer many benefits to the road safety, traffic mobility, and energy efficiency. However, real-time optimal control of CVs is still an open challenge, due to the nonlinear vehicle dynamics, non-convex fuel consumption model, and highly dynamic uncertain inter-vehicle interactions. To tackle these issues, a novel real-time optimal control approach that balances the computational efficiency and solution optimality is proposed for the purpose of onboard application. To this end, the pseudospectral collocation method is integrated with a sequential convex programming approach to develop two new optimization algorithms, which are implemented within a model predictive control (MPC) framework to allow for real-time generation of optimal merging speed profiles. One algorithm leverages the line search technique to improve convergence, and the other benefits from the trust region method for better computational efficiency. The optimality and convergence process of both proposed algorithms are investigated by comparing their solutions with a popular non-linear solver. Furthermore, simulation results show that the proposed methods outperform the benchmark in terms of computational cost, fuel consumption, and traffic efficiency. In particular, the proposed fuel-economy merging rule can save 57.1% fuel consumption on average on four different traffic volumes. Meanwhile, the proposed optimal control algorithms can reduce 2.2% travel time on average comparing to the “first-in-first-out” merging rule.

33 ADVANCED PROPULSION SYSTEMS↗

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗

Performing a multi-unit level-3 PSA with MACCS

MACCS (MELCOR Accident Consequence Code System), WinMACCS, and MelMACCS now facilitate a multi-unit consequence analysis. MACCS evaluates the consequences of an atmospheric release of radioactive gases and aerosols into the atmosphere and is most commonly used to perform probabilistic safety assessments (PSAs) and related consequence analyses for nuclear power plants (NPPs). WinMACCS is a user-friendly preprocessor for MACCS. MelMACCS extracts source-term information from a MELCOR plot file. The current development can combine an arbitrary number of source terms, representing simultaneous releases from a multi-unit facility, into a single consequence analysis. The development supports different release signatures, fission product inventories, and accident initiation times for each unit. The treatment is completely general except that the model is currently limited to collocated units. A major practical consideration for performing a multi-unit PSA is that a comprehensive treatment for more than two units may involve an intractable number of combinations of source terms. This paper proposes and evaluates an approach for reducing the number of calculations to be tractable, even for sites with eight or ten units. The approximation error introduced by the approach is acceptable and is considerably less than other errors and uncertainties inherent in a Level 3 PSA.

42 ENGINEERING↗

Rain footprints on C-band synthetic aperture radar images of the ocean - Revisited

It is well known that rain leaves footprints on the sea surface that sometimes become visible on synthetic aperture radar (SAR) images. Rain cells can easily be detected on SAR images at all radar frequencies when they are associated with a downdraft pattern. But rain cells are not always associated with downdraft and rain can also occur in other forms, as stratified rain, rain bands, and squall lines. It turns out that radar signatures of rain at C-band are much more complex than at L- or X-band radar and that it is particularly difficult to identify unambiguously rain events on C-band SAR images acquired over the ocean. This is because C-band lies in the transition region where raindrops impinging onto the sea surface can increase (usually) or decrease the backscattered radar power and where volume scattering and attenuation by rain drops in the atmosphere are not always negligible (at very high rain rates). In order to get an insight into the physical mechanisms causing the C-band radar signatures of rain, we first revisit results obtained from historic laboratory and field experiments and multi-frequency/multi-polarization SAR data acquired during the SIR-C/X-SAR spaceshuttle mission in 1994. Then we analyze several C-band SAR images acquired by the European satellites Envisat and Sentinel-1A, and the Canadian satellite Radarsat-2 and compare them, whenever possible, with quasi-coincident and collocated weather radar images. The observational data show that, at low to medium rain rates, the main physical mechanism causing C-band radar signatures of rain is Bragg scattering at ring waves generated by the rain drops impinging onto the sea surface, which increase the radar backscatter. However, areas of increased radar backscatter are often accompanied by adjacent areas of decreased radar backscatter, which is due to attenuation of the Bragg waves by turbulence also generated by the impinging rain drops. Furthermore, we present a full-polarimetric Radarsat-2 SAR image of a rain cell together with a polarimetric decomposition analysis, which shows that the C-band radar signature of a rain cell is caused by surface scattering. The observation show that radar signatures of rain cells often contain segments, where the co-polarized as well the cross-polarized radar backscatter are strongly enhanced, which indicates non-Bragg scattering contributions to the scattering process. Furthermore, the polarimetric decomposition analysis shows that the C-band radar signature of a rain cell is dominated by surface scattering. Possible mechanisms, like scattering at splash products, are discussed. Whether the normalized radar cross section (NRCS) due to rain is increased or decreased depends on rain rate, wind speed, incidence angle, and history of the rain event. At low to moderate wind speeds (< 10 ms - 1 ) and low to medium high rain rates (< 50 mm h - 1 ), the NRCS is usually increased by up to 8 dB, and at high wind speeds (> 10 m s - 1 ) and low to high rain rates (but < 50 mm h - 1 ), the NRCS is usually decreased by up to 3 dB.

54 ENVIRONMENTAL SCIENCES↗

Use of physics to improve solar forecast: Part III, impacts of different cloud types

Cloud-type impacts present a great challenge to solar forecasting due to diverse and complex cloud-radiation interactions. This third part of our paper sequence seeks to address this challenge by quantifying the forecast accuracies under eight cloud types: cumulus (Cu), stratified clouds (St), altocumulus (Ac), altostratus (As), cirrostratus/anvil (Cr), cirrus (Ci), congestus (Co), deep convective clouds (Dc) across four physics-informed persistence models reported in Part I. To generalize the cloud impacts, the eight cloud types are further grouped into three cloud categories based on their common features: weak convective clouds, stratiform clouds, and strong convective clouds. Here, the decade-long (2001 ~ 2014) collocated measurements of irradiances and cloud types at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Program South Great Plain (SGP) Central Facility site are used for model evaluation. Results reveal a clear performance hierarchy for global horizontal irradiance (GHI) and direct normal irradiance (DNI): best for weak convective clouds and cirrus, intermediate for stratiform clouds, and worst for strong convective clouds. Performance for diffuse horizontal irradiance (DHI) is less influenced by cloud types. Cloud albedo dominates all three irradiances for Dc, while both cloud albedo and cloud fraction are influential for other cloud types. A 12 %~33 % improvement in accuracy at 6-hour lead time compared to the benchmark smart model confirms the effectiveness of incorporating physics into the models for various cloud types; further improvements are expected by directly integrating cloud type information into forecasting models by modifying the physical formulation of cloud-radiation interaction, and/or using more advanced machine learning models.

14 SOLAR ENERGY↗

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗

Prediction of the Thermal Runaway Limit and Optimal Operation of Heat Transfer-Limited, Fixed-Bed Reactor Systems

Here, we derive a new prediction for thermal runaway starting from the alpha model for fixed-bed reactor systems. This method accounts for thermal resistance internal to the reactor tube and the radial temperature gradients that result. To showcase our method, we compare its predictions to other common criteria for thermal runaway using o-xylene oxidation as the example chemistry. Even in systems where internal heat transfer is negligible, the empirical practical design criterion for thermal runaway is inaccurate. For cases where internal heat transfer is relevant, our runaway limit is more stringent than limits derived from simpler 1-D models. To augment our work, we optimize the product yield with the thermal runaway constraint using orthogonal collocation. Using the alpha model, the results illustrate that the thermal runaway limit can be accurately determined using either numerical or analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Computation of Rotationally-Averaged Pump–Probe X-ray and Electron Diffraction Signals

We develop a new algorithm for the computation of the rotationally averaged elastic molecular diffraction signal for the cases of perpendicular or parallel pump–probe geometries. The algorithm first collocates the charge density from an arbitrary ab initio wave function onto a Becke quadrature grid [A. Becke, J. Chem. Phys. 1988, 88, 2457], providing a high-fidelity multiresolution representation of the charge density. A double sum is then performed over the Becke grid points, and the interaction between points computed using the scattering kernels of Williamson and Zewail [J. C. Williamson and A. H. Zewail, J. Phys. Chem. 1994, 98, 2766]. These kernels analytically average over the molecular orientations with the cos 2 γ selection factor appropriate for one-photon dipole absorption in a perpendicular pump–probe geometry. We demonstrate that the method is converged with small grids containing <500 points/atom. We implement the algorithm on a GPU for increased efficiency and emonstrate the algorithm for molecules with up to a few dozen atoms. We explore the accuracy of the independent atom model (IAM) by comparison with our new and more accurate method. We also investigate the possibility of detecting signatures of electronic transitions in polyatomic pump–probe diffraction experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Contribution of Speciated Monoterpenes to Secondary Aerosol in the Eastern North Atlantic

Extension of laboratory results suggests that monoterpenes, a subset of marine biogenic volatile organic compounds, might play a competitive role in secondary aerosol formation in the remote marine atmosphere, where the oxidation of dimethyl sulfide has traditionally been thought to dominate. However, the current assessment of the role of monoterpenes in secondary aerosol formation in the remote marine atmosphere is limited by the scarcity of speciated monoterpene measurements and the complete absence of collocated measurements of particle chemical composition and monoterpene speciation. Here, in this paper, we present measurements of gas-phase volatile organic compounds, with particular focus on speciated monoterpenes, and commensurate measurements of aerosol chemical composition in the Eastern North Atlantic. The average monoterpene concentration in periods of marine air was 14 ± 10 ppt, and the dominant isomer was β-pinene (46 ± 10%), with other contributions from α-pinene, limonene, and β-ocimene. The total monoterpene concentration in marine air was greater than that of isoprene by a factor of 3.9 on average, in contrast to prior research and potentially due to the large β-pinene contribution. Monoterpenes were significantly smaller in concentration than dimethyl sulfide, which was also highlighted by the substantial contribution of sulfate to non-refractory submicron aerosol. The mean and interquartile range of the sulfate to organic aerosol mass ratio in clean marine air were 2.7 (0.7–3.4). Finally, a simple box model analysis showed that the BVOC precursors, dimethyl sulfide, methanethiol, monoterpenes, and isoprene, could sustain the measured sulfate to organic aerosol ratio in clean marine air.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Control of strong light-matter coupling using the capacitance of metamaterial nanocavities

Metallic nanocavities with deep subwavelength mode volumes can lead to dramatic changes in the behavior of emitters placed in their vicinity. The resulting collocation and interaction often leads to strong coupling. We present for the first time experimental evidence that the Rabi splitting is directly proportional to the electrostatic capacitance associated with the metallic nanocavity. As a result, the system analyzed consists of different metamaterial geometries with the same resonance wavelength coupled to intersubband transitions in quantum wells.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Projecting Changes in Societally Impactful Northeastern U.S. Snowstorms

The northeastern United States is vulnerable to many impacts from snowfall-producing winter cyclones that are amplified by the proximity of population centers to storm tracks. Historically, climatic snowfall assessments have centered around seasonal means even though local impacts typically occur at scales of hours to days. To detect snowstorms at the event level, an objective algorithm is defined based on the Regional Snowfall Index. The metric collocates storm snowfall with population to produce statistics of snowstorms with societal impacts. When applied to the Community Earth System Model Large Ensemble, broad declines in snowstorm frequency are projected by the later 21 st century. These decreases are primarily due to a warmer atmosphere less conducive to snowfall as the predominant precipitation type. However, we note reductions are less significant for major events, since more hostile thermodynamic environments are partially offset by increased precipitation associated with cyclones that dynamically drive high-impact snowstorms.

58 GEOSCIENCES↗