Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Regularization”

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 523 records · Page 29

High resolution dual detector volume-of-interest cone beam breast CT - Demonstration with a bench top system

Purpose: In this study, we used a small field high resolution detector in conjunction with a full field flat panel detector to implement and investigate the dual detector volume-of-interest (VOI) cone beam breast computed tomography (CBCT) technique on a bench-top system. The potential of using this technique to image small calcifications without increasing the overall dose to the breast was demonstrated. Significant reduction of scatter components in the high resolution projection image data of the VOI was also shown. Methods: With the regular flat panel based CBCT technique, exposures were made at 80 kVp to generate an air kerma of 6 mGys at the isocenter. With the dual detector VOI CBCT technique, a high resolution small field CMOS detector was used to scan a cylindrical VOI (2.5 cm in diameter and height, 4.5 cm off-center) with collimated x-rays at four times of regular exposure level. A flat panel detector was used for full field scan with low x-ray exposures at half of the regular exposure level. The low exposure full field image data were used to fill in the truncated space in the VOI scan data and generate a complete projection image set. The Feldkamp-Davis-Kress (FDK) filtered backprojection algorithm was used to reconstruct high resolution images for the VOI. Two scanning techniques, one breast centered and the other VOI centered, were implemented and investigated. Paraffin cylinders with embedded thin aluminum (Al) wires were imaged and used in conjunction with optically stimulated luminescence (OSL) dose measurements to demonstrate the ability of this technique to image small calcifications without increasing the mean glandular dose (MGD). Results: Using exposures that produce an air kerma of 6 mGys at the isocenter, the regular CBCT technique was able to resolve the cross-sections of Al wires as thin as 254 {mu}m in diameter in the phantom. For the specific VOI studied, by increasing the exposure level by a factor of 4 for the VOI scan and reducing the exposure level by a factor of 2 for the full filed scan, the dual-detector CBCT technique was able to resolve the cross-sections of Al wires as thin as 152 {mu}m in diameter. The CNR evaluated for the entire Al wire cross-section was found to be improved from 5.5 in regular CBCT to 14.4 and 16.8 with the breast centered and VOI centered scanning techniques, respectively. Even inside VOI center, the VOI scan resulted in significant dose saving with the dose reduced by a factor of 1.6 at the VOI center. Dose saving outside the VOI was substantial with the dose reduced by a factor of 7.3 and 7.8 at the breast center for the breast centered and VOI centered scans, respectively, when compared to full field scan at the same exposure level. The differences between the two dual detector techniques in terms of dose saving and scatter reduction were small with VOI scan at 4x exposure level and full field scan at 0.5x exposure level. The MGDs were only 94% of that from the regular CBCT scan. Conclusions: For the specific VOI studied, the dual detector VOI CBCT technique has the potential to provide high quality images inside the VOI with MGD similar to or even lower than that of full field breast CBCT. It was also found that our results were compromised by the use of inadequate detectors for the VOI scan. An appropriately selected detector would better optimize the image quality improvement that can be achieved with the VOI CBCT technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MO-DE-207A-02: A Feature-Preserving Image Reconstruction Method for Improved Pancreaticlesion Classification in Diagnostic CT Imaging

Purpose: To develop a feature-preserving model based image reconstruction (MBIR) method that improves performance in pancreatic lesion classification at equal or reduced radiation dose. Methods: A set of pancreatic lesion models was created with both benign and premalignant lesion types. These two classes of lesions are distinguished by their fine internal structures; their delineation is therefore crucial to the task of pancreatic lesion classification. To reduce image noise while preserving the features of the lesions, we developed a MBIR method with curvature-based regularization. The novel regularization encourages formation of smooth surfaces that model both the exterior shape and the internal features of pancreatic lesions. Given that the curvature depends on the unknown image, image reconstruction or denoising becomes a non-convex optimization problem; to address this issue an iterative-reweighting scheme was used to calculate and update the curvature using the image from the previous iteration. Evaluation was carried out with insertion of the lesion models into the pancreas of a patient CT image. Results: Visual inspection was used to compare conventional TV regularization with our curvature-based regularization. Several penalty-strengths were considered for TV regularization, all of which resulted in erasing portions of the septation (thin partition) in a premalignant lesion. At matched noise variance (50% noise reduction in the patient stomach region), the connectivity of the septation was well preserved using the proposed curvature-based method. Conclusion: The curvature-based regularization is able to reduce image noise while simultaneously preserving the lesion features. This method could potentially improve task performance for pancreatic lesion classification at equal or reduced radiation dose. The result is of high significance for longitudinal surveillance studies of patients with pancreatic cysts, which may develop into pancreatic cancer. The Senior Author receives financial support from Siemens GmbH Healthcare.

60 APPLIED LIFE SCIENCES↗

Complement to thermodynamics of dyonic Taub-NUT-AdS spacetime

We examine the thermodynamics of Euclidean dyonic Taub-Nut/Bolt-AdS 4 black holes for a variety of horizon geometries to understand how gauge field regularity conditions influence the thermodynamic relations. We find several distinct features that distinguish the NUT-charged case from its dyonic Reissner-Nordstrom counterpart. For the Nut solution, the gauge field vanishes at the horizon and so regularity is ensured. For the Bolt solution we find that the norm of the gauge field is required to vanish at the horizon in order to satisfy both regularity and the first law of thermodynamics. This regularity condition yields a constraint on the electric and magnetic charges and so reduces cohomogeneity of the system; for spherical horizons, the regularity condition removing the Misner string singularity further reduces cohomogeneity. We observe that bolt solutions with increasing electric charge have positive heat capacity, but upon turning on the magnetic charge to make the solution dyonic, we find that the properties of the uncharged one are retained, having both positive and negative heat capacity. We also study the extremal Bolt solution, finding that Misner string disappears at the horizon in the zero temperature limit. We find that the extremal solution has finite-temperature-like behaviour, with the electric potential playing a role similar to temperature.

79 ASTRONOMY AND ASTROPHYSICS↗

Image Reconstruction from Under sampled Fourier Data Using the Polynomial Annihilation Transform

Fourier samples are collected in a variety of applications including magnetic resonance imaging and synthetic aperture radar. The data are typically under-sampled and noisy. In recent years, l 1 regularization has received considerable attention in designing image reconstruction algorithms from under-sampled and noisy Fourier data. The underlying image is assumed to have some sparsity features, that is, some measurable features of the image have sparse representation. The reconstruction algorithm is typically designed to solve a convex optimization problem, which consists of a fidelity term penalized by one or more l 1 regularization terms. The Split Bregman Algorithm provides a fast explicit solution for the case when TV is used for the l1l1 regularization terms. Due to its numerical efficiency, it has been widely adopted for a variety of applications. A well known drawback in using TV as an l 1 regularization term is that the reconstructed image will tend to default to a piecewise constant image. This issue has been addressed in several ways. Recently, the polynomial annihilation edge detection method was used to generate a higher order sparsifying transform, and was coined the “polynomial annihilation (PA) transform.” This paper adapts the Split Bregman Algorithm for the case when the PA transform is used as the l 1 regularization term. In so doing, we achieve a more accurate image reconstruction method from under-sampled and noisy Fourier data. Our new method compares favorably to the TV Split Bregman Algorithm, as well as to the popular TGV combined with shearlet approach.

97 MATHEMATICS AND COMPUTING↗

Distillation of secret-key from a class of compound memoryless quantum sources

We consider secret-key distillation from tripartite compound classical-quantum-quantum (cqq) sources with free forward public communication under strong security criterion. We design protocols which are universally reliable and secure in this scenario. These are shown to achieve asymptotically optimal rates as long as a certain regularity condition is fulfilled by the set of its generating density matrices. We derive a multi-letter formula which describes the optimal forward secret-key capacity for all compound cqq sources being regular in this sense. We also determine the forward secret-key distillation capacity for situations where the legitimate sending party has perfect knowledge of his/her marginal state deriving from the source statistics. In this case regularity conditions can be dropped. Our results show that the capacities with and without the mentioned kind of state knowledge are equal as long as the source is generated by a regular set of density matrices. We demonstrate that regularity of cqq sources is not only a technical but also an operational issue. For this reason, we give an example of a source which has zero secret-key distillation capacity without sender knowledge, while achieving positive rates is possible if sender marginal knowledge is provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The need for application-based adaptation of deformable image registration

Purpose: To utilize a deformable phantom to objectively evaluate the accuracy of 11 different deformable image registration (DIR) algorithms. Methods: The phantom represents an axial plane of the pelvic anatomy. Urethane plastic serves as the bony anatomy and urethane rubber with three levels of Hounsfield units (HU) is used to represent fat and organs, including the prostate. A plastic insert is placed into the phantom to simulate bladder filling. Nonradiopaque markers reside on the phantom surface. Optical camera images of these markers are used to measure the positions and determine the deformation from the bladder insert. Eleven different DIR algorithms are applied to the full and empty-bladder computed tomography images of the phantom (fixed and moving volumes, respectively) to calculate the deformation. The algorithms include those from MIM Software (MIM) and Velocity Medical Solutions (VEL) and nine different implementations from the deformable image registration and adaptive radiotherapy toolbox for Matlab. These algorithms warp one image to make it similar to another, but must utilize a method for regularization to avoid physically unrealistic deformation scenarios. The mean absolute difference (MAD) between the HUs at the marker locations on one image and the calculated location on the other serves as a metric to evaluate the balance between image similarity and regularization. To demonstrate the effect of regularization on registration accuracy, an additional beta version of MIM was created with a variable smoothness factor that controls the emphasis of the algorithm on regularization. The distance to agreement between the measured and calculated marker deformations is used to compare the overall spatial accuracy of the DIR algorithms. This overall spatial accuracy is also utilized to evaluate the phantom geometry and the ability of the phantom soft-tissue heterogeneity to represent patient data. To evaluate the ability of the DIR algorithms to accurately transfer anatomical contours, the rectum is delineated on both the fixed and moving images. A Dice similarity coefficient is then calculated between the contour on the fixed image and that transferred, via the calculated deformation, from the moving to the fixed image. Results: The phantom possesses sufficient soft-tissue heterogeneity to act as a proxy for patient data. Large discrepancies appear between the algorithms and the measured ground-truth deformation. VEL yields the smallest mean spatial error and a Dice coefficient of 0.90. MIM produces the lowest MAD value and the highest Dice coefficient of 0.96, but creates the largest spatial errors. Increasing the MIM smoothness factor above the default value improves the overall spatial accuracy, but the factor associated with the lowest mean error decreases the Dice coefficient to 0.85. Conclusions: Different applications of DIR require disparate balances between image similarity and regularization. A DIR algorithm that is optimized only for its ability to transfer anatomical contours will yield large deformation errors in homogeneous regions, which is problematic for dose mapping. For this reason, these algorithms must be tested for their overall spatial accuracy. The developed phantom is an objective tool for this purpose.

60 APPLIED LIFE SCIENCES↗

SU-G-BRA-11: Tumor Tracking in An Iterative Volume of Interest Based 4D CBCT Reconstruction

Purpose: 4D CBCT can allow evaluation of tumor motion immediately prior to radiation therapy, but suffers from heavy artifacts that limit its ability to track tumors. Various iterative and compressed sensing reconstructions have been proposed to reduce these artifacts, but are costly time-wise and can degrade the image quality of bony anatomy for alignment with regularization. We have previously proposed an iterative volume of interest (I4D VOI) method which minimizes reconstruction time and maintains image quality of bony anatomy by focusing a 4D reconstruction within a VOI. The purpose of this study is to test the tumor tracking accuracy of this method compared to existing methods. Methods: Long scan (8–10 mins) CBCT data with corresponding RPM data was collected for 12 lung cancer patients. The full data set was sorted into 8 phases and reconstructed using FDK cone beam reconstruction to serve as a gold standard. The data was reduced in way that maintains a normal breathing pattern and used to reconstruct 4D images using FDK, low and high regularization TV minimization (λ=2,10), and the proposed I4D VOI method with PTVs used for the VOI. Tumor trajectories were found using rigid registration within the VOI for each reconstruction and compared to the gold standard. Results: The root mean square error (RMSE) values were 2.70mm for FDK, 2.50mm for low regularization TV, 1.48mm for high regularization TV, and 2.34mm for I4D VOI. Streak artifacts in I4D VOI were reduced compared to FDK and images were less blurred than TV reconstructed images. Conclusion: I4D VOI performed at least as well as existing methods in tumor tracking, with the exception of high regularization TV minimization. These results along with the reconstruction time and outside VOI image quality advantages suggest I4D VOI to be an improvement over existing methods. Funding support provided by CPRIT grant RP110562-P2-01.

60 APPLIED LIFE SCIENCES↗

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗

The sami galaxy survey: revisiting galaxy classification through high-order stellar kinematics

Recent cosmological hydrodynamical simulations suggest that integral field spectroscopy can connect the high-order stellar kinematic moments h{sub 3} (∼skewness) and h{sub 4} (∼kurtosis) in galaxies to their cosmological assembly history. Here, we assess these results by measuring the stellar kinematics on a sample of 315 galaxies, without a morphological selection, using two-dimensional integral field data from the SAMI Galaxy Survey. Proxies for the spin parameter (λ{sub R{sub e}}) and ellipticity (ϵ{sub e}) are used to separate fast and slow rotators; there exists a good correspondence to regular and non-regular rotators, respectively, as also seen in earlier studies. We confirm that regular rotators show a strong h{sub 3} versus V/σ anti-correlation, whereas quasi-regular and non-regular rotators show a more vertical relation in h{sub 3} and V/σ. Motivated by recent cosmological simulations, we develop an alternative approach to kinematically classify galaxies from their individual h{sub 3} versus V/σ signatures. Within the SAMI Galaxy Survey, we identify five classes of high-order stellar kinematic signatures using Gaussian mixture models. Class 1 corresponds to slow rotators, whereas Classes 2–5 correspond to fast rotators. We find that galaxies with similar λ{sub R{sub e}}--ϵ{sub e} values can show distinctly different h{sub 3}--V/σ signatures. Class 5 objects are previously unidentified fast rotators that show a weak h{sub 3} versus V/σ anti-correlation. From simulations, these objects are predicted to be disk-less galaxies formed by gas-poor mergers. From morphological examination, however, there is evidence for large stellar disks. Instead, Class 5 objects are more likely disturbed galaxies, have counter-rotating bulges, or bars in edge-on galaxies. Finally, we interpret the strong anti-correlation in h{sub 3} versus V/σ as evidence for disks in most fast rotators, suggesting a dearth of gas-poor mergers among fast rotators.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

Multilevel-in-width training for deep neural network regression

A common challenge in regression is that for many problems, the degrees of freedom required for a high-quality solution also allows for overfitting. Regularization is a class of strategies that seek to restrict the range of possible solutions so as to discourage overfitting while still enabling good solutions, and different regularization strategies impose different types of restrictions. In this paper, we present a multilevel regularization strategy that constructs and trains a hierarchy of neural networks, each of which has layers that are wider versions of the previous network's layers. Here, we draw intuition and techniques from the field of Algebraic Multigrid (AMG), traditionally used for solving linear and nonlinear systems of equations, and specifically adapt the Full Approximation Scheme (FAS) for nonlinear systems of equations to the problem of deep learning. Training through V-cycles then encourage the neural networks to build a hierarchical understanding of the problem. We refer to this approach as multilevel-in-width to distinguish from prior multilevel works which hierarchically alter the depth of neural networks. The resulting approach is a highly flexible framework that can be applied to a variety of layer types, which we demonstrate with both fully connected and convolutional layers. We experimentally show with PDE regression problems that our multilevel training approach is an effective regularizer, improving the generalize performance of the neural networks studied.

97 MATHEMATICS AND COMPUTING↗

Warped $AdS_6 × S^2$ in Type IIB supergravity I: local solutions

We investigate the existence of solutions with 16 residual supersymmetries to Type IIB supergravity on a space-time of the form $AdS_6×S^2$ warped over a two-dimensional Riemann surface $Σ$. The $SO(2,5)×SO(3)$ isometry extends to invariance under the exceptional Lie superalgebra $F$(4). In the present paper, we construct the general Ansatz compatible with these symmetries, derive the corresponding reduced BPS equations, and obtain their complete local solution in terms of two locally holomorphic functions $A_±$ on $Σ$, subject to certain positivity and regularity conditions. Globally, ($A_+,A_-$) are allowed to be multiple-valued on $Σ$ and be holomorphic sections of a holomorphic bundle over $Σ$ with structure group contained in $SU$(1,1)×C. Globally regular solutions are expected to provide the near-horizon geometry of (p,q) 5-brane and 7-brane webs which are holographic duals to five-dimensional conformal field theories. A preliminary analysis of the positivity and regularity conditions will be presented here, leaving the construction of globally regular solutions to a subsequent paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

Factorization of generalized γ-generating matrices

The class of γ-generating matrices and its subclasses of regular and singular γ-generating matrices were introduced by D. Z. Arov in [8], where it was shown that every γ-generating matrix admits an essentially unique regular–singular factorization. The class of generalized γ-generating matrices was introduced in [14, 20]. In the present paper, subclasses of singular and regular generalized –generating matrices are introduced and studied. As the main result, a theorem of existence of the regular–singular factorization for a rational generalized γ-generating matrix is proved.

97 MATHEMATICS AND COMPUTING↗

Net energy ratio for the production of steam pretreated biomass-based pellets

In this study, a process model was developed to determine the net energy ratio (NER) for both regular and steam-pretreated pellet production from ligno-cellulosic biomass. NER is a ratio of the net energy output to the total net energy input from non-renewable energy source into the system. Scenarios were developed to measure the effect of temperature and level of steam pretreatment on the NER of both production processes. The NER for the base case at 6 kg h –1 is 1.29 and 5.0 for steam-pretreated and regular pellet production respectively. However, at the large scale NER would improve. The major factor for NER is energy for steam and drying unit. The sensitivity analysis for the model shows that the optimum temperature for steam pretreatment is 200 °C with 50% pretreatment (Steam pretreating 50% feed stock, while the rest is undergoing regular pelletization). Uncertainty result for steam pretreated and regular pellet is 1.35 ± 0.09 and 4.52 ± 0.34 respectively.

09 BIOMASS FUELS↗

A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

02 PETROLEUM↗

Single-walled carbon nanotubes (SWCNTs) inhibit heat shock protein 90 (HSP90) signaling in human lung fibroblasts and keratinocytes

Single-walled carbon nanotubes (SWCNTs) are carbon-based nanomaterials that possess immense industrial potential. Despite accumulating evidence that exposure to SWCNTs might be toxic to humans, our understanding of the mechanisms for cellular toxicity of SWCNTs remain limited. Here, we demonstrated that acute exposure of short (1–3 μm) and regular-length (5–30 μm) pristine, carboxylated or hydroxylated SWCNTs inhibited cell proliferation in human somatic and human stem cells in a cell type-dependent manner. The toxicity of regular-length pristine SWCNT was most evidenced in NP69 > CYT00086 > MCF-10A > MRC-5 > HaCaT > HEK-293 T > HepG2. In contrast, the short pristine SWCNTs were relatively less toxic in most of the cells being tested, except for NP69 which is more sensitive to short pristine SWCNTs as compared to regular-length pristine SWCNTs. Interestingly, carboxylation and hydroxylation of regular-length SWCNTs, but not the short SWCNTs, significantly reduced the cytotoxicity. Exposure of SWCNTs also induced caspase 3 and 9 activities, mitochondrial membrane depolarization, and significant apoptosis and necrosis in MRC-5 embryonic lung fibroblasts. In contrast, SWCNTs inhibited the proliferation of HaCaT human keratinocytes without inducing cell death. Further analyses by gene expression profiling and Connectivity Map analysis showed that SWCNTs induced a gene expression signature characteristic of heat shock protein 90 (HSP90) inhibition in MRC-5 cells, suggesting that SWCNTs may inhibit the HSP90 signaling pathway. Indeed, exposure of MRC-5 cells to SWCNTs results in a dose-dependent decrease in HSP90 client proteins (AKT, CDK4 and BCL2) and a concomitant increase in HSP70 expression. In addition, SWCNTs also significantly inhibited HSP90-dependent protein refolding. Finally, we showed that ectopic expression of HSP90, but not HSP40 or HSP70, completely abrogated the cytotoxic effects of SWCNTs, suggesting that SWCNT-induced cellular toxicity is HSP90 dependent. In summary, our findings suggest that the toxic effects of SWCNTs are mediated through inhibition of HSP90 in human lung fibroblasts and keratinocytes. - Highlights: • SWCNTs inhibited cell proliferation in human somatic and human stem cells. • SWCNTs induced cytotoxicity and growth inhibition in MRC-5 and HaCat cells respectively. • SWCNTs induced a gene signature characteristic of HSP90 inhibition. • SWCNT-induced toxicity in MRC-5 and HaCat is HSP90-dependent.

60 APPLIED LIFE SCIENCES↗