Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluid”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Code Comparison Study Fosters Confidence in the Numerical Simulation of Enhanced Geothermal Systems

Numerical simulation has become a standard analytical tool for scientists and engineers to evaluate the potential and performance of enhanced geothermal systems. A variety of numerical simulators developed by industry, universities, and national laboratories are currently available and being applied to better understand enhanced geothermal systems at the field scale. To yield credible predictions and be of value to site operators, numerical simulators must be able to accurately represent the complex coupled processes induced by producing geothermal systems, such as fracture aperture changes due to thermal stimulation, fracture shear displacement with fluid injection, rate of thermal depletion of reservoir rocks, and permeability alteration with mineral precipitation or dissolution. A suite of numerical simulators was exercised on a series of test problems that considered coupled thermal, hydraulic, geomechanical, and geochemical (THMC) processes. Problems were selected and designed to isolate selected coupled processes, to be executed on workstation class computers, and have simple but illustrative metrics for result comparisons. This paper summarizes the initial suite of seven benchmark problems, describes the code comparison activities, provides example results for problems and documents the capabilities of currently available numerical simulation codes to represent coupled processes that occur during the production of geothermal resources. Code comparisons described in this paper use the ISO (International Organization for Standardization) standard ISO-13538 for proficiency testing of numerical simulators. This approach was adopted for a recent code comparison study within the radiation transfer-modeling field of atmospheric sciences, which was focused on canopy reflectance models. This standard specifies statistical methods for analyzing laboratory data from proficiency testing schemes to demonstrate that the measurement results do not exhibit evidence of an unacceptable level of bias.

benchmark problems↗
DOE OSTI2013

Modeling Solar Wind Flow with the Multi-Scale Fluid-Kinetic Simulation Suite

Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. We have enhanced the code with additional physical treatments for the transport of turbulence and acceleration of pickup ions in the interplanetary space and at the termination shock. In this article, we present the results of our numerical simulation of the solar wind (SW) interaction with the local interstellar medium (LISM) in different time-dependent and stationary formulations. Numerical results are compared with the Ulysses, Voyager, and OMNI observations. Finally, the SW boundary conditions are derived from in-situ spacecraft measurements and remote observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Declaration of the Chemical Composition of Waste-products for the German Repository Konrad Using Substance Vectors - 20500

The waste acceptance criteria of the German repository Konrad require beside a radiological declaration a compositional declaration to be able to avoid negative impact on the groundwater by a presumed leakage into the groundwater after complete corrosion of the containers in the future. The declaration of the composition has to specify compounds and substances not element analyses, because the chemical behavior and toxicity especially of organic compounds varies strongly The relevant compounds have to be described and quantified using defined substance codes provided by the operator. For simple waste types like pure concrete or steel components, the definition of the composition is rather simple. The composition of more complex waste streams is generally declared using producer specific material vectors, which can be derived for defined waste types (e.g. grouted fluids, ash from combustion, super-compacted mixed waste etc.). Additional substance vectors exist for containers, drums and grout, which together with the waste composition yield the total composition of the waste packages. The substance vectors for homogeneous waste streams like cementitious demolition waste, solidified fluids and ash are mostly based on analyses provided by the waste owner. In some cases, analyses from non-radioactive analogues or data sheets can also be used to demonstrate the composition. For each qualified treatment of radioactive waste streams the material composition is described and balanced on the basis of specific material vectors, which are composed of material codes from the material and cask list provided by the Federal company for radioactive waste disposal. For mixed waste streams like super-compacted or grouted mixed waste or materials simply packed in containers an analyses cannot be used. In these cases, the compositional information from the bookkeeping can be used for the entire waste stream to calculate a mean composition including a compositional range (min / max composition). Only waste packages with an approved material description can be delivered to Konrad. Therefore, the material declaration of all waste streams and the resulting waste products is crucial for all waste owners in Germany to be able to have enough waste packages ready, when the repository is going to open in 2027. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2001

Adjustable Autonomy Testbed

The Adjustable Autonomy Testbed (AAT) is a simulation-based testbed located in the Intelligent Systems Laboratory in the Automation, Robotics and Simulation Division at NASA Johnson Space Center. The purpose of the testbed is to support evaluation and validation of prototypes of adjustable autonomous agent software for control and fault management for complex systems. The AA T project has developed prototype adjustable autonomous agent software and human interfaces for cooperative fault management. This software builds on current autonomous agent technology by altering the architecture, components and interfaces for effective teamwork between autonomous systems and human experts. Autonomous agents include a planner, flexible executive, low level control and deductive model-based fault isolation. Adjustable autonomy is intended to increase the flexibility and effectiveness of fault management with an autonomous system. The test domain for this work is control of advanced life support systems for habitats for planetary exploration. The CONFIG hybrid discrete event simulation environment provides flexible and dynamically reconfigurable models of the behavior of components and fluids in the life support systems. Both discrete event and continuous (discrete time) simulation are supported, and flows and pressures are computed globally. This provides fast dynamic simulations of interacting hardware systems in closed loops that can be reconfigured during operations scenarios, producing complex cascading effects of operations and failures. Current object-oriented model libraries support modeling of fluid systems, and models have been developed of physico-chemical and biological subsystems for processing advanced life support gases. In FY01, water recovery system models will be developed.

Malin, Jane T.↗
NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗
DOE OSTI2022

Latent Space Simulation for Carbon Capture Design Optimization

The CO2 capture efficiency in solvent-based carbon capture systems (CCSs) critically depends on the gas-solvent interfacial area (IA), making maximization of IA a foundational challenge in CCS design. While the IA associated with a particular CCS design can be estimated via a computational fluid dynamics (CFD) simulation, using CFD to derive the IAs associated with numerous CCS designs is prohibitively costly. However, previous works such as Deep Fluids (Kim et al., 2019) show that large simulation speed-ups and low error are achievable by replacing CFD simulators with neural-network (NN) surrogates that mimic the CFD simulation process. This raises the possibility for a fast, accurate replacement for a CFD simulator, and thus computationally feasible IA-based CCS-design optimization. As such, here, we explore whether an existing NN-surrogate approach (and variants we develop) can successfully be applied to our complex carbon-capture CFD simulations, with the ultimate goal of obtaining a fast and accurate simulator for our CCS-design application. Our experiments build on the Deep Fluids approach and find that resulting surrogates can produce large speed ups (4000x) while maintaining IA relative errors as low as 4% on unseen CCS configurations (interpolating between configurations seen during training). Thus, despite less faithfulness to the underlying physics of the problem, NN surrogates may be a promising tool for our CCS design optimization problem. Notably, though, the Deep Fluids approach has limitations for our application (such as model non-transferability to CCS-packing changes), which we discuss. We conclude with potential directions for future work that may, like innovations we introduced here (e.g., transformer-based dynamics prediction), improve performance on our complicated dataset of CCS CFD simulations.

AI Surrogate Simulation, Carbon Capture, design op↗
DOE OSTI2023

Reducing the Flow Maldistribution in Heat Exchangers through a Novel Polymer Manifold: Numerical Evaluation

The maldistribution of working fluid is one of the issues in heat exchangers that causes a reduction in performance of not only the heat exchanger but also the entire HVAC system. One of the methods to reduce such maldistribution is to improve manifold design to evenly distribute the flow. In the present work, an advanced maldistribution reduction manifold, which was based on a preliminary maldistribution reduction manifold, was designed to further improve the flow distribution in the heat exchanger. In the design, spiral baffles are used to create vortices in the tubes to regulate the flow in each tube. The design also keeps the tubes away from the manifold inlet to avoid direct flow from the inlet. Due to the complexity, the design of the advanced maldistribution reduction manifold is for AM only, which cannot be fabricated by traditional manufacturing. To evaluate the design, a computational fluid dynamic model is developed to study flow distribution in heat exchanger manifolds. The simulation results reveal that the relative standard deviation of the tubes in the advanced maldistribution reduction design is half of the preliminary maldistribution reduction design and about 1/20 of the reference design.

42 ENGINEERING↗
DOE OSTI2025

In situ Raman investigation of Dy complexation in Cl-bearing aqueous solutions at 20–300 °C

Raman spectroscopy provides a versatile tool for in situ characterization of aqueous rare earth elements (REE) speciation at the molecular level. Complexation of REE with ligands such as Cl – and OH – is of particular interest for understanding the mobility of REE in NaCl-bearing hydrothermal fluids responsible for enriching REE to economic levels in nature. Raman spectroscopic studies of REE speciation in Cl-bearing aqueous fluids are primarily conducted at ambient temperature, whereas natural systems indicate temperatures of >100–600 °C. In this study, the speciation of Dy in acidic chloride-bearing hydrothermal solutions was investigated using confocal Raman spectroscopy with a new capillary Raman heating stage at 20–300 °C. Background solutions (pure water, NaCl-solutions) and solutions with 0.14–1.8 mol kg –1 dissolved DyCl 3 were sealed in quartz capillary cells. Comparison of the spectra for Dy chloride solutions with those for background solutions and the spectra for reference Dy-bearing solids was used to identify Raman bands specific to Dy–O and Dy–Cl bonds. The Raman band for the Dy–O stretching mode of hydrated Dy 3+ aqua ions was measured at 365–384 cm –1 and a Raman band for the Dy–Cl stretching modes of Dy chloride complexes was measured near 240 cm –1 . The Dy–O band decreases systematically with temperature, whereas the Dy–Cl band systematically increases, indicating a systematic increase in the stability of Dy chloride complexes with temperature. Here, this study provides the framework for expanding the use of in situ Raman spectroscopy to investigate the speciation of REE in aqueous solutions to hydrothermal conditions.

58 GEOSCIENCES↗
NASA NTRS1995

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗
DOE OSTI2014

Discrete models of fluids: spatial averaging, closure and model reduction

We consider semidiscrete ODE models of single-phase fluids and two-fluid mixtures. In the presence of multiple fine-scale heterogeneities, the size of these ODE systems can be very large. Spatial averaging is then a useful tool for reducing computational complexity of the problem. The averages satisfy exact balance equations of mass, momentum, and energy. These equations do not form a satisfactory continuum model because evaluation of stress and heat flux requires solving the underlying ODEs. To produce continuum equations that can be simulated without resolving microscale dynamics, we recently proposed a closure method based on the use of regularized deconvolution. Here we continue the investigation of deconvolution closure with the long term objective of developing consistent computational upscaling for multiphase particle methods. The structure of the fine-scale particle solvers is reminiscent of molecular dynamics. For this reason we use nonlinear averaging introduced for atomistic systems by Noll, Hardy, and Murdoch-Bedeaux. We also consider a simpler linear averaging originally developed in large eddy simulation of turbulence. We present several simple but representative examples of spatially averaged ODEs, where the closure error can be analyzed. Based on this analysis we suggest a general strategy for reducing the relative error of approximate closure. For problems with periodic highly oscillatory material parameters we propose a spectral boosting technique that augments the standard deconvolution and helps to correctly account for dispersion effects. We also conduct several numerical experiments, one of which is a complete mesoscale simulation of a stratified two-fluid flow in a channel. In this simulation, the operation count per coarse time step scales sublinearly with the number of particles.

Multiscale modeling, discrete models, fluid flow↗
DOE OSTI2017

Pressure evolution equation for the particulate phase in inhomogeneous compressible disperse multiphase flows

An analytical expression describing the unsteady pressure evolution of the dispersed phase driven by variations in the carrier phase is presented. In this article, the term “dispersed phase” represents rigid particles, droplets, or bubbles. Letting both the dispersed and continuous phases be inhomogeneous, unsteady, and compressible, the developed pressure equation describes the particle response and its eventual equilibration with that of the carrier fluid. The study involves impingement of a plane traveling wave of a given frequency and subsequent volume-averaged particle pressure calculation due to a single wave. The ambient or continuous fluid's pressure and density-weighted normal velocity are identified as the source terms governing the particle pressure. Analogous to the generalized Faxén theorem, which is applicable to the particle equation of motion, the pressure expression is also written in terms of the surface average of time-varying incoming flow properties. The surface average allows the current formulation to be generalized for any complex incident flow, including situations where the particle size is comparable to that of the incoming flow. Further, the particle pressure is also found to depend on the dispersed-to-continuous fluid density ratio and speed of sound ratio in addition to dynamic viscosities of both fluids. The model is applied to predict the unsteady pressure variation inside an aluminum particle subjected to normal shock waves. The results are compared against numerical simulations and found to be in good agreement. Furthermore, it is shown that, although the analysis is conducted in the limit of negligible flow Reynolds and Mach numbers, it can be used to compute the density and volume of the dispersed phase to reasonable accuracy. Finally, analogous to the pressure evolution expression, an equation describing the time-dependent particle radius is deduced and is shown to reduce to the Rayleigh-Plesset equation in the linear limit.

42 ENGINEERING↗
NASA NTRS1982

Computational fluid dynamics - The coming revolution

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Graves, R. A., Jr.↗
DOE OSTI2020

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗
DOE OSTI2015

On the rogue waves propagation in non-Maxwellian complex space plasmas

The implications of the non-Maxwellian electron distributions (nonthermal/or suprathermal/or nonextensive distributions) are examined on the dust-ion acoustic (DIA) rogue/freak waves in a dusty warm plasma. Using a reductive perturbation technique, the basic set of fluid equations is reduced to a nonlinear Schrödinger equation. The latter is used to study the nonlinear evolution of modulationally unstable DIA wavepackets and to describe the rogue waves (RWs) propagation. Rogue waves are large-amplitude short-lived wave groups, routinely observed in space plasmas. The possible region for the rogue waves to exist is defined precisely for typical parameters of space plasmas. It is shown that the RWs strengthen for decreasing plasma nonthermality and increasing superthermality. For nonextensive electrons, the RWs amplitude exhibits a bit more complex behavior, depending on the entropic index q. Moreover, our numerical results reveal that the RWs exist with all values of the ion-to-electron temperature ratio σ for nonthermal and superthermal distributions and there is no limitation for the freak waves to propagate in both two distributions in the present plasma system. But, for nonextensive electron distribution, the bright- and dark-type waves can propagate in this case, which means that there is a limitation for the existence of freak waves. Our systematic investigation should be useful in understanding the properties of DIA solitary waves that may occur in non-Maxwellian space plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗
DOE OSTI2016

Conventional and Advanced Separations in Mass Spectrometry-Based Metabolomics: Methodologies and Applications

Metabolomics is the quantitative analysis of all metabolites in a given sample. Due to the chemical complexity of the metabolome, optimal separations are required for comprehensive identification and quantification of sample constituents. This chapter provides an overview of both conventional and advanced separations methods in practice for reducing the complexity of metabolite extracts delivered to the mass spectrometer detector, and covers gas chromatography (GC), liquid chromatography (LC), capillary electrophoresis (CE), supercritical fluid chromatography (SFC) and ion mobility spectrometry (IMS) separation techniques coupled with mass spectrometry (MS) as both uni-dimensional and as multi-dimensional approaches.

capillary electrophoresis – mass spectrometry (CE-↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2014

Hybrid immersed interface-immersed boundary methods for AC dielectrophoresis

Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface–immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of the hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗