Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Integrated High-Speed Torque Control System for a Robotic Joint

A control system for achieving high-speed torque for a joint of a robot includes a printed circuit board assembly (PCBA) having a collocated joint processor and high-speed communication bus. The PCBA may also include a power inverter module (PIM) and local sensor conditioning electronics (SCE) for processing sensor data from one or more motor position sensors. Torque control of a motor of the joint is provided via the PCBA as a high-speed torque loop. Each joint processor may be embedded within or collocated with the robotic joint being controlled. Collocation of the joint processor, PIM, and high-speed bus may increase noise immunity of the control system, and the localized processing of sensor data from the joint motor at the joint level may minimize bus cabling to and from each control node. The joint processor may include a field programmable gate array (FPGA).

Davis, Donald R.↗
NASA NTRS2016

Parallel Aircraft Trajectory Optimization with Analytic Derivatives

Trajectory optimization is an integral component for the design of aerospace vehicles, but emerging aircraft technologies have introduced new demands on trajectory analysis that current tools are not well suited to address. Designing aircraft with technologies such as hybrid electric propulsion and morphing wings requires consideration of the operational behavior as well as the physical design characteristics of the aircraft. The addition of operational variables can dramatically increase the number of design variables which motivates the use of gradient based optimization with analytic derivatives to solve the larger optimization problems. In this work we develop an aircraft trajectory analysis tool using a Legendre-Gauss-Lobatto based collocation scheme, providing analytic derivatives via the OpenMDAO multidisciplinary optimization framework. This collocation method uses an implicit time integration scheme that provides a high degree of sparsity and thus several potential options for parallelization. The performance of the new implementation was investigated via a series of single and multi-trajectory optimizations using a combination of parallel computing and constraint aggregation. The computational performance results show that in order to take full advantage of the sparsity in the problem it is vital to parallelize both the non-linear analysis evaluations and the derivative computations themselves. The constraint aggregation results showed a significant numerical challenge due to difficulty in achieving tight convergence tolerances. Overall, the results demonstrate the value of applying analytic derivatives to trajectory optimization problems and lay the foundation for future application of this collocation based method to the design of aircraft with where operational scheduling of technologies is key to achieving good performance.

aircraft↗
NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗
NASA NTRS2019

Low-Thrust Trajectory Design for a Cislunar CubeSat Leveraging Structures from the Bicircular Restricted Four-Body Problem

The upcoming Lunar IceCube (LIC) mission will deliver a 6U CubeSat to a low lunar orbit via a ride-share opportunity during NASAs Artemis 1 mission. This presents a challenging trajectory design scenario, as the vast change in energy required to transfer from the initial deployment state to the destination orbit is compounded by the limitations of the LICs low-thrust engine. This investigation addresses these challenges by developing a trajectory design framework that utilizes dynamical structures available in the Bicircular Restricted Four-Body Problem (BCR4BP) along with a robust direct collocation algorithm. Maps are created that expedite the selection of invariant manifold paths from a periodic staging orbit in the BCR4BP that offer favorable connections between the LIC transfer phases. Initial guesses assembled from these maps are passed to a direct collocation algorithm that corrects them in the BCR4BP while including the variable low-thrust acceleration of the spacecraft engine. Results indicate that the ordered motion provided by the BCR4BP and the robustness of direct collocation combine to offer an efficient and adaptable framework for designing a baseline trajectory for the LIC mission.

Pritchett, Robert↗
NASA NTRS2019

Low-Thrust Trajectory Design for A Cislunar CubeSat Leveraging Structures from the Bicircular Restricted Four-Body Problem

The upcoming Lunar IceCube (LIC) mission will deliver a 6U CubeSat to a low lunar orbit via a ride-share opportunity during NASAs Artemis 1 mission. This presents a challenging trajectory design scenario, as the vast change in energy required to transfer from the initial deployment state to the destination orbit is compounded by the limitations of the LICs low-thrust engine. This investigation addresses these challenges by developing a trajectory design framework that utilizes dynamical structures available in the Bicircular Restricted Four-Body Problem (BCR4BP) along with a robust direct collocation algorithm. Maps are created that expedite the selection of invariant manifold paths from a periodic staging orbit in the BCR4BP that offer favorable connections between the LIC transfer phases. Initial guesses assembled from these maps are passed to a direct collocation algorithm that corrects them in the BCR4BP while including the variable low-thrust acceleration of the spacecraft engine. Results indicate that the ordered motion provided by the BCR4BP and the robustness of direct collocation combine to offer an efficient and adaptable framework for designing a baseline trajectory for the LIC mission.

Pritchett, Robert↗
NASA NTRS2020

A Machine Learning-Based Cloud Detection and Thermodynamic Phase Classification Algorithm using Passive Spectral Observations

We trained two Random Forest (RF) machine-learning models for cloud mask and cloud thermodynamic phase detection using spectral observations from VIIRS on Suomi NPP (SNPP). Observations from CALIOP were carefully selected to provide reference labels. The two RF models were trained for all-day and daytime-only conditions using a 4-year collocated VIIRS/CALIOP dataset from 2013 to 2016. Due to the orbit difference, the collocated CALIOP and SNPP VIIRS training samples cover a broad viewing zenith angle range, which is a great benefit to overall model performance. The all-day model uses 3 VIIRS infrared (IR) bands (8.6,11, and 12 μm) and the daytime model uses 5 Near-IR (NIR) and Shortwave-IR (SWIR) bands (0.86, 1.24, 1.38, 1.64 and 2.25 μm) together with the 3 IR bands to detect clear, liquid water, and ice cloud pixels. Up to 7 surface types, namely, ocean/water, forest, cropland, grassland, snow/ice, barren/desert, and shrubland, were considered separately to enhance performance for both models. Detection of cloudy pixels and thermodynamic phase with the two RF models were compared against collocated CALIOP products from 2017. It is shown that, with a conservative screening process that excludes the most challenging cloudy pixels for passive remote sensing, the two RF models have high accuracy rates in comparison with the CALIOP reference for both cloud detection and thermodynamic phase. Other existing SNPP VIIRS and Aqua MODIS cloud mask and phase products are also evaluated, with results showing that the two RF models and the MODIS MYD06 optical property phase product are the top 3 algorithms with respect to lidar observations during the daytime. During the nighttime, the RF all-day model works best for both cloud detection and phase, in particular for pixels over snow/ice surfaces. The present RF models can be extended to other similar passive instruments if training samples can be collected from CALIOP or other lidars. However, the quality of reference labels and potential sampling issues that may impact model performance would need further attention.

cloud detection↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

A defining feature of the discontinuous Galerkin (DG) method for ODE is that the piecewise polynomial solution can have a jump discontinuity at the beginning of each step. Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠 − 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK-DG, collocation, and CG schemes.

numerical methods↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗
NASA NTRSDate not supplied

Creating Seasonal Climatologies of Aerosol Lidar Ratios Over Ocean Using MODIS AOD Constrained Retrievals and GOCART Model Simulations

The current CALIPSO algorithms assign one lidar ratio (i.e., extinction-to-backscatter ratio; LR) value globally for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms through the development of regional and seasonal LR climatologies. In this work, aerosol LRs are inferred through CALIOP backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type (e.g., marine). The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. In this poster, we show twelve-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine aerosols and how they correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). Near land masses, smaller SSVFs and larger LRs are found (due to the influence of over-land aerosols). In the remote ocean regions (likely less impacted by over-land aerosols), larger SSVFs and smaller LRs are found. The developed relationship between the GOCART model SSVFs and MODIS AOD constrained LRs is used to create model-assisted seasonal LR maps. Additionally, we show maps of inferred LRs from constrained retrievals using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) product and comparisons with those from the MODIS analyses. The technique demonstrated here benefits CALIPSO in the near-term, but similar methods can also be applied to the next generation space-based elastic backscatter lidars with collocated passive sensors, such as those of the upcoming NASA Atmosphere Observing System (AOS).

Travis Toth↗
NASA NTRSDate not supplied

Mapping Aerosol Lidar Ratios Over Ocean Using Constrained Retrievals and A Global Aerosol Model

The current NASA Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) algorithms assign one lidar ratio (i.e., extinction-to-backscatter ratio; LR) value globally for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms through the development of regional and seasonal LR climatologies. In this work, aerosol LRs are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type (e.g., marine). The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. In this presentation, we show that the twelve-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine aerosols correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). For example, smaller SSVFs (< 65%) and larger LRs (> 55 sr), are found near land masses. This is indicative of the influence of over-land aerosols (e.g., pollution and biomass burning smoke). In the remote oceans (i.e., regions likely less impacted by non-sea salt aerosols), larger SSVFs (> 95%) and smaller LRs (< 25 sr) are found. The developed relationship between the GOCART SSVFs and MODIS AOD constrained LRs (polynomial fit intersect values of ~56 sr for SSVF of 0% and ~21 sr for SSVF of 100%) is used to produce model-assisted climatological LR maps on seasonal scales. Additionally, we show maps of inferred LRs from constrained retrievals using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) product and comparisons with those from the MODIS analyses. The technique demonstrated in this study not only benefits CALIPSO in the near-term, but similar methods can be applied to future spaceborne elastic backscatter lidars with collocated passive sensors, such as those associated with the upcoming NASA Atmosphere Observing System (AOS).

Travis Toth↗
NASA NTRSDate not supplied

Mapping Aerosol Lidar Ratios Over Ocean using MODIS AOD Constrained Retrievals and GOCART Model Simulations

After 17 years, the NASA Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission ceased science operations in August 2023. For the final CALIPSO data products release (Version 5), the CALIPSO project seeks to improve the accuracy of its aerosol extinction by advancing knowledge of aerosol lidar ratios (i.e., extinction-to-backscatter ratios; LRs) for various aerosol types. The current algorithm assigns one LR value globally for each of the seven tropospheric aerosol types. The CALIPSO team aims to improve the retrieval algorithm through the development of regional and seasonal LR climatologies for the same aerosol types. In this study, aerosol LRs are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth (AOD) from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data over oceans during daytime. This analysis is subsampled for those profiles that are cloud-free and contain only one CALIOP-classified aerosol type. The CALIOP profiles are then collocated with aerosol volume fractions obtained through Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulations. This presentation will reveal findings that the 12-year (June 2006-August 2018) mean spatial distributions of inferred aerosol LRs for CALIOP-classified marine and dusty marine aerosols correspond inversely with patterns of GOCART sea salt volume fraction (SSVF). For example, smaller SSVFs (< 65%) and larger LRs (> 55 sr), are found near land masses (Fig. 1). This indicates the influence of advected anthropogenic aerosols (e.g., pollution and biomass burning smoke). In the remote oceans (i.e., regions likely less impacted by non-sea salt aerosols), the SSVFs are larger (> 95%) and the LRs are smaller (< 25 sr) (Fig. 1). A polynomial fit of the MODIS AOD constrained LRs to the corresponding GOCART SSVFs (intersect values of ~58 sr for SSVF of 0% and ~21 sr for SSVF of 100%) is further used to produce model-assisted climatological LR maps on seasonal scales. Additionally, we will show results of a LR validation analysis for which we compare the revised CALIPSO AODs obtained by applying the seasonal/regional constrained LRs against CALIPSO Version 4.51 Ocean Derived Column Optical Depth (ODCOD). While the majority of the presentation will focus on LRs for CALIOP-classified marine and dusty marine aerosols, an overview of LR results will show preliminary results for other aerosol types over ocean, such as dust and elevated smoke. The technique demon-strated in this study highlights the benefits not only to the final planned CALIPSO data release in 2025, but similar methods can be applied to future spaceborne elastic backscatter lidars with collocated passive sensors (e.g., such as those associated with NASA’s proposed Atmosphere Observing System).

Travis D Toth↗
DOE OSTI2016

Radiation Transport in Random Media With Large Fluctuations

Strongly heterogeneous media arise in several applications that include radiation shields, nuclear fuel, BWR moderators, clouds, planetary and stellar atmospheres, turbulent gases and plasmas. Neutral and charged particle transport computations in such media have relied heavily on formulating transport equations with spatially random coefficients (physical data) and developing solution methods to deal with the additional stochastic dimensions. Attempts at developing approximate closures that yield only low order statistical information (e.g., mean and variance of the flux) have proved to be highly restrictive under real physics conditions or rely on techniques that require fluctuation amplitudes to be small for robustness. Recently, stochastic spectral methods such as polynomial chaos and stochastic collocation have been developed for aleatoric uncertainty quantification and sensitivity analysis, and successfully applied in radiation transport work. Advances in UQ techniques have tended to focus on efficiently handling large numbers of uncertain variables but the rigorous stochastic basis of the approach also promotes its use in situations where stochasticity is due to spatial heterogeneity and the associated uncertainty is large. Here we apply these techniques to radiation transport in media with spatially randomly varying cross sections without restriction on fluctuation amplitudes. Specifically, we represent the cross sections as a lognormal spatial random process with specified mean, variance and covariance function and use a Karhunen-Loeve (KL) decomposition to generate cross section realizations that are strictly positive. Woodcock Monte Carlo (WMC) is then used to simulate transport using random sampling of cross sections and deterministic sampling based on a stochastic collocation technique. Numerical results for the mean and variance in the scalar flux and leakage currents are obtained for weak and strong random variations. We focus in this work on demonstrating the approach and defer application to specific problems to a future investigation. We have demonstrated that a combination of stochastic spectral representation and Woodcock Monte Carlo simulation allows efficient computation of radiation transport in strongly random media. The use of a lognormal distribution to represent the fluctuations in medium properties facilitates enforcement of strict positivity of cross section realizations, a condition that has challenged standard approaches that assume Gaussian fluctuations. The fact that the lognormal random process is a memoryless nonlinear transformation of a Gaussian process further enables efficient and accurate reconstruction of the random cross section realizations using a Karhunen-Loeve representation of the cross section. Finally, the robustness of this model enables fluctuations of arbitrarily large amplitude to be studied, using both random sampling Monte Carlo and more efficient stochastic collocation techniques. The present work is being extended to allow reconstruction of probability density functions of output variables using polynomial chaos representation, and to accommodate general cross section covariance functions. Also, techniques are being explored to more judiciously select ceiling cross sections to enhance the efficiency of Woodcock Monte Carlo for transport in random media. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Evaluating Grid Strength under Uncertain Renewable Generation

The increasing displacement of synchronous generators with renewable resources such as wind and solar via power electronic interfaces causes a reduction in short-circuit strength and weak grid issues. The variation and uncertainty of renewable energy increase challenges for identifying weak grid conditions. This paper proposes an efficient method to analyze the impact of uncertain renewable energy on grid strength. The proposed method uses the probabilistic collocation method (PCM) to approximate the results of grid strength assessment under uncertain renewable generation, in order to reduce computational burden without compromising result accuracy when compared with traditional Monte Carlo simulation (MCS). To improve the accuracy of the approximation results, the proposed method integrates the K-means clustering technique with PCM to select the approximation samples of input variables. The efficacy of the proposed method is demonstrated by comparison with MCS on the modified IEEE 9-bus system and modified IEEE 39-bus system with multiple renewable generators.

grid strength↗
DOE OSTI2015

Seismicity on the western Greenland Ice Sheet: Surface fracture in the vicinity of active moulins

We analyzed geophone and GPS measurements collected within the ablation zone of the western Greenland Ice Sheet during a ~35 day period of the 2011 melt season to study changes in ice deformation before, during, and after a supraglacial lake drainage event. During rapid lake drainage, ice flow speeds increased to ~400% of winter values, and icequake activity peaked. At times >7 days after drainage, this seismicity developed variability over both diurnal and longer periods (~10 days), while coincident ice speeds fell to ~150% of winter values and showed nightly peaks in spatial variability. Approximately 95% of all detected seismicity in the lake basin and its immediate vicinity was triggered by fracture propagation within near-surface ice (<330 m deep) that generated Rayleigh waves. Icequakes occurring before and during drainage frequently were collocated with the down flow (west) end of the primary hydrofracture through which the lake drained but shifted farther west and outside the lake basin after the drainage. We interpret these results to reveal vertical hydrofracture opening and local uplift during the drainage, followed by enhanced seismicity and ice flow on the downstream side of the lake basin. This region collocates with interferometric synthetic aperture radar-measured speedup in previous years and could reflect the migration path of the meltwater supplied to the bed by the lake. The diurnal seismic signal can be associated with nightly reductions in surface melt input that increase effective basal pressure and traction, thereby promoting elevated strain in the surficial ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

How does increasing horizontal resolution in a global climate model improve the simulation of aerosol-cloud interactions?

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25 degrees, with the meteorology nudged towards the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the northern hemisphere mid-latitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m⁻² (-30%) compared to the 2° model, leading to a 0.26 W m⁻² reduction (-15%) in the global annual average AIF.

Ma, Po -Lun↗
DOE OSTI2015

Chemical Imaging of Ambient Aerosol Particles: Observational Constraints on Mixing State Parameterization

A new parameterization for quantifying the mixing state of aerosol populations has been applied for the first time to samples of ambient particles analyzed using spectro-microscopy techniques. Scanning transmission x-ray microscopy/near edge x-ray absorption fine structure (STXM/NEXAFS) and computer controlled scanning electron microscopy/energy dispersive x-ray spectroscopy (CCSEM/EDX) were used to probe the composition of the organic and inorganic fraction of individual particles collected on June 27th and 28th during the 2010 Carbonaceous Aerosols and Radiative Effects (CARES) study in the Central Valley, California. The first field site, T0, was located in downtown Sacramento, while T1 was located near the Sierra Nevada Mountains. Mass estimates of the aerosol particle components were used to calculate mixing state metrics, such as the particle-specific diversity, bulk population diversity, and mixing state index, for each sample. Both microscopy imaging techniques showed more changes over these two days in the mixing state at the T0 site than at the T1 site. The STXM data showed evidence of changes in the mixing state associated with a build-up of organic matter confirmed by collocated measurements and the largest impact on the mixing state was due to an increase in soot dominant particles during this build-up. The CCSEM/EDX analysis showed the presence of two types of particle populations; the first was dominated by aged sea salt particles and had a higher mixing state index (indicating a more homogeneous population), the second was dominated by carbonaceous particles and had a lower mixing state index.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

Superconvergence results of Legendre spectral projection methods for Volterra integral equations of second kind

In this paper, Legendre spectral projection methods are applied for the Volterra integral equations of second kind with a smooth kernel. We prove that the approximate solutions of the Legendre Galerkin and Legendre collocation methods converge to the exact solution with the order O(n{sup −r}) in L{sup 2}-norm and order O(n{sup −r+(1/2)}) in infinity norm, and the iterated Legendre Galerkin solution converges with the order O(n{sup −2r}) in both L{sup 2}-norm and infinity norm, whereas the iterated Legendre collocation solution converges with the order O(n{sup −r}) in both L{sup 2}-norm and infinity norm, n being the highest degree of Legendre polynomials employed in the approximation and r being the smoothness of the kernels. We have also considered multi-Galerkin method and its iterated version, and prove that the iterated multi-Galerkin solution converges with the order O(n{sup −3r}) in both infinity and L{sup 2} norm. Numerical examples are given to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗