Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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

Measuring Electron Diffusion and Constraining the Neutral Current $\pi^0$ Background for Single-Photon Events in MicroBooNE

Liquid Argon Time Projection Chambers (LArTPCs) are a rising technology in the field of experimental neutrino physics. LArTPCs use ionization electrons and scintillation light to reconstruct neutrino interactions with exceptional calorimetric and position resolution capabilities. Here, I present two analyses conducted in the MicroBooNE LArTPC at Fermilab: a measurement of the longitudinal electron diffusion coefficient, $D_L$, in the MicroBooNE detector and a constraint of the systematic uncertainty on MicroBooNE's single-photon analysis due to the dominant neutral current (NC) \piz{} background. Longitudinal electron diffusion modifies the spatial and timing resolution of the detector, and measuring it will help correct for these effects. Furthermore, current measurements of $D_L$ in liquid argon are sparse and in tension with one another, making the MicroBooNE measurement especially valuable. We report a measurement of $3.74^{+0.28}_{-0.29}$ cm$^2$/s. MicroBooNE is searchin g for si ngle-photon events as a potential explanation for the MiniBooNE low-energy excess (LEE) of electron neutrino-like events, which has been interpreted as evidence for low-mass sterile neutrinos. However, this search is overwhelmed by a large NC \piz{} background. By performing a sideband selection of NC \piz{} events, we apply a data-driven rate constraint to the single-photon analysis to reduce the systematic uncertainties. At present, this constraint improves the single-photon analysis' median sensitivity to the LEE-like signal from 0.9$\sigma$ to $1.5\sigma$. This sensitivity is expected to improve significantly as more data become available. Both of these measurements will not only benefit MicroBooNE, but also inform future LArTPC experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Developing novel electrodes with ultralow catalyst loading for high-efficiency hydrogen production in proton exchange membrane electrolyzer cells

Hydrogen plays more crucial roles for decarbonizing the planets and meeting the climate challenges because of its high energy density and zero-emission. It can be produced with proton exchange membrane electrolyzer cells (PEMECs) driven by sustainable and renewable energy resources. Although PEMECs have a number of advantages, including high purity production, quick response, and the ability to operate at high pressure facilitating the gas delivering, their performance and cost greatly hinder their commercial-scale applications. To achieve high-efficiency and cost-reduced hydrogen production in PEMECs, we proposed thin engineered liquid/gas diffusion layers (LGDLs) and associated electrodes, i.e., catalyst-coated LGDLs (CCLGDLs), over conventional porous transport layers (PTLs) and catalyst-coated membranes (CCMs). The research approaches in this project are based on material synthesis, in-situ and ex-situ characterizations, component design and treatment, numerical modeling, and cost analysis. The thin and tunable LGDLs (TT-LGDLs) and CCLGDLs were successfully developed with great performance improvement as demonstrated in lab-scale, bench-scale, and system-scale electrolyzer tests. The electrode thickness was reduced from 370 µm to less than 100 µm with simplified fabrication processes. With the catalytically enhanced Ir-based catalyst coating, the as-developed CCLGDLs with a catalyst loading of 0.34 mg Ir /cm 2 achieved a cell performance of 1.77 V at 2 A cm -2 , exhibiting the catalyst mass activity enhanced by >20 times with significant catalyst saving over conventional catalyst cell design. In-situ PEMEC characterizations, including the current distribution mapping and high-speed and multiscale visualizations, were conducted for a deeper understanding of mass transport and electrochemical reactions within an electrolyzer with LGDLs and CCLGDLs. A 2D cell model was developed and validated for the enhanced performance on TT-LGDL through reducing ohmic losses due to nonuniform hydration and water transport. Further, the cost analysis results have shown a path to move beyond equivalency and surpass costs associated with the project baseline. In this project, the design and fabrication of TT-LGDLs and CCLGDLs will contribute to the performance enhancement, manufacturing simplification, and cost reduction for PEMECs and other energy conversion devices, thus shortening their pathways towards commercialization. This project also provides a good foundation for furthering the in-situ reaction interface research.

08 HYDROGEN↗
DOE OSTI2021

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Preparation of pH ‐responsive cellulose nanofibril/sodium alginate based hydrogels for drug release

Abstract In this work, hydrogels composed of cellulose nanofibril (CNF) and sodium alginate (SA) were prepared as a pH‐responsive release system by controlling the amount of CNF and SA. Scanning electron microscopy, thermogravimetric analysis, and swelling rate methods were used to characterize. The results showed that when the mass ratio of CNF to SA was 2:1, CNF/SA hydrogels with better morphology were obtained. Due to the deprotonation or protonation of calcium alginate, as the pH changing from 1.5 to 11.0, the swelling ratio increased from 7 times to 19 times compared with its dry weight. The drug release behaviors of the CNF/SA hydrogels using ibuprofen (IBU) as a model drug were investigated in vitro. The release of IBU was faster in alkaline conditions and slower in acidic conditions. The drug release mechanism under pH stimulation could be interpreted by the Fickian diffusion and erosion effect based on Korsmeyer‐Peppas model.

Li, Yuhang↗
DOE OSTI2022

High-performance all-solid-state Li–S batteries enabled by an all-electrochem-active prelithiated Si anode

The development of all-solid-state Li–S batteries has been greatly impeded by dendrite growth and dendrite penetration, which are both related to the Li metal anode. As a compromised alternative, we report lithium silicide (Li x Si) as a dendrite-free and high-capacity anode with Li source. Here, the as-synthesized Li x Si is relatively soft, highly electronically conductive, and with a high Li diffusivity. These distinctive properties make Li x Si anode viable as an “all-electrochem-active” electrode (consisted of 100 wt.% Li x Si). Compared with the typical composite electrode, the all-electrochem-active electrode not merely maximizes the electrode-level energy density but also minimizes the electrolyte-related interfacial degradation. Li x Si symmetric cell demonstrates a reversible cycling at 4 mA cm –2 for over 320 h. Stress change and morphological evolution of the Li x Si electrode are investigated upon dealloying/alloying. When paired with a S cathode (active mass loading of 3 mg cm –2 ), Li x Si–S full cell shows a good cycling behavior over 500 cycles and rate performance (69% capacity retained at 1.2C) even at 25 °C.

25 ENERGY STORAGE↗
DOE OSTI2022

Catalyst layer formulations for slot-die coating of PEM fuel cell electrodes

The series of electrodes were fabricated by the scalable and manufacturable slot-die coating method for proton exchange membrane fuel cell (PEMFC) application. The inks with different amounts of solids were studied by rheological methods in order to establish a coating window with minimum manufacturing defects. The obtained electrodes were characterized by SEM, AFM, and optical microscopy, which showed that they were uniform and homogeneous with minimum defects. The electrochemical evaluation of the manufactured gas diffusion electrodes (GDE) showed that the main characteristics of the electrodes, like electrochemical surface area, proton resistivity, and double layer capacitance, were found to be close for all samples confirming the reproducibility of the slot-die process. Additionally, we studied the effects of membrane thickness on the performance of the GDE membrane electrode assemblies and determined that a decrease in membrane thickness favored the performance. Furthermore, the obtained results clearly demonstrated the applicability and feasibility of the approach for the Manufacturing of catalyst layers for the fuel cell application with potential for future mass production.

08 HYDROGEN↗
DOE OSTI2022

Local stress within a granular molecular solvent matrix, a mechanism for individual ion hydration

The activity of water is important in many natural systems and directed processes. To better understand water activity, the hydration of individual ions is approximated for a recently developed speciation-based solution model. Here a matrix of values for hydration of complete salts (anion and cation) was used as source data to calculate individual ion hydration values. Some electrolytes were excluded from calculations of individual ion hydration due to the prevalence of ion pairing and resulting liberation of water upon ion pairing. Evaluation of the periodic trends emerging from individual ion hydration values suggests that ionic volume, relative to the volume of molecular water, dictates both modeled hydration requirements (related to vapor-liquid-equilibrium-derived activity) as well as properties historically described by hydrodynamic radius (ion diffusion and solution viscosity). When water is conceptualized as a granular three-dimensional matrix, solutes may function as substitutional defects (dilation or constriction) inducing local compressive or tensile stress. These results have implications on the mechanism of the Hofmeister series and series reversals.

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

Multiwavelength Observations of 2HWC J1928+177: Dark Accelerator or New TeV Gamma-Ray Binary?

2HWC J1928+177 is a Galactic TeV gamma-ray source detected by the High Altitude Water Cherenkov (HAWC) Observatory up to ~56 TeV. The HAWC source, later confirmed by the High Energy Stereoscopic System, still remains unidentified as a dark accelerator since there is no apparent supernova remnant or pulsar wind nebula detected in the lower-energy bands. The radio pulsar PSR J1928+1746, coinciding with the HAWC source position, has no X-ray counterpart. In this work, our SED modeling shows that inverse Compton scattering in the putative pulsar wind nebula can account for the TeV emission only if the unseen nebula is extended beyond $r\sim 4^{\prime} $. Alternatively, TeV gamma-rays may be produced by hadronic interactions between relativistic protons from an undetected supernova remnant associated with the radio pulsar and a nearby molecular cloud G52.9+0.1. NuSTAR and Chandra observations detected a variable X-ray point source within the HAWC error circle, potentially associated with a bright infrared (IR) source. The X-ray spectra can be fitted with an absorbed power-law model with N H = (9 ± 3) × 10 22 cm –2 and Γ X = 1.6 ± 0.3 and exhibit long-term X-ray flux variability over the last decade. If the X-ray source, possibly associated with the IR source (likely an O star), is the counterpart of the HAWC source, it may be a new TeV gamma-ray binary powered by collisions between the pulsar wind and stellar wind. Follow-up X-ray observations are warranted to search for diffuse X-ray emission and determine the nature of the HAWC source.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Improving Thermal Stability of Perovskite Solar Cells by Suppressing Ion Migration Using Copolymer Grain Encapsulation

Thermal stability of organic–inorganic hybrid perovskites (OIHPs) remains as one of the critical challenges against the stable operation of perovskite solar cells (PSCs) in direct sunlight with elevated temperatures. Here, we show that the addition of a polystyrene-co-polyacrylonitrile (SAN) copolymer can significantly enhance thermal stability of OIHPs and improve the stability of the corresponding PSCs by suppressing the migration of organic cations in OHIP. The methylammonium lead iodide (MAPI) with SAN incorporated within the perovskite layer featured a superior thermal stability compared to pure MAPI without SAN, only displaying an average of 5–15% decrease in PCE even after continuous thermal aging for 24 h at 100 °C. The secondary ion mass spectrometry revealed that the thermal degradation of the pure MAPI was largely associated with MA+ out-migration. Conducting atomic force microscopy analysis further indicated that the incorporated SAN led to a suppression of ionic currents present at the grain boundaries of the perovskite film, which was understood by high immiscibility between SAN and MA+ components as confirmed by the experimentally estimated Flory–Huggins parameter between them. This study newly identifies a potential of using polymer grain encapsulation for enhancing thermal stability of OIHPs and their solar cell performance by suppressing the out-diffusion of cationic organic components.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Editors’ Choice—Flooded by Success: On the Role of Electrode Wettability in CO 2 Electrolyzers that Generate Liquid Products

Economic operation of carbon dioxide (CO 2 ) electrolyzers generating liquid products will likely require high reactant conversions and product concentrations, conditions anticipated to challenge existing gas diffusion electrodes (GDEs). Notably, electrode wettability will increase as lower surface tension products (e.g., formic acid, alcohols) are introduced into electrolyte streams, potentially leading to flooding. To understand the hydraulically stable operating envelopes in mixed aqueous-organic liquid domains, we connect intrinsic electrode wettability descriptors to operating parameters such as electrolyte flow rate and current. We first measure contact angles of water-organic product dilutions on polytetrafluoroethylene (PTFE) and graphite surfaces as planar analogues for GDE components. We then use material balances around the reactive gas-liquid-solid interface to calculate product mass fractions as functions of water sweep rate and current. Product composition maps visualize the extent to which changes in cell performance influence capillary pressure, a determinant of GDE saturation. Analyses suggest that formic acid mixtures pose little risk for GDE flooding across a wide range of conditions, but effluents containing <30% alcohol by mass may cause flooding. This study reveals opportunities to integrate microstructural features and oleophobic surface treatments into GDEs to repel aqueous-organic mixtures and expand the window of stable operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2022

Challenges to Reliable Production Nitrogen Doping of Nb for SRF Accelerating Cavities

Over the last several years, alloying of the surface layer of niobium SRF cavities has been demonstrated to beneficially lower the superconducting RF surface resistance. Nitrogen, titanium, and oxygen have all been demonstrated as effective alloying agents, occupying interstitial sites in the niobium lattice within the RF penetration depth and even deeper, when allowed to thermally diffuse into the surface at appropriate temperatures. The use of nitrogen for this function has been often termed ’nitrogen doping’ and is being applied in the LCLS-II and LCLS-II HE projects. We report characterization studies of the distribution of nitrogen into the exposed niobium surface and how such distribution is affected by the quality of the vacuum furnace environment in which the doping takes place, and the complexity of nitride crystal growth on different grain orientations of surface niobium. Using state-of-the-art quantification methods by dynamic secondary ion mass spectrometry (SIMS) depth profiling in niobium, we identify several phenomena involving furnace-sourced contamination. We also highlight a potential issue with N₂ flow constraints from the flange ’caps’ used during heat treatments.

Reece, Charles↗
DOE OSTI2023

Nonequilibrium Behavior in Isoprene Secondary Organic Aerosol

Recent studies have shown that instantaneous gas-particle equilibrium partitioning assumptions fail to predict SOA formation, even at high relative humidity (~85%), and photochemical aging seems to be one driving factor. In this study, we probe the minimum aging time scale required to observe nonequilibrium partitioning of semivolatile organic compounds (SVOCs) between the gas and aerosol phase at ~50% RH. Seed isoprene SOA is generated by photo-oxidation in the presence of effloresced ammonium sulfate seeds at <1 ppbv NO x , aged photochemically or in the dark for 0.3–6 h, and subsequently exposed to fresh isoprene SVOCs. Our results show that the equilibrium partitioning assumption is accurate for fresh isoprene SOA but breaks down after isoprene SOA has been aged for as short as 20 min even in the dark. Modeling results show that a semisolid SOA phase state is necessary to reproduce the observed particle size distribution evolution. The observed nonequilibrium partitioning behavior and inferred semisolid phase state are corroborated by offline mass spectrometric analysis on the bulk aerosol particles showing the formation of organosulfates and oligomers. Here, the unexpected short time scale for the phase transition within isoprene SOA has important implications for the growth of atmospheric ultrafine particles to climate-relevant sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sulphate resistance of low‐clinker engineered cementitious composites examined by MicroXRF imaging

Abstract Engineered cementitious composites (ECC) are a class of high‐performing fibre‐reinforced cementitious materials recognised for their increased ductility and durability compared to conventional cement‐based materials, owing to their autogenously controlled tight crack widths, even when subjected to high strains. To reduce ECC's environmental impact, this research examines the use of a low‐clinker binder − limestone‐calcined clay cement (LC3) − as an alternative to portland cement (PC), along with fly ash to further reduce the clinker proportion and the embodied CO 2 of the formulations. In conventional concrete, LC3 hydrates to a denser microstructure resulting from the synergistic reaction between limestone and calcined clay. At the lower water contents typical of ECC and with the presence of fly ash, the influence of the binder composition on the microstructure is difficult to anticipate. To examine the influence of these compositional variables on microstructure, permeability and durability, the sulphate resistance of LC3‐based ECC is explored. Specifically, the ECC‐LC3 blends are designed with high clinker replacement rate of 75% by mass of binder and contain either conventional fly ash or reclaimed fly ash at 50% by mass of binder. Expansion of ECC‐LC3 samples subjected to standard sodium sulphate test conditions was measured up to 12 months and the depth of penetration of sulphates into the ECC‐LC3 of varying compositions was quantified using micro‐X‐Ray Fluorescence (microXRF) imaging and modelling. The expansion results show that the ECC‐LC3 formulations performed better than the PC samples and can provide adequate resistance to external sulphate attack, even when reclaimed fly ashes are used in place of the conventional ash. In addition, the shallow penetration of sulphate into these cementitious composites demonstrates the low diffusion coefficients values that were determined using the quantitative data from MicroXRF imaging.

Microscopy↗
DOE OSTI2023

Perivascular network segmentations derived from high-field MRI and their implications for perivascular and parenchymal mass transport in the rat brain

A custom segmentation workflow was applied to ex vivo high-field MR images of rat brains acquired following in vivo intraventricular contrast agent infusion to generate maps of the perivascular spaces (PVS). The resulting perivascular network segmentations enabled analysis of perivascular connections to the ventricles, parenchymal solute clearance, and dispersive solute transport within PVS. Numerous perivascular connections between the brain surface and the ventricles suggest the ventricles integrate into a PVS-mediated clearance system and raise the possibility of cerebrospinal fluid (CSF) return from the subarachnoid space to the ventricles via PVS. Assuming rapid solute exchange between the PVS and CSF spaces primarily by advection, the extensive perivascular network decreased the mean clearance distance from parenchyma to the nearest CSF compartment resulting in an over 21-fold reduction in the estimated diffusive clearance time scale, irrespective of solute diffusivity. This corresponds to an estimated diffusive clearance time scale under 10 min for amyloid-beta which suggests that the widespread distribution of PVS may render diffusion an effective parenchymal clearance mechanism. Additional analysis of oscillatory solute dispersion within PVS indicates that advection rather than dispersion is likely the primary transport mechanism for dissolved compounds greater than 66 kDa in the long (> 2 mm) perivascular segments identified here, although dispersion may be significant for smaller compounds in shorter perivascular segments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗
DOE OSTI2024

Complete Development of Critical Capabilities for TRISO Fission Product Source Term Calculations and Quantify Mechanisms for Pd Penetration of SiC

Overall fission product (FP) release will be an important consideration for the licensing and deployment of advanced reactors utilizing tristructural isotropic (TRISO) fuels. This work focuses on enhancing and applying the BISON models needed to predict FP transport within TRISO particles and particle failure probability, both of which factor directly into release predictions. Specifically, this report details (1) the development of the models needed to predict palladium (Pd) conservation at the engineering scale and the application of those models to characterize Pd fluxes for input into a mechanistic multiscale model for Pd penetration; (2) the refinement of sorption mass transfer models and the development of models for trapping in porous layers, which were applied and compared to particle scans from AGR-2 to provide proof of concept for a method of particle-scale validation that may reduce uncertainties compared to compact-scale validation using data from integral effects tests; (3) the development of a failure-statistics-informed, mesh-independent methodology for applying smeared cracking, enabling further study of the localized multiphysics behaviors associated with cascading particle failure mechanisms; and (4) the preliminary characterization of those coupled multiphysics particle failure behaviors using smeared, nonretentive diffusivities to provide a baseline for future study and to guide ongoing engineering applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗