This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid mome...This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid momentum equations are solved by using the fractional step method,while the energy equation is solved by employing the second-order Adams-Bashforth method.For efficient load balancing between the CPU cores,the coupling between particles and fluid is obtained by applying the body force in the fluid equations,which depends on the solid volume fraction of particles contained in each grid cell,and then by linearly interpolating the particle temperature and velocity on the fluid grid cell(in place of the delta function commonly used in literature).Several test cases from the literature are studied,and good agreement is observed between the simulation results and the literature.Finally,a scaling study on multiple core machines is performed,demonstrating the proposed IBM's capabilities for a significant reduction in processing time.展开更多
The directional drifting of particles/molecules with broken symmetry has received increasing attention. Through molecular dynamics simulations, we investigate the effects of various solvents on the time-dependent dire...The directional drifting of particles/molecules with broken symmetry has received increasing attention. Through molecular dynamics simulations, we investigate the effects of various solvents on the time-dependent directional drifting of a particle with broken symmetry. Our simulations show that the distance of directional drift of the asymmetrical particle is reduced while the ratio of the drift to the mean displacement of the particle is enhanced with increasing mass, size, and interaction strength of the solvent atoms in a short time range. Among the parameters considered, solvent atom size is a particularly influential factor for enhancing the directional drift of asymmetrical particles, while the effects of the interaction strength and the mass of the solvent atoms are relatively weaker. These findings are of great importance to the understanding and control of the Brownian motion of particles in various physical, chemical, and biological processes within finite time spans.展开更多
Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O...Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O(1000)cores due to the cost of solving at least one large three-dimensional pressure Poisson problem per time step.The parallelization may perform moderately well only or even poorly sometimes despite using an efficient algebraic multigrid preconditioner[38].We present an accurate and scalable solver using a direction splitting algorithm[12]to transform all three-dimensional parabolic/elliptic problems(and in particular the elliptic pressure Poisson problem)into a sequence of three one-dimensional parabolic sub-problems,thus improving its scalability up to multiple thousands of cores.We employ this algorithm to solve mass and momentum conservation equations in flows laden with fixed non-spherical rigid bodies.We consider the presence of rigid bodies on the(uniform or non-uniform)fixed Cartesian fluid grid by modifying the diffusion and divergence stencils on the impacted grid node near the rigid body boundary.Compared to[12],we use a higher-order interpolation scheme for the velocity field to maintain a secondorder stress estimation on the particle boundary,resulting in more accurate dimensionless coefficients such as drag C_(d)and lift C_(l).We also correct the interpolation scheme due to the presence of any nearby particle to maintain an acceptable accuracy,making the solver robust even when particles are densely packed in a sub-region of the computational domain.We present classical validation tests involving a single or multiple(up to O(1000))rigid bodies and assess the robustness,accuracy and computing speed of the solver.We further show that the Direction Splitting solver is∼5 times faster on 5120 cores than our solver[38]based on a classical projection method[5].展开更多
A direct drive actuator (DDA) with direct drive valves (DDVs) as the control device is an ideal solution for a flight actuation system. This paper presents a novel triple-redundant voice coil motor (TRVCM) used ...A direct drive actuator (DDA) with direct drive valves (DDVs) as the control device is an ideal solution for a flight actuation system. This paper presents a novel triple-redundant voice coil motor (TRVCM) used for redundant DDVs. The TRVCM features electrical/mechanical hybrid triple-redundancy by securing three stators along with three moving coils in the same frame. A permanent magnet (PM) Halbach array is employed in each redundant VCM to simplify the system structure. A back-to-back design between neighborly redundancies is adopted to decouple the magnetic flux linkage. The particle swarm optimization (PSO) method is implemented to optimize design parameters based on the analytical magnetic circuit model. The optimization objective function is defined as the acceleration capacity of the motor to achieve high dynamic performance. The optimal geometric parameters are verified with 3D magnetic field finite element analysis (FEA). A research prototype has been developed for experimental purpose. The experimental results of magnetic field density and force output show that the proposed TRVCM has great potential of applications in DDA systems.展开更多
文摘This paper investigates an efficient immersed boundary method(IBM)for multiple-core CPU machines with local grid refinement for the calculation of heat transfer between fluids and finite-sized particles.The fluid momentum equations are solved by using the fractional step method,while the energy equation is solved by employing the second-order Adams-Bashforth method.For efficient load balancing between the CPU cores,the coupling between particles and fluid is obtained by applying the body force in the fluid equations,which depends on the solid volume fraction of particles contained in each grid cell,and then by linearly interpolating the particle temperature and velocity on the fluid grid cell(in place of the delta function commonly used in literature).Several test cases from the literature are studied,and good agreement is observed between the simulation results and the literature.Finally,a scaling study on multiple core machines is performed,demonstrating the proposed IBM's capabilities for a significant reduction in processing time.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11175230 and 11474299)Innovation Program of the Shanghai Municipality Education Commission, China (Grant No. 14ZZ095)Shanghai Supercomputer Center and Supercomputing Center of the Chinese Academy of Sciences
文摘The directional drifting of particles/molecules with broken symmetry has received increasing attention. Through molecular dynamics simulations, we investigate the effects of various solvents on the time-dependent directional drifting of a particle with broken symmetry. Our simulations show that the distance of directional drift of the asymmetrical particle is reduced while the ratio of the drift to the mean displacement of the particle is enhanced with increasing mass, size, and interaction strength of the solvent atoms in a short time range. Among the parameters considered, solvent atom size is a particularly influential factor for enhancing the directional drift of asymmetrical particles, while the effects of the interaction strength and the mass of the solvent atoms are relatively weaker. These findings are of great importance to the understanding and control of the Brownian motion of particles in various physical, chemical, and biological processes within finite time spans.
基金support of the University of British Columbia via its Four Year Doctoral Fellowship programThe authors greatly appreciate the financial support of the Natural Sciences and Engineering Research Council of Canada(NSERC)via Anthony Wachs’s Discovery Grant RGPIN-2016-06572+1 种基金This research was enabled by the support provided by Compute Canada(http://www.computecanada.ca)through Anthony Wachs’s 2020,2021,and 2022 Resources for Research Groups allocation qpf-764-abThis research was also supported in part through computational resources and services provided by Advanced Research Computing at the University of British Columbia.
文摘Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O(1000)cores due to the cost of solving at least one large three-dimensional pressure Poisson problem per time step.The parallelization may perform moderately well only or even poorly sometimes despite using an efficient algebraic multigrid preconditioner[38].We present an accurate and scalable solver using a direction splitting algorithm[12]to transform all three-dimensional parabolic/elliptic problems(and in particular the elliptic pressure Poisson problem)into a sequence of three one-dimensional parabolic sub-problems,thus improving its scalability up to multiple thousands of cores.We employ this algorithm to solve mass and momentum conservation equations in flows laden with fixed non-spherical rigid bodies.We consider the presence of rigid bodies on the(uniform or non-uniform)fixed Cartesian fluid grid by modifying the diffusion and divergence stencils on the impacted grid node near the rigid body boundary.Compared to[12],we use a higher-order interpolation scheme for the velocity field to maintain a secondorder stress estimation on the particle boundary,resulting in more accurate dimensionless coefficients such as drag C_(d)and lift C_(l).We also correct the interpolation scheme due to the presence of any nearby particle to maintain an acceptable accuracy,making the solver robust even when particles are densely packed in a sub-region of the computational domain.We present classical validation tests involving a single or multiple(up to O(1000))rigid bodies and assess the robustness,accuracy and computing speed of the solver.We further show that the Direction Splitting solver is∼5 times faster on 5120 cores than our solver[38]based on a classical projection method[5].
基金supported by National Science Foundation for Distinguished Young Scholars of China(No.50825502)National Natural Science Foundation of China(No.51105016)
文摘A direct drive actuator (DDA) with direct drive valves (DDVs) as the control device is an ideal solution for a flight actuation system. This paper presents a novel triple-redundant voice coil motor (TRVCM) used for redundant DDVs. The TRVCM features electrical/mechanical hybrid triple-redundancy by securing three stators along with three moving coils in the same frame. A permanent magnet (PM) Halbach array is employed in each redundant VCM to simplify the system structure. A back-to-back design between neighborly redundancies is adopted to decouple the magnetic flux linkage. The particle swarm optimization (PSO) method is implemented to optimize design parameters based on the analytical magnetic circuit model. The optimization objective function is defined as the acceleration capacity of the motor to achieve high dynamic performance. The optimal geometric parameters are verified with 3D magnetic field finite element analysis (FEA). A research prototype has been developed for experimental purpose. The experimental results of magnetic field density and force output show that the proposed TRVCM has great potential of applications in DDA systems.