SEARCH · Engineering Papers

Results for “Conservation equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2022

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Path integral games with de Sitter α-vacua

The α-vacua are a 1-parameter family of quantum field vacua in de Sitter space which are invariant under the isometry group SO(1, d). In this work give a path integral construction of the de Sitter α-vacua. We explain that these states can be prepared by acting on the Bunch-Davies vacuum with a certain non-local charge operator. While most conserved charges live on a single codimension-1 manifold, we show that this particular charge lives on a pair of two codimension-1 manifolds which are antipodal mirrors of each other. The rules for the manipulation of this charge as an insertion in the path integral are explained. We further explain how this charge can be used to solve for the wavefunctionals of the α-vacua at $\mathcal{I}$ | (in the regime that α is small) by deforming the equator of de Sitter space to $\mathcal{I}$ + /$\mathcal{I}$ – .

Global Symmetries↗

DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Lepton flavor asymmetries: from the early Universe to BBN

Large primordial lepton flavor asymmetries with almost vanishing total baryon-minus-lepton number can evade the usual BBN and CMB constraints if neutrino oscillations lead to perfect flavor equilibration. Solving the momentum averaged quantum kinetic equations (QKEs) describing neutrino oscillations and interactions, we perform the first systematic investigation of this scenario, uncovering a rich flavor structure in stark contradiction to the assumption of simple flavor equilibration. We find (i) a particular direction in flavor space, ∆ne ≃ – 2/3 (– 1)∆n μ for normal (inverted) neutrino mass hierarchy, in which the flavor equilibration is efficient and primordial asymmetries are essentially unconstrained, (ii) a minimal washout factor, ∆$n_{e}^{2}$| BBN ≤ 0.03 (0.016) ∑ α ∆$n_{α}^{2}$| ini yielding a conservative estimate for the allowed primordial asymmetries in a generic flavor direction, and (iii) particularly strong or weak washout if one of the initial flavor asymmetries vanishes due to non-adiabatic muon- or electron-driven MSW transitions. These results open up the possibility of a first-order QCD phase transition facilitated by large lepton asymmetries as well as baryogenesis from large and compensated ∆n e = ∆n μ asymmetries. Our systematic approach of deriving momentum averaged QKEs includes collision terms beyond the damping approximation, energy transfer between the neutrino and electron-photon plasma, and provides a fast and reliable way to investigate the impact of primordial lepton asymmetries at the time of BBN. We publicly release the Mathematica code COFLASY-M on https://github.com/mariofnavarro/COFLASY which solves the QKEs numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Using Likwid and Byfl to Benchmark Hardware Performance

This paper outlines a benchmarking study conducted during my internship at LANL, focusing on CPU (Computer Processing Unit) and program performance assessment. The primary goal was to gather memory access data using three methods across five polybench kernels The data gathered would then be used to compare and contrast to one another and calculate operational intensity for performance comparisons. Benchmarking tools like Byfl and Likwid were employed, with Byfl offering hardware-independent data through LLVM compiler communication and Likwid directly interacting with computer hardware. The study considered various benchmarking factors, including optimization levels, Big O notation ((n)), CPU diversity and specific kernel equations. Big O notation was utilized to simplify code complexity, with detailed breakdwons of operations and memory components for each polybench application. Specific O(n) equations enabled nuanced kernel compariosns, facilitating the identification of performance variations. CPU efficiency assessments were conducted using Likwid tests on two CPUs. The central focus on code optimization aimed at achieving higher speeds and reduced memory usage through streamlined code. Future work propsoes creating a roofline model, synthesizing benchmarking data into a comprehensive data graph to assist in optimizing code and improving hardware performance. The potential impact on the laboratory or national mission was underscored, emphasizing the importance of optimizing applications and hardware to conserve resources and accelerate program execution. The specific relevance to LANL’s operations in math-intensive fields such as Nuclear Fission, Space Exploration, and Nanotechnology highlights the necessity of efficient benchmarking for resource conservation and proram speed. Overall, this study contributes to the understanding of CPU and program performance, providing insights for future optimization efforts in a laboratory setting

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

CurvilinearGrids.jl: A Julia package for curvilinear coordinate transformations

Finite-difference discretizations of partial differential equations are widespread throughout the scientific community. Oftentimes finite-differences are used to compute spatial gradients of fields on a discrete grid, which is typically a uniform or rectilinear Cartesian mesh. Arbitrary multidimensional geometry is difficult to discretize directly with finite differences, however, due to non-uniform grid spacing and non-orthogonality. Curvilinear coordinate transformations can be used as a strategy to enable arbitrary geometry. While these curvilinear transformations are straightforward, the governing PDEs require additional terms (metrics) and must adhere to strict conservation laws; these criteria complicate the application of the transformation and require careful implementation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Energy Intensity Baselining and Tracking Guidance

Each company joining the U.S. Department of Energy’s (DOE’s) Better Buildings, Better Plants Program (Better Plants) commits to establishing an energy consumption and energy intensity (EI) baseline and to tracking its energy performance over a 10-year period against that baseline. The baseline must reflect a company’s energy consumption over a 12-month period, covering all its U.S.-based operations. Energy consumption is calculated by fuel type in terms of primary energy (also known as source energy). EI is broadly defined as the amount of energy consumed per unit of output produced. For this guidance document and for the program, the term energy performance represents an evaluation of a facility’s capacity to use energy efficiently. Metrics used to assess a facility’s energy performance can include EI, energy consumption, improvements in EI, etc. Establishing an energy baseline and tracking system is a critical first step in effectively managing energy use. Developing a baseline can help a company understand energy use within the corporation and give it a point of comparison to evaluate future efforts to improve energy performance. It can also support efforts to validate a company’s energy management activities, improve comparative analyses when using benchmarks, and help in predicting future energy needs. In addition, a company that normalizes its performance data can determine highly defensible measures of energy savings generated through implemented energy efficiency projects. Establishing a baseline and tracking energy performance is also a requirement for ISO 50001 certification. Although basic energy data can be collected through utility bills, most manufacturers will have to perform additional analyses to develop accurate and robust energy baselines and tracking systems. Energy is consumed in many ways within the manufacturing sector and can come from multiple sources. Energy is sometimes generated and sold to other parties or captured and reused on-site. External events can exert a significant impact on a facility or company’s energy use independent of any purposeful efforts to improve energy efficiency. Operational changes, such as production shifts—which may be inevitable for some companies over the 10-year period covered by the program—can also make a big difference in energy use. Since Better Plants asks companies to account for all their U.S.-based operations, mergers, acquisitions, and divestitures can also have significant implications for a company’s energy metrics. This document aims to demystify the sometimes complex baselining process. It devotes special attention to the task of normalizing and adjusting energy consumption to account for external factors, such as weather and production changes. A key recommendation is that companies use regression analysis to normalize their energy consumption data whenever possible. Regression analysis is a statistical technique that estimates the dependence of a variable (i.e., energy use in the context of Better Plants) on one or more independent variables such as ambient temperature, while controlling for the influence of other variables at the same time. A properly developed regression analysis can provide a reliable estimate of energy savings resulting from energy improvement strategies and projects by accounting for the effects of variables such as annual production levels and weather. DOE has developed a companion Energy Performance Indicator software tool (EnPI) to simplify the baselining process. This tool can run regression models, calculate changes in EI at the facility level, and automatically compile facility-level data into a corporate-wide metric. Note that although the relevant equations used to calculate EI are provided in this document, the EnPI tool will automatically perform most calculations for the user. Additionally, Better Plants Partners (Partners) can call on their Technical Account Manager (TAM) to help them establish a baseline and assist with the necessary calculations to track progress.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

The theory of transport in helical spin-structure crystals

Abstract We study helical structures in spin-spiral single crystals. In the continuum approach for the helicity potential energy the simple electronic band splits into two non-parabolic bands. For low exchange integrals, the lower band is described by a surface with a saddle shape in the direction of the helicity axis. Using the Boltzmann equation with the relaxation due to acoustic phonons, we discover the dependence of the current on the angle between the electric field and helicity axis leading to the both parallel and perpendicular to the electric field components in the electroconductivity. The latter can be interpreted as a planar Hall effect. In addition, we find that the transition rates depend on an electron spin allowing the transition between the bands. The electric conductivities exhibit nonlinear behaviors with respect to chemical potential µ . We explain this effect as the interference of the band anisotropy, spin conservation, and interband transitions. The proposed theory with the spherical model in the effective mass approximation for conduction electrons can elucidate nonlinear dependencies that can be identified in experiments. We find the excellent agreement between the theoretical and experimental data for parallel resistivity depending on temperature at the phase transition from helical to ferromagnetic state in a M n P single crystal. In addition, we predict that the perpendicular resistivity abruptly drops to zero in the ferromagnetic phase.

Physics↗

DOE OSTI2022

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Machine-assisted discovery of integrable symplectic mappings

Integrable systems possess a hidden symmetry associated with the existence of conserved quantities known as integrals of motion. These systems play an important role in understanding general dynamics in accelerators and have potential for future designs. This work will cover two automated methods for finding integrable symplectic maps of the plane. The first algorithm is based on the observation that the evolution of an integrable system in phase space is confined to a lower-dimensional submanifold of a specific type. The second algorithm relies on an analysis of dynamical variables. Both methods rediscover some of the famous McMillan-Suris integrable mappings and ultra-discrete Painlev\'e equations. Over 100 new integrable families are presented and analyzed, some of which are isolated in the space of parameters, while others are families with one parameter (or the ratio of parameters) being either continuous or discrete. In addition, the newly discovered maps are related to a general 2D symplectic map through the use of discrete perturbation theory. A method is proposed for constructing smooth near-integrable dynamical systems based on mappings with polygon invariants.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Machine-Assisted Discovery of Integrable Symplectic Mappings

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

"Hall" transport of liquid crystal solitons in Couette flow

Topology establishes a unifying framework for a diverse range of scientific areas including particle physics, cosmology, and condensed matter physics. One of the most fascinating manifestations of topology in the context of condensed matter is the topological Hall effect, and its relative: the Skyrmion Hall effect. Skyrmions are stable vortex-like spin configurations in certain chiral magnets, and when subject to external electric currents can drift in the transverse direction to the current. These quasi-particles are characterised by a conserved topological charge which in the Skyrmion Hall effect plays the role of electric charges in the ordinary Hall effect. Recently, it has been shown that liquid crystals endowed with chiral properties serve as an ideal testbed for the fundamental investigation of topological solitons, including their two- and three-dimensional realisations. Here, we show experimentally and numerically that three-dimensional solitons aka "torons" exhibit a Hall-like effect when driven by shear flows: the torons are deflected in the direction perpendicular to the shear plane. The experimental results are rationalised in terms of the dynamic Ericksen-Leslie equations, which predict the emergence of the transverse component of the net mass flow, the magnitude of which scales as the 3rd power of the shear rate. The perturbation analysis highlights an interplay of the viscous and chiral elastic torques as the mechanism for the emergence of net transverse currents. Numerical simulations demonstrate, however, that torons are not merely dragged by the flow but move with their own transverse speed, much larger than the average flow velocity in the transverse direction. Our findings may enable responsive microfluidic applications relying on soft topological solitons.

FOS: Physical sciences↗

DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗

DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗