scholarly journals A Monotonic Method of Split Particles

2021 ◽  
Author(s):  
Yanilkin Yury ◽  
Shmelev Vladimir ◽  
Kolobyanin Vadim

The problem of correct calculation of the motion of a multicomponent (multimaterial) medium is the most serious problem for Lagrangian–Eulerian and Eulerian techniques, especially in multicomponent cells in the vicinity of interfaces. There are two main approaches to solving the advection equation for a multicomponent medium. The first approach is based on the identification of interfaces and determining their position at each time step by the concentration field. In this case, the interface can be explicitly distinguished or reconstructed by the concentration field. The latter algorithm is the basis of widely used methods such as VOF. The second approach involves the use of the particle or marker method. In this case, the material fluxes of substances are determined by the particles with which certain masses of substances bind. Both approaches have their own advantages and drawbacks. The advantages of the particle method consist in the Lagrangian representation of particles and the possibility of” drawbacks. The main disadvantage of the particle method is the strong non-monotonicity of the solution caused by the discrete transfer of mass and mass-related quantities from cell to cell. This paper describes a particle method that is free of this drawback. Monotonization of the particle method is performed by spliting the particles so that the volume of matter flowing out of the cell corresponds to the volume calculated according to standard schemes of Lagrangian–Eulerian and Eulerian methods. In order not to generate an infinite chain of spliting, further split particles are re-united when certain conditions are met. The method is developed for modeling 2D and 3D gas-dynamic flows with accompanying processes, in which it is necessary to preserve the history of the process at Lagrangian points.

Author(s):  
Ethan Corle ◽  
Matthew Floros ◽  
Sven Schmitz

The methods of using the viscous vortex particle method, dynamic inflow, and uniform inflow to conduct whirl-flutter stability analysis are evaluated on a four-bladed, soft-inplane tiltrotor model using the Rotorcraft Comprehensive Analysis System. For the first time, coupled transient simulations between comprehensive analysis and a vortex particle method inflow model are used to predict whirl-flutter stability. Resolution studies are performed for both spatial and temporal resolution in the transient solution. Stability in transient analysis is noted to be influenced by both. As the particle resolution is refined, a reduction in simulation time-step size must also be performed. An azimuthal time step size of 0.3 deg is used to consider a range of particle resolutions to understand the influence on whirl-flutter stability predictions. Comparisons are made between uniform inflow, dynamic inflow, and the vortex particle method with respect to prediction capabilities when compared to wing beam-bending frequency and damping experimental data. Challenges in assessing the most accurate inflow model are noted due to uncertainty in experimental data; however, a consistent trend of increasing damping with additional levels of fidelity in the inflow model is observed. Excellent correlation is observed between the dynamic inflow predictions and the vortex particle method predictions in which the wing is not part of the inflow model, indicating that the dynamic inflow model is adequate for capturing damping due to the induced velocity on the rotor disk. Additional damping is noted in the full vortex particle method model, with the wing included, which is attributed to either an interactional aerodynamic effect between the rotor and the wing or a more accurate representation of the unsteady loading on the wing due to induced velocities.


2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.


2007 ◽  
Vol 534-536 ◽  
pp. 573-576
Author(s):  
Eugene Olevsky

The directions of further developments in the modeling of sintering are pointed out, including multi-scale modeling of sintering, on-line sintering damage criteria, particle agglomeration, sintering with phase transformations. A true multi-scale approach is applied for the development of a new meso-macro methodology for modeling of sintering. The developed macroscopic level computational framework envelopes the mesoscopic simulators. No closed forms of constitutive relationships are assumed for the parameters of the material. When a time-step of the calculations is finished for one macroscopic element, the mesostructures of the next element are restored from the initial state according to the history of loading. The model framework is able to predict the final dimensions of the sintered specimen on a global scale and identify the granular structure in any localized area for prediction of the material properties.


2018 ◽  
Vol 2018 (1) ◽  
pp. 000409-000414
Author(s):  
David Bernard

Abstract As advanced packaging continues to develop to support novel and emerging technologies, the need for, ideally non-destructive, test and inspection continues to be vital to ensure the quality and assurance of functionality, wherever the package may go. This is made ever more difficult as the package complexity increases, whilst the feature sizes within continue to decrease. X-ray technology has long been an important part of the non-destructive inspection protocol over the history of advanced packaging and will continue to need to play a more important part in the future. This paper will review the advances made in both 2D and 3D X-ray inspection over recent years and the new opportunities that are now starting to be available, especially in 3D, or CT, inspection, that will enable this 120-year-old technology to remain relevant to and supportive of the needs of advanced packaging. To highlight the above, a case study will be presented on the faults that 2D and CT X-ray analysis can find in LEDs during their manufacture. LEDs are a good example of the remarkable developments in packaging and technology over the last 20 years, where the use of higher powers, smaller sized features and increased reliability requirements intensify the need for higher quality, more consistent production output. Flaws cannot be accepted, especially as higher usage powers mean higher operating temperatures which, in turn, then requires very good thermal conductivity in the package to move heat away from key areas. Without good heat dissipation then heat stresses at the interfaces can cause delamination or die fractures, so reducing LED lifetimes. The presence of voids, particularly at the die to package interface, creates air gaps that reduces heat transfer efficiency. As many LEDs are potted, or encapsulated, the only non-destructive test option to check for voiding and other faults is by using 2D and CT X-ray analysis.


Author(s):  
Masahiro Kondo

AbstractIn general, mechanical energy monotonically decreases in a physically consistent system, constructed with conservative force and dissipative force. This feature is important in designing a particle method, which is a discrete system approximating continuum fluid with particles. When the discretized system can be fit into a framework of analytical mechanics, it will be a physically consistent system which prevents instability like particle scattering along with unphysical mechanical energy increase. This is the case also in incompressible particle methods. However, most incompressible particle methods do not satisfy the physical consistency, and they need empirical relaxations to suppress the system instability due to the unphysical energy behavior. In this study, a new incompressible particle method with the physical consistency, moving particle full-implicit (MPFI) method, is developed, where the discretized interaction forces are related to an analytical mechanical framework for the systems with dissipation. Moreover, a new pressure evaluation technique based on the virial theorem is proposed for the system. Using the MPFI method, static pressure, droplet extension, standing wave and dam break calculations were conducted. The capability to predict pressure and motion of incompressible free surface flow was presented, and energy dissipation property depending on the particle size and time step width was studied through the calculations.


Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. S65-S77 ◽  
Author(s):  
Hu Jin ◽  
George A. McMechan ◽  
Bao Nguyen

We have developed a new method of extracting angle-domain common-image gathers (ADCIGs) from prestack reverse time migration (RTM) that has minimal intermediate storage requirements. To include multipathing, we applied an imaging condition for prestack RTM that uses multiple excitation image times. Instead of saving the full-source snapshots at all time steps, we picked and saved only a few of the highest amplitude arrivals, and their corresponding excitation times, of the source wavefield at each grid point, and we crosscorrelated with the receiver wavefield. When extracting the ADCIGs from RTM, we calculated the source propagation direction from the gradient of the excitation times. The result was that the source time snapshots do not have to be saved or reconstructed during RTM or while extracting ADCIGs. We calculated the receiver propagation direction from Poynting vectors during the receiver extrapolation at each time step and the reflector normal direction by the phase-gradient method. With a new strategy that uses three direction vectors (the source and receiver propagation directions as well as the reflector normal direction), we provided more reliable ADCIGs that are free of low-wavenumber artifacts than any two of them do separately, when the migration velocity model was near to the correct velocity model. The 2D and 3D synthetic tests demonstrated the successful application of the new algorithms with acceptable accuracy, improved storage efficiency, and without an input/output bottleneck.


2014 ◽  
Vol 69 (3-4) ◽  
pp. 145-154 ◽  
Author(s):  
Ke-Pei Men ◽  
Kai Zhao

According to the statistical data, a total of 23 M ≥ 8 earthquakes occurred in Mainland China from 1303 to 2012. The seismic activity of M ≥ 8 earthquakes has showed an obvious self-organized orderliness. It should be remarked especially that there were three ordered pairs of M ≥8 earthquakes occurred in West China during 1902 - 2001, of which the time interval in each pair of two earthquakes was four years. This is a unique and rare earthquake example in earthquake history of China and the world. In the guidance of the information forecasting theory of Wen-Bo Weng, based on previous research results, combining ordered analysis with complex network technology, this paper focuses on the summary of the ordered network structure of M ≥ 8 earthquakes, supplements new information, constructs and further optimizes the 2D- and 3D-ordered network structure of M ≥ 8 earthquakes to make prediction research. At last, a new prediction opinion is presented that the future ordered pair of great earthquakes will probably occur around 2022 and 2026 in Mainland China.


Sign in / Sign up

Export Citation Format

Share Document