3D Radiative Transfer Equation Coupled with Heat Conduction Equation with Realistic Boundary Conditions Applied on Complex Geometries

This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem is considered steady-state but solved at each time iteration of the transient conduction problem. The discrete ordinate method along with the decentered streamline-upwind Petrov-Galerkin method is developed. Since specular reflection is considered on borders, a very accurate algorithm has been developed for calculation of partition ratio coefficients of incident solid angles to the several reflected solid angles. The developed algorithms are tested on a paraboloid-shaped geometry used for example on concentrated solar power technologies.

Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Simulation of incompressible multiphase flows with complex geometry using etching multiblock method

The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.

On Integrable Conditions of Generalized Almost Complex Structures

Generalized complex geometry is a new kind of geometrical structure which contains complex and symplectic geometry as its special cases. This paper gives the equivalence between the integrable conditions of a generalized almost complex structure in big bracket formalism and those in the general framework.

Configuration Optimization of Two-Rail Sliders on Dynamic Characteristic

An optimal design methodology for the configuration of two rail slider was proposed to get better dynamic performance. The taper length, taper height and the rail width of the reading/writing head are considered as design variables. The complex geometry method is utilized as the search scheme in the optimization process. Optimization results show that the new slider has better dynamic characteristics and is more stable than the original designed slider. The optimization process also demonstrates that the optimum model and optimum method is effective.

Φ-INEQUALITIES AND LAWS OF LARGE NUMBERS OF HARDY MARTINGALETRANSFORMS

Using stopping time method we proved the Φ-inequalities, pointwise convergence, strong and weak laws of large numbers of Hardy martingale transforms with values in complex Banach spaces, and applying them to give several characterizations of AUMD spaces.

3D NUMERICAL SOLUTION OF AERO-NOISE WITH HIGH-ORDER DISCONTINUOUS GALERKIN METHOD

The flow-induced noise is simulated with a hybrid method.Firstly,a steady-state background flow field is given by solving Reynolds averaged Navier-Stokes(RANS)equations with finite volume(FV)method on structured grid.Then the linearized Euler equations(LEE)can be constructed based on the resulted background flow field,where the source term on the right hand side is computed using stochastic noise generation and radiation(SNGR)method.Finally,the unsteady acoustic field is obtained through solving LEE using high-order discontinuous Galerkin(DG)method on unstructured grid,where the parallel computing based on mesh partitioning and a″Quadrature-Free Implementation″method for high-order DG are employed to accelerate the computation.In order to demonstrate the sound propagation in detail,a visualization method for high-order schemes is also developed here.Moreover,in order to test the validation and the accuracy,a 3D cavity test in comparison with the experimental data is displayed first in this paper,then a 3D high-lift wing is also simulated to demonstrate its capability for very complex geometries.

Computational implementation of a GIS developed tool for prediction of dynamic ground movement and deformation due to underground extraction sequence

In the last century, there has been a significant development in the evaluation of methods to predict ground movement due to underground extraction. Some remarkable developments in three-dimensional computational methods have been supported in civil engineering, subsidence engineering and mining engineering practice. However, ground movement problem due to mining extraction sequence is effectively four dimensional （4D）. A rational prediction is getting more and more important for long-term underground mining planning. Hence, computer-based analytical methods that realistically simulate spatially distributed time-dependent ground movement process are needed for the reliable long-term underground mining planning to minimize the surface environmental damages. In this research, a new computational system is developed to simulate four-dimensional （4D） ground movement by combining a stochastic medium theory, Knothe time-delay model and geographic information system （GIS） technology. All the calculations are implemented by a computational program, in which the components of GIS are used to fulfill the spatial-temporal analysis model. In this paper a tight coupling strategy based on component object model of GIS technology is used to overcome the problems of complex three-dimensional extraction model and spatial data integration. Moreover, the implementation of computational of the interfaces of the developed tool is described. The GIS based developed tool is validated by two study cases. The developed computational tool and models are achieved within the GIS system so the effective and efficient calculation methodology can be obtained, so the simulation problems of 4D ground movement due to underground mining extraction sequence can be solved by implementation of the developed tool in GIS.

Three-Dimensional Meshfree Analysis of Interlocking Concrete Blocks for Step Seawall Structure

This study adapts the flexible characteristic of meshfree method in analyzing three-dimensional(3D)complex geometry structures,which are the interlocking concrete blocks of step seawall.The elastostatic behavior of the block is analysed by solving the Galerkin weak form formulation over local support domain.The 3D moving least square(MLS)approximation is applied to build the interpolation functions of unknowns.The pre-defined number of nodes in an integration domain ranging from 10 to 60 nodes is also investigated for their effect on the studied results.The accuracy and efficiency of the studied method on 3D elastostatic responses are validated through the comparison with the solutions of standard finite element method(FEM)using linear shape functions on tetrahedral elements and the well-known commercial software,ANSYS.The results show that elastostatic responses of studied concrete block obtained by meshfree method converge faster and are more accurate than those of standard FEM.The studied meshfree method is effective in the analysis of static responses of complex geometry structures.The amount of discretised nodes within the integration domain used in building MLS shape functions should be in the range from 30 to 60 nodes and should not be less than 20 nodes.

Developing Hybrid cell-edge and cell-node Dissipative Compact Scheme for Complex Geometry Flows

Developing high resolution finite difference scheme and enabling the use of this scheme on complex geometry are the aims of this study.High resolution has been achieved by Dissipative Compact Schemes(DCS),however,according to the recent research,applications of DCS on complex geometry may have serious problem for that the Geometric Conservation Law(GCL)is not satisfied,and this may cause numerical instability.To cope with this problem,a new scheme named Hybrid cell-edge and cell-node Dissipative Compact Scheme(HDCS)has been formulated.The formulation of the HDCS contains two steps.First,a new central compact scheme is formulated for the purpose of conveniently fulfilling the GCL,and then dissipation is added on the central scheme by high-order dissipative interpolation of cell-edge variables.The solutions of Euler and Navier-Stokes equations show that the HDCS can be applied successfully on complex geometry,while the DCS may suffer numerical instabilities.Moreover,high resolution of the HDCS may be observed in the test of scattering of acoustic waves by multiple cylinders.

Hybrid Diffuse and Sharp Interface Immersed Boundary Methods for Particulate Flows in the Presence of Complex Boundaries

A coupling framework that leverages the advantages of the diffuse and sharp interface immersed boundary(IB)methods is presented for handling the interaction among particles and particles with the static complex geometries of the environment.In the proposed coupling approach,the curvilinear IB method is employed to represent the static complex geometries,a variant of the direct forcing IB method is proposed for simulating particles,and the discrete element method is employed for particle-particle and particle-wall collisions.The proposed approach is validated using several classical benchmark problems,which include flow around a sphere,sedimentation of a sphere,collision of two sedimenting spheres,and collision between a particle and a flat wall,with the present predictions showing an overall good agreement with the results reported in the literature.The capability of the proposed framework is further demonstrated by simulating the interaction between multiple particles and a wall-mounted cylinder,and the particle-laden turbulent flow over periodic hills.The proposed method provides an efficient way to simulate particle-laden turbulent flows in environments with complex boundaries.

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional （3D） complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.

Multiderivative Combined Dissipative Compact Scheme Satisfying Geometric Conservation Law III: Characteristic-Wise Hybrid Method

Based on newly developed weight-based smoothness detectors and non-linear interpolations designed to capture discontinuities for the multiderivative com-bined dissipative compact scheme(MDCS),hybrid linear and nonlinear interpolations are proposed to form hybrid MDCS.These detectors are derived from the weights used for the nonlinear interpolations and can provide suitable switches between the linear and the nonlinear schemes to realize the characteristics for the hybrid MDCS of capturing discontinuities and maintaining high resolution in the region without large discontinuities.To save computational cost,the nonlinear scheme with characteris-tic decomposition is only applied in the detected discontinuities region by specially designed hybrid strategy.Typical tests show that the hybrid MDCS is capable of cap-turing discontinuities and maintaining high resolution power for the smooth region at the same time.With the satisfaction of the geometric conservative law(GCL),the MDCS is further applied on curvilinear mesh to present its promising capability of handling pragmatic simulations.

Positivity-Preserving Runge-Kutta Discontinuous Galerkin Method on Adaptive Cartesian Grid for Strong Moving Shock

In order to suppress the failure of preserving positivity of density or pres-sure,a positivity-preserving limiter technique coupled with h-adaptive Runge-Kutta discontinuous Galerkin(RKDG)method is developed in this paper.Such a method is implemented to simulate flows with the large Mach number,strong shock/obstacle interactions and shock diffractions.The Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is also presented.This ap-proach directly uses the cell solution polynomial of DG finite element space as the interpolation formula.The method is validated by the well documented test ex-amples involving unsteady compressible flows through complex bodies over a large Mach numbers.The numerical results demonstrate the robustness and the versatility of the proposed approach.