Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pharmaceuticals”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of the cyanobacterial amino acid metabolism with a precise genome-scale metabolic reconstruction of Anabaena sp. UTEX 2576

Filamentous cyanobacteria such as Anabaena sp. UTEX 2576 (a.k.a., Nostoc sp. PCC 7120) are a sustainable platform for commodity and specialty chemical production. Anabaena is a model cyanobacterium adopted to study production of nitrogen-containing metabolites useful for chemical, cosmetic, and pharmaceutical industries. Therefore, a precise description of the Anabaena metabolic network is desired to analyze their amino-acid metabolism and elucidate potential applications of these cyanobacteria as host organisms for biotechnological production. Secondary metabolite production depends on the metabolic availability of amino acids. We provide new insight on the biotechnological utilization of Anabaena after predicting phycocyanobilin and amino-acid production rates, using a genome-scale metabolic model (iDN1004). This metabolic reconstruction is a highly comprehensive representation of the global metabolism of Anabaena, which also contains experimental biomass equations and constraints under photoautotrophic and photodiazotrophic growth. Modelling results ranked proteinogenic amino acids based on predicted metabolic fluxes through amino-acid producing reactions. From these, l-aspartate, glycine, l-serine, l-valine, l-alanine, l-threonine, and l-leucine were selected as the best branching points to conduct metabolic engineering, where l-aspartate, l-serine, l-valine and glycine serve as potential precursors of the secondary metabolites Schizokinen (a siderophore), Sphingosine (a ceramide), Lyngbyatoxin A (an alkaloid), and Shinorine (a sunscreen).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High-level co-production of 3-hydroxypropionic acid and 1,3-propanediol from glycerol: Metabolic engineering and process optimization

3-Hydroxypropionic acid (3-HP) and 1,3-propanediol (1,3-PDO) are value-added chemicals with versatile applications in the chemical, pharmaceutical, and food industries. Nevertheless, sustainable production of 3-HP and 1,3-PDO is often limited by the lack of efficient strains and suitable fermentation configurations. We report attempts have been made to improve the co-production of both metabolites through metabolic engineering of Escherichia coli and process optimization. First, the 3-HP and 1,3-PDO co-biosynthetic pathways were recruited and optimized in E. coli, followed by coupling the pathways to the transhydrogenase-mediated cofactor regeneration systems that increased cofactor availability and product synthesis. Next, pathway rebalancing and block of by-product formation significantly improved 3-HP and 1,3-PDO net titer. Subsequently, glycerol flux toward 3-HP and 1,3-PDO synthesis was maximized by removing metabolic repression and fine-tuning the glycerol oxidation pathway. Lastly, the combined fermentation process optimization and two-stage pH-controlled fed-batch fermentation co-produced 140.50 g/L 3-HP and 1,3-PDO, with 0.85 mol/mol net yield.

1,3-propanediol↗
DOE OSTI2026

Functionalized benzylamines from commercial kraft lignin

Benzylamines are key intermediates in pharmaceuticals, agrochemicals, and polymers, but their conventional production relies on benzyl chloride - a petroleum-derived compound with high toxicity and energy demands. Lignin, accounting for up to 30% of plant biomass, is the largest renewable source of aromatic carbon on Earth. However, its highly complex and recalcitrant structure poses a major barrier to efficient conversion into high-value chemicals. Here, in this study, we developed a catalytic approach to convert commercial kraft lignin into phenolic benzylamines through selective depolymerization and subsequent functionalization. We systematically evaluated the effects of three alcohol solvents, formic acid (FA), and a ruthenium-on‑carbon (Ru/C) catalyst on monophenol yield and selectivity. Up to 6.5 wt% monophenol yield was achieved using methanol (MeOH), FA, and Ru/C at 300 °C for 2 h. Quantum thermodynamic simulations based on the COSMO-RS model confirmed the superior solvation and reactivity of the MeOH + FA system, rationalizing observed product yield. The purified monophenolic products, primarily guaiacol and alkyl guaiacols, were then converted into functionalized benzylamines with >90% yield via a multicomponent Mannich reaction under mild conditions. Techno economic analysis (TEA) and life cycle assessment (LCA) underscore the importance of improving lignin depolymerization yields and expanding biorefinery scale. Solvent-only configurations outperform other options in both cost and emissions, with the methanol-only case performing the best ($\$$105 /kg and 26 kg CO 2e /kg) at a large-scale facility. This study establishes a scalable, bio-based pathway for producing benzylamines from commercial kraft lignin, advancing lignin valorization and offering a sustainable alternative to produce petrochemical-based benzylamines.

Catalytic reduction↗
DOE OSTI2022

Rare earth element separations by high-speed counter-current chromatography

We report that following the initial development of High-Speed Counter-Current Chromatography (HSCCC) in the 1960s, several studies have explored its applicability in the separation of rare earth elements (REEs). More recently, however, HSCCC publications have transitioned towards the separation of natural products or pharmaceuticals, leaving the application for REEs largely unexplored from a practical standpoint. Herein, we expand upon prior work in this field by evaluating the suitability of HSCCC to separation of a subset of non-radioactive REEs (Nd, Sm, Eu, Tb, and Y) at 10 -4 mol levels using di-(2-ethylhexyl)phosphoric acid (HDEHP) in n-heptane as the stationary phase and hydrochloric acid as the mobile phase. First, the effect of flow rate on the stationary phase volume retention ratio and resolution of Nd/Sm/Eu subgroup was evaluated followed by optimization of step-gradient elution profiles resulting in additional recovery of Tb and Y within a seven-hour window. The five REEs were separated at the baseline resolution level or above. Elution profiles obtained from multiple runs across two independently operated columns and across independent runs were cross analyzed. Reproducibility in elution profiles point to future applications in radioelement separation chemistry, where both chemical and radiochemical purity are of importance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Use of Fisher's Ratio assisted multivariate curve resolution- alternating least squares for discovery-based analysis using ultrahigh pressure liquid chromatography-high resolution mass spectrometry

Non-targeted analysis of complex chemical mixtures can be difficult considering the convoluted nature of the matrix and the potential unknown chemical differences between samples or classes of samples. Ultrahigh pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOF) is an ideal technique to probe chemical differences for a wide variety of samples. While UHPLC-QTOF can discover minute chemical differences down to low part per billion (ppb) concentrations with a high degree of confidence, the application of high-resolution mass spectrometry can yield massive amounts of information (∼ 10 gb per sample) that cannot be analyzed manually. Therefore, the application of chemometric techniques is mandatory for the interrogation of complex samples. Fisher's ratio (FR) assisted multivariate curve resolution-alternating least squares (MCR-ALS) was used to the discover and identify the chemical differences between two classes of materials: 1) a pond water matrix and 2) the matrix spiked with a pharmaceutical standard mix containing 17 compounds. Thirteen of the seventeen spiked compounds were discovered using FR analysis, and then five were successfully deconvoluted using MCR-ALS wherein the number of curves chosen were automatically determined using singular value decomposition (SVD). In conclusion, the use of an automated FR assisted MCR-ALS will aid in discovering trace levels of chemical components without the need for the researcher to provide potentially biased input which will aid in non-targeted workflow.

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

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Overcoming barriers to medium-chain fatty alcohol production

Medium-chain fatty alcohols ( mc FaOHs) are aliphatic primary alcohols containing six to twelve carbons that are widely used in materials, pharmaceuticals, and cosmetics. Microbial biosynthesis has been touted as a route to less-abundant chain-length molecules and as a sustainable alternative to current petrochemical processes. Several metabolic engineering strategies for producing mc FaOHs have been demonstrated in the literature, yet processes continue to suffer from poor selectivity and mc FaOH toxicity, leading to reduced titers, rates, and yields of the desired compounds. Here, this opinion examines the current state of microbial mc FaOH biosynthesis, summarizing engineering efforts to tailor selectivity and improve product tolerance by implementing engineering strategies that circumvent or overcome mc FaOH toxicity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

High-throughput screening tools facilitate calculation of a combined exposure-bioactivity index for chemicals with endocrine activity

High-throughput and computational tools provide a new opportunity to calculate combined bioactivity of exposure to diverse chemicals acting through a common mechanism. We used high throughput in vitro bioactivity data and exposure predictions from the U.S. EPA’s Toxicity and Exposure Forecaster (ToxCast and ExpoCast) to estimate combined estrogen receptor (ER) agonist activity of non-pharmaceutical chemical exposures for the general U.S. population. High-throughput toxicokinetic (HTTK) data provide conversion factors that relate bioactive concentrations measured in vitro (µM), to predicted population geometric mean exposure rates (mg/kg/day). These data were available for 22 chemicals with ER agonist activity and were estimated for other ER bioactive chemicals based on the geometric mean of HTTK values across chemicals. For each chemical, ER bioactivity across ToxCast assays was compared to predicted population geometric mean exposure at different levels of in vitro potency and model certainty. Dose additivity was assumed in calculating a Combined Exposure-Bioactivity Index (CEBI), the sum of exposure/bioactivity ratios. Combined estrogen bioactivity was also calculated in terms of the percent maximum bioactivity of chemical mixtures in human plasma using a concentration-addition model. Estimated CEBIs vary greatly depending on assumptions used for exposure and bioactivity. In general, CEBI values were <1 when using median of the estimated general population chemical intake rates, while CEBI were ≥1 when using the upper 95th confidence bound for those same intake rates for all chemicals. Concentration-addition model predictions of mixture bioactivity yield comparable results. Based on current in vitro bioactivity data, HTTK methods, and exposure models, combined exposure scenarios sufficient to influence estrogen bioactivity in the general population cannot be ruled out. Future improvements in screening methods and computational models could reduce uncertainty and better inform the potential combined effects of estrogenic chemicals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗
DOE OSTI2021

Tumor Oxygenation by Myo-Inositol Trispyrophosphate Enhances Radiation Response

Tumor hypoxia is a major limiting factor for successful radiation therapy outcomes, with hypoxic cells being up to 3-fold more radiation resistant than normoxic cells; tumor hypoxia creates a tumor microenvironment that is hostile to immune response. Thus, pharmaceutical-induced tumor oxygenation before radiation therapy represents an interesting method to enhance the efficacy of radiation therapy. Myo-inositol trispyrophosphate (ITPP) triggers a decrease in the affinity of oxygen to hemoglobin, which leads to an increased release of oxygen upon tissue demand, including in hypoxic tumors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Senolysis induced by 25-hydroxycholesterol targets CRYAB in multiple cell types

Cellular senescence is a driver of many age-related pathologies. There is an active search for pharmaceuticals termed senolytics that can mitigate or remove senescent cells in vivo by targeting genes that promote the survival of senescent cells. We utilized single-cell RNA sequencing to identify CRYAB as a robust senescence-induced gene and potential target for senolysis. Using chemical inhibitor screening for CRYAB disruption, we identified 25-hydroxycholesterol (25HC), an endogenous metabolite of cholesterol biosynthesis, as a potent senolytic. We then validated 25HC as a senolytic in mouse and human cells in culture and in vivo in mouse skeletal muscle. Thus, 25HC represents a potential class of senolytics, which may be useful in combating diseases or physiologies in which cellular senescence is a key driver.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Influence of alveolar mixing and multiple breaths of aerosol intake on particle deposition in the human lungs

Predictive dosimetry models play an important role in assessing health effect of inhaled particulate matter and in optimizing delivery of inhaled pharmaceutical aerosols. Here, in this study, the commonly used 1D Multiple-Path Particle Dosimetry model (MPPD) was improved by including a mechanistically based model component for alveolar mixing of particles and by extending the model capabilities to account for multiple breaths of aerosol intake. These modifications increased the retained fraction of particles and consequently particle deposition predictions in the deep lung during tidal breathing. Comparison with an existing dataset (J. Aerosol Sci., 99:27-39, 2016) obtained under two breathing conditions referred to as slow and fast breathing showed significant differences in 1 µm particle deposition between predictions (slow: 30±1%, fast: 21±1%) and measurements (slow: 43±9%, fast: 30±5%) when the prior version of MPPD (single breath and no mixing, J. Aerosol Sci., 151:105647, 2021) was used. Adding a mixing model and multiple breaths moved the predictions (slow: 34±2%, fast:25±2%) closer to the range of deposition measurements. For 2.9 µm particles, predictions from both the original (slow: 70±2%, fast: 57±2%) and the revised MPPD model (slow: 71±2%, fast: 59±3%) compared well with experiments (slow: 67±8%, fast: 58±10%). This was expected as suspended fraction of 2.9 µm particles was small and thus the addition of alveolar mixing and multi breath capability only slightly increased the retained fraction for particles of this size and greater. The revised 1D model improves dose predictions in the deep lung and support human risk assessment from exposure to airborne particles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

UV/chlorine advanced oxidation processes for degradation of contaminants in water and wastewater: A comprehensive review

Advanced oxidation processes (AOPs) have recently gained increasing attention as promising alternatives to conventional water and wastewater treatment technologies. As an AOP, there has been a growing interest in UV/chlorine due to its multiple disinfection modes and ability to maintain chlorine residuals in drinking water distribution systems. Compared to the popular UV/H 2 O 2 AOP treatment process, UV/chlorine AOP possesses several advantages including a higher radical production efficiency and a lower power demand. The process has recently been applied for the degradation of a wide range of natural and synthetic organic as well as inorganic compounds including pharmaceutical and personal care products, natural organic matter, herbicides and insecticides, taste and odor, and different types of water and wastewaters. The focus of the present work is to review and evaluate the efficiency of the UV/chlorine process on the oxidation of these pollutants. Here a critical overview is provided of the fundamentals of UV/chlorine, controlling factors of the process performance, and a comparison with other UV-based AOPs. The formation of disinfection by-products (DBPs), one of the concerns associated with UV/chlorine treatment, is discussed and strategies to tackle this obstacle are provided for future studies and applications.

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

Ionic liquid mediated one-step perpendicular assembly in block copolymer thin films for ultrafiltration membrane applications

The rapid rise of oil and gas, petrochemical, food processing, and pharmaceutical industries has added a complex mixture of contaminants like oil, particulates, metals, and other organic compounds to wastewater. Multipurpose membranes with precise pore structure can offer a solution to this problem and block copolymers (BCP) can be used to achieve such structures. Achieving vertical assembly of BCP domains provides a perfect template for developing uniform through film channels for separation and transport of material, besides being useful for the rectification of defects in photolithography patterns. Producing such morphologies may require extensive film/substrate processing and is not always feasible for scale-up. With this article, we demonstrate a facile solution casting method to induce vertical domain assembly in as-cast diblock copolymer thin films in the presence of an ionic liquid (IL) additive that preferentially segregates to one block and neutralizes interfacial interactions. We also show the tunability of domain sizes by controlling the additive concentration. These vertically aligned morphologies are important for the development of next-generation lithographic techniques and form excellent templates for ultrafiltration membranes with uniform pore sizes. In conclusion, this article demonstrates how IL additives can be used to obtain stable non-equilibrium morphologies in as-cast BCP films and the effect of BCP molecular mass, block volume fractions, and IL content on self-assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

In vivo evolution of lactic acid hyper-tolerant Clostridium thermocellum

Lactic acid (LA) has several applications in the food, cosmetics and pharmaceutical industries, as well as in the production of biodegradable plastic polymers, namely polylactides. Industrial production of LA is essentially based on microbial fermentation. Recent reports have shown the potential of the cellulolytic bacterium Clostridium thermocellum for direct LA production from inexpensive lignocellulosic biomass. However, C. thermocellum is highly sensitive to acids and does not grow at pH < 6.0. Improvement of LA tolerance of this microorganism is pivotal for its application in cost-efficient production of LA. In the present study, the LA tolerance of C. thermocellum strains LL345 (wild-type fermentation profile) and LL1111 (high LA yield) was increased by adaptive laboratory evolution. At large inoculum size (10 %), the maximum tolerated LA concentration of strain LL1111 was more than doubled, from 15 g/L to 35 g/L, while subcultures evolved from LL345 showed 50–85 % faster growth in medium containing 45 g/L LA. Gene mutations (pyruvate phosphate dikinase, histidine protein kinase/phosphorylase) possibly affecting carbohydrate and/or phosphate metabolism have been detected in most LA-adapted populations. Although improvement of LA tolerance may sometimes also enable higher LA production in microorganisms, C. thermocellum LA-adapted cultures showed a yield of LA, and generally of other organic acids, similar to or lower than parental strains. Based on its improved LA tolerance and LA titer similar to its parent strain (LL1111), mixed adapted culture LL1630 showed the highest performing phenotype and could serve as a framework for improving LA production by further metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Trithiol ligand provides tumor-targeting 191 Pt-complexes with high molar activity and promising in vivo properties

The Auger electron-emitting radionuclide 191 Pt is a promising candidate for radiopharmaceutical therapy. Herein, we explored novel labeling methods for 191 Pt using thiol-containing ligands to improve the in vivo stability and targeting ability of 191 Pt-labeled complexes. We synthesized dithiol-containing N 2 S 2 and NS 2 ligands, and a trithiol ligand, and then compared their radiochemical reactivity with 191 Pt. [ 191 Pt]Pt-trithiol was synthesized and its biodistribution was evaluated in mice and compared with free 191 Pt. Finally, a 191 Pt-trithiol complex targeting prostate-specific membrane antigen (PSMA): [ 191 Pt]Pt-trithiol-PSMA was developed and evaluated in mice bearing tumor xenografts and compared with a 191 Pt-complex labeled via monothiol-containing Cys ([ 191 Pt]Pt-Cys-PSMA). A comparison of N 2 S 2 , NS 2 , and trithiol showed that the trithiol ligand is the best for producing 191 Pt-labeled compounds in high yield and as a single peak in preparative HPLC. Notably, the trithiol ligand made 191 Pt-labeled compounds and precursors separatable, achieving 191 Pt-labeled products with a high molar activity: 200–400 mCi/μmol (7.4–14.8 GBq/μmol) at EOS. Additionally, [ 191 Pt]Pt-trithiol and [ 191 Pt]Pt-trithiol-PSMA were stable in vivo with rapid clearance compared with free 191 Pt and [ 191 Pt]Pt-Cys-PSMA. [ 191 Pt]Pt-trithiol-PSMA resulted in a low uptake in most normal organs and a high uptake in the kidneys and prostate cancer with PSMA expression. Furthermore, this study demonstrated that a labeling method with trithiol for Pt radionuclides achieves 191 Pt-labeled products with high molar activity. 191 Pt-trithiol-PSMA showed promising in vivo stability and tumor-targeting specificity, which should facilitate the pharmaceutical development of Pt radionuclides for radiopharmaceutical therapy, especially Auger electron cancer therapy.

Auger emitters↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗