Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mobile computing”

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 1,405 records · Page 78

Compact mobile-reader system for two-way wireless communication, tracking and status monitoring for transport safety and security

A system for monitoring a plurality radio frequency identification tags is described. The system uses at least one set of radio frequency identification tags. Each tag is attached to a container and includes several sensors for detecting physical conditions of said container. The system includes at least one autonomous intermediate reader in wireless communication with the frequency identification tags. The intermediate reader includes external wireless communication system, intermediate reader logic controller, and a self-contained rechargeable power supply. The system uses a central status reporting system in communication the intermediate reader.

97 MATHEMATICS AND COMPUTING↗

Identification and recovery of rare-earth permanent magnets from waste electrical and electronic equipment

Highlights: • Different WEEE were analyzed and the rare-earth permanent magnets were collected. • Chemical composition evolution of Nd-Fe-B alloys over the years is illustrated. • Loudspeakers from laptops and loudspeakers from flat TV contained Nd-Fe-B magnets. • The amount of Nd-Fe-B could double if the loudspeakers streams are included. - Abstract: Nd-Fe-B permanent magnets are a strategic material for a number of emerging technologies. They are a key component in the most energy efficient electric motors and generators, thus, they are vital for energy technologies, industrial applications and automation, and future forms of mobility. Rare earth elements (REEs) such as neodymium, dysprosium and praseodymium are also found in waste electrical and electronic equipment (WEEE) in volumes that grow with the technological evolution, and are marked as critical elements by the European Commission due to their high economic importance combined with significant supply risks. Recycling could be a good approach to compensate for the lack of rare earths (REs) on the market. However, less than 1% of REs are currently being recycled, mainly because of non-existing collection logistics, lack of information about the quantity of RE materials available for recycling and recycling-unfriendly product designs. To improve these lack of information, different waste streams of electrical and electronic equipment from an industrial recycling plant were analyzed in order to localize, identify and collect RE permanent magnets of the Nd-Fe-B type. This particular type of magnets were mainly found in hard disk drives (HDDs) from laptops and desktop computers, as well as in loudspeakers from compact products such as flat screen TVs, PC screens, and laptops. Since HDDs have been investigated thoroughly by many authors, this study focusses on other potential Nd-Fe-B resources in electronic waste. The study includes a systematic survey of the chemical composition of the Nd-Fe-B magnets found in the selected waste streams, which illustrates the evolution of the Nd-Fe-B alloys over the years. The study also provides an overview over the types of magnets integrated in different waste electric and electronic equipment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Rainbow Coverage Path Planning for a Patrolling Mobile Robot With Circular Sensing Range

In this paper, we proposed a coverage path (CP) planning approach for a mobile robot moving in a dynamically changing environment. Our family of algorithms begins by finding the first path around the center of interior obstacles/disks of a given target region (TR) through the Graham scan algorithm. The first path is then used as a foundational path to generate a collision free, first order differentiable, and observable CP through a complete set of input and output transformation algorithms which together form the rainbow CP planning approach. The last algorithm of the rainbow CP planning approach finds a sufficient number of observation points on the CP needed to observe the TR. The novelty of our rainbow CP planning approach is that it partitions the TR into different shapes with different properties needed to obtain complete coverage while achieving first order path differentiability. The main technical contributions of the proposed approach is to provide a holistic solution that segments any TR, uses triangulation to determine the line of sights and observation points, and computes the collision-free CP within a quadratic runtime. The proposed method can be readily generalized to address problems of higher dimensions, and it is scalable with respect to the size of TR and the number of robots. Computer simulations are used to illustrate the effectiveness and correctness of the proposed approach.

54 ENVIRONMENTAL SCIENCES↗

Audits Made Simple

A company just got notified there is a big external audit coming in 3 months. Getting ready for an audit can be challenging, scary, and full of surprises. This Gold Paper describes a typical audit from notification of the intent to audit through disposition of the final report including Best Practices, Opportunities for Improvement (OFI), and issues that must be fixed. Good preparation can improve the chances of success. Ensuring the auditors understand the environment and requirements is paramount to success. It helps the auditors understand that the enterprise really does think that security is important. Understanding and following a structured process ensures a smooth audit process. Ensuring follow-up on OFIs and issues in a structured fashion will also make the next audit easier. It is important to keep in mind that the auditors will use the previous report as a starting point. Now the only worry is the actual audit and subsequent report and how well the company has done. “Virtualization, mobilization, and cloud technology have created new points of entry into business, leaving them vulnerable to covert cyber-attacks,” states a report by Ernst & Young. Executives are struggling with the potential for a data breach, and resolving this business risk is high on their agenda. While most audit committee members are financially solid, they lack knowledge of technology issues, making understanding of this risk difficult. Success with an audit helps management understand the risk management is accepting, highlight gaps in security, and generate ideas on how to improve. These results are an important aspect of managing the business risk.

97 MATHEMATICS AND COMPUTING↗

POINT: Partially Observable Imitation Network for Traffic Signal Control

Smart traffic signals bring together transportation infrastructure and advance technologies to improve the mobility and efficiency of urban transportation network. Adaptive traffic signal control studies can be categorized into modeling-based approaches and learning-based approaches. In order to take advantages of these two systems, this study developed an offline-online combined Partial Observable Imitation Network for Traffic signal control (POINT). In the offline system, the traffic signal timing optimization problem was formulated as a Mixed Integer Nonlinear Programming (MINLP) given complete traffic information, i.e., second-by-second speeds and locations of all vehicles. Furthermore, the objective of MINLP is to minimize total travel delays considering individual vehicle trajectories under Connected Vehicle (CV) environment. The calculated optimal solutions under various traffic conditions were considered as the ”expert” decisions. In the online system, an imitation neural network model was developed to learn the ”expert” signal plans generated from offline system. Given partial observable traffic conditions in real time, e.g., the aggregate-level of traffic volume, the POINT model can compute the signal timing parameters in the online system. The numerical results demonstrated that the proposed method outperformed other state-of-the-art signal control method under high and unbalanced traffic demand levels in terms of reducing travel delays and queue length.

33 ADVANCED PROPULSION SYSTEMS↗

Estimating the vibration level of an L-shaped beam using power flow techniques

The response of one component of an L-shaped beam, with point force excitation on the other component, is estimated using the power flow method. The transmitted power from the source component to the receiver component is expressed in terms of the transfer and input mobilities at the excitation point and the joint. The response is estimated both in narrow frequency bands, using the exact geometry of the beams, and as a frequency averaged response using infinite beam models. The results using this power flow technique are compared to the results obtained using finite element analysis (FEA) of the L-shaped beam for the low frequency response and to results obtained using statistical energy analysis (SEA) for the high frequencies. The agreement between the FEA results and the power flow method results at low frequencies is very good. SEA results are in terms of frequency averaged levels and these are in perfect agreement with the results obtained using the infinite beam models in the power flow method. The narrow frequency band results from the power flow method also converge to the SEA results at high frequencies. The advantage of the power flow method is that detail of the response can be retained while reducing computation time, which will allow the narrow frequency band analysis of the response to be extended to higher frequencies.

Cuschieri, J. M.↗

Evaluation of imaging quality for flat-panel detector based low dose C-arm CT system

The image quality associated with the extent of the angle of gantry rotation, the number of projection views, and the dose of X-ray radiation was investigated in flat-panel detector (FPD) based C-arm cone-beam computed tomography (CBCT) system for medical applications. A prototype CBCT system for the projection acquisition used the X-ray tube (A-132, Varian inc.) having rhenium-tungsten molybdenum target and flat panel a-Si X-ray detector (PaxScan 4030CB, Varian inc.) having a 397 x 298 mm active area with 388 μm pixel pitch and 1024 x 768 pixels in 2 by 2 binning mode. The performance comparison of X-ray imaging quality was carried out using the Feldkamp, Davis, and Kress (FDK) reconstruction algorithm between different conditions of projection acquisition. In this work, head-and-dental (75 kVp/20 mA) and chest (90 kVp/25 mA) phantoms were used to evaluate the image quality. The 361 (30 fps x 12 s) projection data during 360 deg. gantry rotation with 1 deg. interval for the 3D reconstruction were acquired. Parke weighting function were applied to handle redundant data and improve the reconstructed image quality in a mobile C-arm system with limited rotation angles. The reconstructed 3D images were investigated for comparison of qualitative image quality in terms of scan protocols (projection views, rotation angles and exposure dose). Furthermore, the performance evaluation in image quality will be investigated regarding X-ray dose and limited projection data for a FPD based mobile C-arm CBCT system. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Computational Assessment of Naturally Occurring Neutron and Photon Background Radiation Produced by Extraterrestrial Sources

In this study, a computational assessment of the variation in terrestrial neutron and photon background from extraterrestrial sources is presented. The motivation of this assessment is to evaluate the practicality of developing a tool or database to estimate background in real time (or near–real time) during an experimental measurement or to even predict the background for future measurements. The extraterrestrial source focused on during this assessment is naturally occurring galactic cosmic rays (GCRs). The MCNP6 transport code was used to perform the computational assessment. However, the GCR source available in MCNP6 was not used. Rather, models developed and maintained by NASA were used to generate the GCR sources. The largest variation in both neutron and photon background spectra was found to be caused by changes in elevation on Earth's surface, which can be as large as an order of magnitude. All other perturbations produced background variations on the order of a factor of 3 or less. The most interesting finding was that ~80% and 50% of terrestrial background neutrons and photons, respectively, are generated by interactions in Earth's surface and other naturally occurring and man-made objects near a detector of particles from extraterrestrial sources and their progeny created in Earth's atmosphere. In conclusion, this assessment shows that it will be difficult to estimate the terrestrial background from extraterrestrial sources without a good understanding of a detector's surroundings. Therefore, estimating or predicting background during a measurement environment like a mobile random search will be difficult.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Low Pressure Gas Effects On The Potency Of An Electron Beam On Ceramic Fabric Materials For Space Welding

This investigation was undertaken to evaluate if molten metal or electron beam impingement could damage or burn through the fabric of the astronauts Extravehicular Mobility Unit (EMU) during electron beam welding exercises performed in space. An 8 kilovolt electron beam with a current in the neighborhood of 100 milliamps from the Ukrainian space welding "Universal Hand Tool" burned holes in Nextel AF-62 ceramic cloth designed to withstand temperatures up to 1427 C. The burnthrough time was on the order of 8 seconds at standoff distances between UHT and cloth ranging from 6 to 24 inches. At both closer (2") and farther (48") standoff distances the potency of the beam against the cloth declined and the burnthrough time went up significantly. Prior to the test it had been expected that the beam would lay down a static charge on the cloth and be deflected without damaging the cloth. The burnthrough is thought to be an effect of partial transmission of beam power by a stream of positive ions generated by the high voltage electron beam from contaminant gas in the "vacuum" chamber. A rough quantitative theoretical computation appears to substantiate this possibility.

Nunes, Arthur C., Jr.↗

Carbon Monoxide Accumulation in the Extravehicular Mobility Unit

Introduction: Life support technology in large closed systems like submarines and space stations catalyzes carbon monoxide (CO) to carbon dioxide, which is easily removed. However, in a small system like the Extravehicular Mobility Unit (EMU), spacesuit, CO from exogenous (contaminated oxygen (O (sub 2) supply) and endogenous (human metabolism) sources will accumulate in the free suit volume. The free volume becomes a sink for CO that is rebreathed by the astronaut. The accumulation through time depends on many variables: the amount absorbed by the astronaut, the amount produced by the astronaut (between 0.28 and 0.34 moles per hour per kilogram)[1], the amount that enters the suit from contaminated O (sub 2), the amount removed through suit leak, the free volume of the suit, and the O (sub 2) partial pressure[2], just to list a few. Contamination of the EMU O (sub 2) supply with no greater than 1 part per million CO was the motivation for empirical measurements from CO pulse oximetry (SpCO) as well as mathematical modeling of the EMU as a rebreather for CO. Methods: We developed a first-order differential mixing equation as well as an iterative method to compute CO accumulation in the EMU. Pre-post measurements of SpCO (Rad-57, Masimo Corporation) from EMU ground training and on-orbit extravehicular activities (EVAs) were collected. Results: Initial modeling without consideration of the astronaut as a sink but only the source of CO showed that after 8 hours breathing 100 percent O (sub 2) with a 10 milliliter per minute (760 millimeters Hg at 21 degrees Centigrade standard) suit leak, an endogenous production rate of 0.23 moles per hour per kilogram for a 70 kilogram person with 42 liters (1.5 cubic feet) free suit volume resulted in a peak CO partial pressure (pCO) of 0.047 millimeters Hg at 4.3 pounds per square inch absolute (222 millimeters Hg). Preliminary results based on a 2008 model[3] with consideration of the astronaut as a sink and source of CO suggests that most of the rebreathed CO stays bound to hemoglobin and myoglobin and; therefore, pCO only increased to 0.002 millimeters Hg in the EMU. Hemoglobin saturation after 8 hours was an insignificant 0.4 percent compared to about 4 percent for cigarette smokers in the general population. This preliminary modeling result supplements 11 pre-post index finger SpCO measurements from EMU ground training (mean 0.5 percent versus 1 percent, probability equal to 0.41 from paired t-test) and 10 on-orbit pre-post EVAs (mean 1.5 percent versus 1.1 percent, probability equal to 0.17 from paired t-test) that showed no consistent increase, at least no increase outside the accuracy of the oximeter (1 percent display resolution with plus or minus 3 percent Standard Deviation). Discussion: Simulations continue, but a preliminary conclusion is that rebreathing endogenous CO accumulating in the EMU is not a serious medical issue during EVA. The absence of CO poisoning signs or symptoms following hundreds of EVAs is also good empirical evidence that corroborates the limited SpCO measurements and preliminary modeling results.

Conkin, J.↗

Privacy Vulnerability of Published Anonymous Mobility Traces

Mobility traces of people and vehicles have been collected and published to assist the design and evaluation of mobile networks, such as large-scale urban sensing networks. Although the published traces are often made anonymous in that the true identities of nodes are replaced by random identifiers, the privacy concern remains. This is because in real life, nodes are open to observations in public spaces, or they may voluntarily or inadvertently disclose partial knowledge of their whereabouts. Thus, snapshots of nodes’ location information can be learned by interested third parties, e.g., directly through chance/engineered meetings between the nodes and their observers, or indirectly through casual conversations or other information sources about people. In this paper, we investigate how an adversary, when equipped with a small amount of the snapshot information termed as side information, can infer an extended view of the whereabouts of a victim node appearing in an anonymous trace. Our results quantify the loss of victim nodes’ privacy as a function of the nodal mobility, the inference strategies of adversaries, and any noise that may appear in the trace or the side information. Generally, our results indicate that the privacy concern is significant in that a relatively small amount of side information is sufficient for the adversary to infer the true identity (either uniquely or with high probability) of a victim in a set of anonymous traces. For instance, an adversary is able to identify the trace of 30%-50% of the victims when she has collected 10 pieces of side information about a victim.

97 MATHEMATICS AND COMPUTING↗

Fission product distribution in irradiated safety-tested and as-irradiated AGR-2 TRISO particles

In this work, two tristructural isotropic (TRISO)-coated nuclear fuel particles were examined by electron probe microanalysis (EPMA) as part of the Advanced Gas Reactor program. The compacts’ average irradiation temperatures ranged from approximately 1260 to 1290 °C. One particle was examined in the as-irradiated condition, while the other was subject to 1600 °C post-irradiation safety testing. This study was undertaken to test a newly-developed EPMA technique to determine fission product masses in TRISO particles on a layer-by-layer basis, and to compare fission product distributions between an as-irradiated and safety-tested particle. Fission product concentration profiles were collected along two radii in each particle, with measured concentrations used to compute the fission product mass in each TRISO particle layer. These measured masses were then compared to those predicted from ORIGEN modeling calculations. Data collected from these measurements show that for these two particles, masses determined via EPMA were within ± 20% of the calculated masses for the rare-earth elements, Mo, Zr, Cs, I, and Pd. Elements that tend to be less homogeneously distributed include Sr, Te, Eu, Ag, and possibly Ba. Measured Ag masses differed by more than 40% from the calculated mass. Lanthanides other than Eu remain primarily within the fuel kernel in the as-irradiated particle but in the safety-tested particle these element masses were divided approximately equally between the kernel and kernel periphery. In both particles, the majority of Sr and Eu accumulated in the carbon-rich kernel periphery, although in the safety-tested particle, Sr and Eu accumulated farther from the fuel kernel than occurred with irradiation alone. A greater mass fraction of mobile elements, such as Cs and I accumulated in the buffer and IPyC in the safety-tested particle as compared to the as-irradiated particle. When fully developed and tested, this mass balance approach to TRISO particle analysis has the potential to provide insight into fuel behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Infrastructure-Based Cooperative Perception at a Traffic Intersection: Overview and Challenges: Preprint

Recent advancement in autonomous driving vehicles and V2X communication has attracted increasing attention towards Intelligent Transportation Systems to build a safe and reliable traffic intersection. However, most of the systems are still at the initial stages and require significant progress to become a reality. This paper presents an overview of NREL Infrastructure Perception and Control (IPC) framework which is an open-source track-data fusion engine which takes input from infrastructure-based perception sensors and cooperatively shared messages from Connected Autonomous Vehicles (CAVs) and Connected Vehicle (CVs) and the challenges associated with deploying such cooperative perception framework at a four-way traffic intersection in the city of Colorado Springs, CO, USA. The sensor data is collected by deploying two radars and two LiDAR sensors on the IPC mobile lab and two radars on diagonally opposite traffic poles at the proposed intersection. The sensor output results imply the need for rapid sensor calibration to bring the collective perception to a common coordinate frame, the importance of time synchronization between the sensors in order to capture accurate spatial and temporal alignment of the objects, and the need for a health monitoring system with fail safe closed-loop detection model for real-time deployment.

camera↗

Renewable Energy and Efficiency Technologies in Scenarios of U.S. Decarbonization in Two Types of Models: Comparison of GCAM Modeling and Sector-Specific Modeling

Energy system projections from analytic models inform actions ranging from short-term and local decisions, such as technology and infrastructure deployment, to global and long-term negotiations and targets. Computational limits require the designers of these models to trade off between coverage and resolution. Some models, such as the Global Change Analysis Model (GCAM), represent all energy sources and uses but at a relatively coarse level of resolution. GCAM balances global supply and demand of all energy carriers by endogenously projecting prices for energy sources and costs of greenhouse gas mitigation while capturing interlinkages between the energy system, water, agriculture and land use, the economy, and the climate. This global model was used to frame the Long-Term Strategy released by the White House in 2021 and has been used to inform national and global economy-wide decarbonization discussions and strategy development for decades. Other models instead focus on a portion of the energy sector with greater detail and resolution. The Regional Energy Deployment System (ReEDS) electricity-sector model, for example, projects capacity expansion with an emphasis on integration of variable renewable energy into the grid of the future. The Transportation Energy and Mobility Pathway Options (TEMPO) transportation-sector model enables analysis of household choices in adoption, charging, and use of electric vehicles. The Scout buildings-sector model supports detailed consideration of the policies and markets that can accelerate the adoption of energy conservation measures in buildings. Such sector-specific models are instrumental in informing technology research, sectoral planning strategies, and sector-specific aspects of greenhouse gas (GHG) mitigation strategies in the United States. These global and sector-specific modeling approaches can complement each other. The global approach ensures consistent, endogenous energy pricing and resource allocation, which can substantially diverge from current conditions in transformative scenarios, while the sector-specific approach facilitates representation of granular details across spatial, temporal, technological, and market dimensions that enable exploration of particular interactions and trade-offs. This report presents the results of recent work to explore the differences and tradeoffs between these approaches by comparing GCAM with the sector-specific ReEDS, TEMPO, and Scout models. The report compares both model structures and results, and discusses their potential relevance and applications.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Ultrasonic imaging in liquid sodium

The fourth generation of nuclear reactor can use liquid sodium as the core coolant. When the reactor is operating, sodium temperatures can reach up to 600 deg. C. During maintenance periods, when the reactor is shut down, the coolant temperature is reduced to 200 deg. C. Because molten sodium is optically opaque, ultrasonic imaging techniques are developed for maintenance activities. Under-sodium imaging aims at i) checking the health of immersed structures. It should also allow ii) to assess component degradation or damage as cracks and shape defects as well as iii) the detection of lost objects. The under-sodium imaging system has to sustain high temperature (up to 300 deg. C) and hostility of the sodium environment. Furthermore, specific constraints such as transducers characteristics or the limited sensor mobility in the reactor vessel have to be considered. This work focuses on developing a methodology for detecting damages such as crack defects with ultrasound devices. Surface-breaking cracks or deep cracks are sought in the weld area, as welds are more subject to defects. Traditional methods enabled us to detect emerging cracks of submillimeter size with sodium-compatible high-temperature transducer. The presented approach relies on making use of prior knowledge about the environment through the implementation of differential imaging and time-reversal techniques. Indeed, this approach allows to detect a change by comparison with a reference measurement and by focusing back to any change in the environment. It is a means of analysis and understanding of the physical phenomena making it possible to design more effective inspection strategies. Difference between the measured signals reveals the acoustic field scattered by a perturbation (a crack for instance), which may occur between periodical measurements. The imaging method relies on the adequate combination of two computed ultrasonic fields, one forward and one adjoint. The adjoint field, which carries the information about the defects, is analogous to a time-reversal operation. One of the interests of the presented method is that the time-reversal operation is not done experimentally but numerically. Numerical simulations have been carried out to validate the practical relevance of this approach. The preliminary numerical results show a nice agreement between the guessed and the actual positions of the defect. After water-tests, in sodium-tests must be done in order to validate the water/sodium transposition. For this purpose, an under-sodium device is under development, which can move the transducers with four degrees of freedom in a 1.5 m{sup 3} sodium pot. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Quantum Technologies for UAS (QTech)

Recent advances in small Unmanned Aerial System (sUAS) technologies lower the barriers for use by both private and commercial entities. However, these advances are also likely to lead to greater vehicle densities, a more heterogenous mix of vehicles and equipment and greater levels of vehicle autonomy, which can increase the chance for communications disruptions. For the safe and secure operation of these vehicles, it is essential to have a robust communications network. This work is focused on harnessing the power of quantum technologies to enable this robust communications network by: (1) utilizing quantum optimization algorithms to design robust network with routing redundancy that can respond adaptively to dynamically changing real-time environment and disruptions, (2) utilize quantum optimization algorithms resource allocation for detection, localization, and tracking of mobile communication disruption agents and (3) utilize quantum key distribution (QKD) to execute secure key sharing in anti-jamming protocols for secure radio frequency (RF) communication. Efforts to map these quantum optimization algorithms to commercially available quantum annealers and soon to be available general-purpose gate-model quantum hardware architectures will be reviewed, and plans for testing the solutions to these algorithms through indoor sUAS flight tests will be discussed. Lastly, efforts to miniaturize and practically deploy Quantum Key Distribution (QKD) hardware, which could ultimately be used to securely exchange encryption keys, in sUAS networks will be reviewed.

Quantum Computing↗

A Deep Gravity model for mobility flows generation

The movements of individuals within and among cities influence critical aspects of our society, such as well-being, the spreading of epidemics, and the quality of the environment. When information about mobility flows is not available for a particular region of interest, we must rely on mathematical models to generate them. In this work, we propose Deep Gravity, an effective model to generate flow probabilities that exploits many features (e.g., land use, road network, transport, food, health facilities) extracted from voluntary geographic data, and uses deep neural networks to discover non-linear relationships between those features and mobility flows. Our experiments, conducted on mobility flows in England, Italy, and New York State, show that Deep Gravity achieves a significant increase in performance, especially in densely populated regions of interest, with respect to the classic gravity model and models that do not use deep neural networks or geographic data. Deep Gravity has good generalization capability, generating realistic flows also for geographic areas for which there is no data availability for training. Finally, we show how flows generated by Deep Gravity may be explained in terms of the geographic features and highlight crucial differences among the three considered countries interpreting the model’s prediction with explainable AI techniques.

97 MATHEMATICS AND COMPUTING↗

Brownian dynamics of confined suspensions of active microrollers

In this study, we develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. Finally, the characteristic height in turn controls the time scale and wavelength for the development of the fingering instability

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗