5,006 research outputs found
Diffusion-Based Coarse Graining in Hybrid Continuum-Discrete Solvers: Theoretical Formulation and A Priori Tests
Coarse graining is an important ingredient in many multi-scale
continuum-discrete solvers such as CFD--DEM (computational fluid
dynamics--discrete element method) solvers for dense particle-laden flows.
Although CFD--DEM solvers have become a mature technique that is widely used in
multiphase flow research and industrial flow simulations, a flexible and
easy-to-implement coarse graining algorithm that can work with CFD solvers of
arbitrary meshes is still lacking. In this work, we proposed a new coarse
graining algorithm for continuum--discrete solvers for dense particle-laden
flows based on solving a transient diffusion equation. Via theoretical analysis
we demonstrated that the proposed method is equivalent to the statistical
kernel method with a Gaussian kernel, but the current method is much more
straightforward to implement in CFD--DEM solvers. \textit{A priori} numerical
tests were performed to obtain the solid volume fraction fields based on given
particle distributions, the results obtained by using the proposed algorithm
were compared with those from other coarse graining methods in the literature
(e.g., the particle centroid method, the divided particle volume method, and
the two-grid formulation). The numerical tests demonstrated that the proposed
coarse graining procedure based on solving diffusion equations is theoretically
sound, easy to implement and parallelize in general CFD solvers, and has
improved mesh-convergence characteristics compared with existing coarse
graining methods. The diffusion-based coarse graining method has been
implemented into a CFD--DEM solver, the results of which are presented in a
separate work (R. Sun and H. Xiao, Diffusion-based coarse graining in hybrid
continuum-discrete solvers: Application in CFD-DEM solvers for particle laden
flows)
Diffusion-Based Coarse Graining in Hybrid Continuum-Discrete Solvers: Applications in CFD-DEM
In this work, a coarse-graining method previously proposed by the authors in
a companion paper based on solving diffusion equations is applied to CFD-DEM
simulations, where coarse graining is used to obtain solid volume fraction,
particle phase velocity, and fluid-particle interaction forces. By examining
the conservation requirements, the variables to solve diffusion equations for
in CFD-DEM simulations are identified. The algorithm is then implemented into a
CFD-DEM solver based on OpenFOAM and LAMMPS, the former being a
general-purpose, three-dimensional CFD solver based on unstructured meshes.
Numerical simulations are performed for a fluidized bed by using the CFD-DEM
solver with the diffusion-based coarse-graining algorithm. Converged results
are obtained on successively refined meshes, even for meshes with cell sizes
comparable to or smaller than the particle diameter. This is a critical
advantage of the proposed method over many existing coarse-graining methods,
and would be particularly valuable when small cells are required in part of the
CFD mesh to resolve certain flow features such as boundary layers in wall
bounded flows and shear layers in jets and wakes. Moreover, we demonstrate that
the overhead computational costs incurred by the proposed coarse-graining
procedure are a small portion of the total costs in typical CFD-DEM simulations
as long as the number of particles per cell is reasonably large, although
admittedly the computational overhead of the coarse graining often exceeds that
of the CFD solver. Other advantages of the present algorithm include more
robust and physically realistic results, flexibility and easy implementation in
almost any CFD solvers, and clear physical interpretation of the computational
parameter needed in the algorithm. In summary, the diffusion-based method is a
theoretically elegant and practically viable option for CFD-DEM simulations
High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications
A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration
Parallel finite element simulation of large ram-air parachutes
In the near future, large ram-air parachutes are expected to provide the capability of delivering 21 ton payloads from altitudes as high as 25,000 ft. In development and test and evaluation of these parachutes the size of the parachute needed and the deployment stages involved make high-performance computing (HPC) simulations a desirable alternative to costly airdrop tests. Although computational simulations based on realistic, 3D, time-dependent models will continue to be a major computational challenge, advanced finite element simulation techniques recently developed for this purpose and the execution of these techniques on HPC platforms are significant steps in the direction to meet this challenge. In this paper, two approaches for analysis of the inflation and gliding of ram-air parachutes are presented. In one of the approaches the point mass flight mechanics equations are solved with the time-varying drag and lift areas obtained from empirical data. This approach is limited to parachutes with similar configurations to those for which data are available. The other approach is 3D finite element computations based on the Navier-Stokes equations governing the airflow around the parachute canopy and Newton's law of motion governing the 3D dynamics of the canopy, with the forces acting on the canopy calculated from the simulated flow field. At the earlier stages of canopy inflation the parachute is modelled as an expanding box, whereas at the later stages, as it expands, the box transforms to a parafoil and glides. These finite element computations are carried out on the massively parallel supercomputers CRAY T3D and Thinking Machines CM-5, typically with millions of coupled, non-linear finite element equations solved simultaneously at every time step or pseudo-time step of the simulation
- …