OPTIMIZATION OF FUSED DEPOSITION MODELING PROCESS PARAMETERS FOR MATERIAL CONSUMPTION AND MANUFACTURING TIME

This paper aims to explore the effect of layer thickness, infill density and build orientation on the material consumption and manufacturing time of specimens printed by the fused deposition modeling process. Specimens in accordance with ASTM Standards were printed by varying the process parameters such as layer thickness, infill density and build orientation. Time required to manufacture the part and amount of material consumed during the process are recorded. Increase in infill density results into increase in material consumption and manufacturing time. Layer thickness and build orientation also impacts manufacturing time and material consumption respectively. With increased application of FDM process, determining the process parameter to decrease the material consumption and manufacturing time shall help the FDM practitioners globally. Present work elucidates the optimization of FDM process parameters to achieve minimum material consumption and manufacturing time.

Numerical Solution of the Two-Phase Stefan Problem in the Enthalpy Formulation with Smoothing the Coefficients

To simulate heat transfer processes with phase transitions, the classical enthalpy model of Stefan is used, accompanied by phase transformations of the medium with absorption and release of latent heat of a change in the state of aggregation. The paper introduces a solution to the two-phase Stefan problem for a one-dimensional quasilinear second-order parabolic equation with discontinuous coefficients. A method for smearing the Dirac delta function using the smoothing of discontinuous coefficients by smooth functions is proposed. The method is based on the use of the integral of errors and the Gaussian normal distribution with an automated selection of the value of the interval of their smoothing by the desired function (temperature). The discontinuous coefficients are replaced by bounded smooth temperature functions. For the numerical solution, the finite difference method and the finite element method with an automated selection of the smearing and smoothing parameters for the coefficients at each time layer are used. The results of numerical calculations are compared with the solution of Stefan’s two-phase self-similar problem --- with a mathematical model of the formation of the ice cover of the reservoir. Numerical simulation of the thawing effect of installing additional piles on the existing pile field is carried out. The temperature on the day surface of the base of the structure is set with account for the amplitude of air temperature fluctuations, taken from the data of the Yakutsk meteorological station. The study presents the results of numerical calculations for concrete piles installed in the summer in large-diameter drilled wells using cement-sand mortars with a temperature of 25 °С. The distributions of soil temperature are obtained for different points in time

The modified high-order Haar wavelet scheme with Runge–Kutta method in the generalized Burgers–Fisher equation and the generalized Burgers–Huxley equation

In this work, we focus on the modified high-order Haar wavelet numerical method, which introduces the third-order Runge–Kutta method in the time layer to improve the original numerical format. We apply the above scheme to two types of strong nonlinear solitary wave differential equations named as the generalized Burgers–Fisher equation and the generalized Burgers–Huxley equation. Numerical experiments verify the correctness of the scheme, which improves the speed of convergence while ensuring stability. We also compare the CPU time, and conclude that our scheme has high efficiency. Compared with the traditional wavelets method, the numerical results reflect the superiority of our format.

Development of the LBM non-isothermal flows with arbitrarily large Mach number

При решении задач динамики жидкостей и газов в области малых скоростей потока и при изотермических условиях с успехом применяется метод решеточных уравнений Больцмана (LBM). Для решения дискретного уравнения Больцмана может быть использован новый метод Particles-on-Demand (PonD), в котором в каждой точке сетки дискретизация функции распределения в пространстве скоростей центрирована относительно текущей скорости потока. В отличие от классического LBM, метод PonD применим не только для задач с малыми скоростями потока и при изотермических условиях. В данной работе реализован метод PonD D1Q5 с итерационным расчетом скорости переноса и явным расчетом первых трех моментов, включая скорости переноса. Показано, что рассмотренная модификация метода PonD хоть и накладывает ограничения на параметры, позволяет проводить расчеты в большем диапазоне допустимых скоростей. The purpose of the paper is to demonstrate applicability of the Particle on Demand (PonD) D1Q5 method with the explicit calculation of the first three moments to problem with high speed of the flow. The standard LBM is applicable for small flow velocities. Thus to overcome this limitation we use PonD. In this work, we use conservative version of PonD - the D1Q5 method with the explicit calculation of the first three moments. Methodology. The Pond over LBM was applied to the Riemann problem in order to demonstrate the advantage of the method. In this work, we choose the case when contact discontinuities could propagate at variable speed. Findings. If the interpolation pattern is fixed relative to the point at which there is a current update of the discrete distribution function, then the transfer step can be written explicitly, thus the scheme is conservative. On the other hand, this imposes additional restrictions on the temperature and the flow rate. But even if the PonD scheme is limited to a fixed interpolation pattern, it is possible to simulate flows with larger values of the Mach number than in the case when the classical method of lattice Boltzmann equations is used. Originality/value. In the described particular case of the PonD method, it is possible to avoid iterations by calculating the temperature and velocity values directly at a new time layer. In this work, we have investigated the properties and the range of applicability (admissible values of temperature and velocity) of such modification of PonD.

APPLICATION OF THE IMPLICIT METHOD OF CHARACTERISTICS ON A REGULAR GRID FOR SOLVING THE GAS-DYNAMIC PROBLEMS

In this work, an implicit method is proposed to numerically solve a system of the onedimensional nonstationary equations of gas dynamics transformed by the method of characteristics. Internal points of the channel for a solid-propellant charge are considered at a preignition period of the solid-propellant rocket engine operation. The use of the implicit method makes it possible to calculate the values of gas-dynamic parameters at nodal points of the regular coordinate grid. Calculations of the gas-dynamic parameters both when integrating over time and along the spatial coordinate are performed with the second order of accuracy. Both subsonic and supersonic flows are studied. It is shown that, when predicting the expected pressure value during the transition from one time layer to another with the second order of accuracy, the twenty-fold efficiency of the implicit method is achieved in comparison with the explicit difference method. The trial calculation is performed.

The effectiveness of spectral decomposition-based layer thickness estimation: A seismic physical modeling example

We have constructed a channel complex model at a scale of 1:10,000 by stacking 3D-printed polylactide layers with negative relief meandering channels. This model was subjected to an ultrasonic common-offset acquisition in a water tank (with the water filling the channels), and the result was treated as a zero-offset 3D acoustic reflection seismogram, receiving a deterministic deconvolution and a poststack migration as data treatment. We then developed an algorithm to yield volumes of estimated two-way time layer thickness from multiple-frequency volumes obtained through the short-time Fourier transform. The estimated thicknesses were compared with the measurements of the physical model obtained through X-ray computed tomography. Despite the strong signal attenuation and imaging issues, the results were rather satisfactory, increasing the confidence in using spectral decomposition for quantitative seismic analysis.

TWO-DIMENSIONAL SPLITTING SCHEMES FOR HYPERBOLIC EQUATIONS

The article considers splitting schemes in geometric directions that approximate the initial-boundary value problem for p-dimensional hyperbolic equation by chain of two-dimensional-one-dimensional problems. Two ways of constructing splitting schemes are considered with an operator factorized on the upper layer, algebraically equivalent to the alternating direction scheme, and additive schemes of total approximation. For the first scheme, the restrictions on the shape of the region G at p=3 can be weakened in comparison with schemes of alternating directions, which are a chain of three-point problems on the upper time layer, the region can be a connected union of cylindrical regions with generators parallel to the axis OX3. In the second case, for a three-dimensional equation of hyperbolic type, an additive scheme is constructed, which is a chain «two-dimensional problem – one-dimensional problem» and approximates the original problem in a summary sense (at integer time steps). For the numerical implementation of the constructed schemes – the numerical solution of two-dimensional elliptic problems – one can use fast direct methods based on the Fourier algorithm, cyclic reduction methods for three-point vector equations, combinations of these methods, and other methods. The proposed two-dimensional splitting schemes in a number of cases turn out to be more economical in terms of total time expenditures, including the time for performing computations and exchanges of information between processors, compared to traditional splitting schemes based on the use of three-point difference problems for multiprocessor computing systems, with different structures of connections between processors type «ruler», «matrix», «cube», with universal switching.

Stability analysis of convection–diffusion equations of different finite-element spaces at discrete times

The convection-controlled diffusion problem is hyperbolic in nature and its solutions tend to have numerical shocks. To solve the problem of time dependence, it is common to use finite-element method in the spatial region, while the algorithm proposed in the past with finite differential procedure is mostly limited to the fixed finite-element grid of the spatial region. We often need to use different finite-element spaces at different times, such as the spread of flame, oil and water frontier problems; so many mathematicians and engineers have set their sights on the use of dynamic finite-element space, but also put across a lot of dynamic finite-element methods in giving the general parabola problem of the variable-mesh finite-element method. The principal purpose of this paper is to adopt different spatial grids for different time layers, and project the approximate solution of the previous time layer to the present time layer to act as the initial value of the current layer to enable us deduce the stability at discrete times.

Nonsplit complex-frequency shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 1: Method

The absorbing boundary condition plays an important role in seismic wave modeling. The perfectly matched layer (PML) boundary condition has been established as one of the most effective and prevalent absorbing boundary conditions. Among the existing PML-type conditions, the complex frequency shift (CFS) PML attracts considerable attention because it can handle the evanescent and grazing waves better. For solving the resultant CFS-PML equation in the time domain, one effective technique is to apply convolution operations, which forms the so-called convolutional PML (CPML). We have developed the corresponding CPML conditions with nonconstant grid compression parameter, and used its combination algorithms specifically with the symplectic partitioned Runge-Kutta and the nearly analytic SPRK methods for solving second-order seismic wave equations. This involves evaluating second-order spatial derivatives with respect to the complex stretching coordinates at the noninteger time layer. Meanwhile, two kinds of simplification algorithms are proposed to compute the composite convolutions terms contained therein.