Discretization Research Articles

An Implicit Transformer-based Fusion Method for Hyperspectral and Multispectral Remote Sensing Image

There is an effective way to enhance the spatial resolution of hyperspectral remote sensing images by fusing them with multispectral remote sensing images. However, most of the existing deep fusion techniques adopt discretized explicit models to approximate the complex continuous nonlinear mapping in the fusion process, leading to limitations in enhancing the fidelity of spatial details. Additionally, existing algorithms commonly utilize discrete methods such as bilinear or bicubic interpolation during the hyperspectral upsampling process, leading to the loss of crucial spatial-spectral features. To this end, this study proposes a novel Implicit Transformer Fusion Generative Adversarial Network (ITF-GAN), which incorporates the continuity perception mechanism of implicit neural representation with the powerful self-attention mechanism of the Transformer architecture, which uses point-to-point implicit functions aiming to efficiently process information in both spatial and spectral dimensions. Besides, a guided implicit neural sampling module is introduced in the hyperspectral image up-sampling process to enhance the coordinated expression of features in the spatial and spectral domains, which improves the spatial resolution and spectral fidelity of the fused image during the upsampling process. A series of fusion experiments including 4x, 8x, and 16x scale factors have shown that ITF-GAN has significant advantages over current popular fusion algorithms in both objective evaluation indicators and subjective visual evaluation.

Application of Extended Eddy-viscosity and Elliptic-Relaxation Approaches to Turbulent Convective Flow in a Partially Divided Cavity

The paper reports on the numerical turbulence model in predicting mass, momentum and heat transfer in a partially divided cavity heated from the side using buoyancy-extended eddy-viscosity and elliptic relaxation approach with the algebraic expressions for the Reynold stress tensor and turbulent heat flux vector. The CDS (central differencing scheme) and LUDS (linear upwind differencing scheme) were used as the discretization method and the governing equations were solved using the finite volume method and Navier-Stokes solver. Validation of the model has been carried out by experimental data of convective flow in the cavity as well as by numerical data DNS (direct numerical simulation). The model agrees very well with the experiment and DNS and it is also able to demonstrate the performance which is comparable to that of the previous advanced second-moment closure model (SMC) in the literature. The results show that the model is suitable for use in simulations of the turbulent convective flow in a cavity with partition and it has the potential to be applied to more complex cavities and a wide range of turbulence levels.

A Second-Order Length-Preserving and Unconditionally Energy Stable Rotational Discrete Gradient Method for Oseen-Frank Gradient Flows

A Second-Order Length-Preserving and Unconditionally Energy Stable Rotational Discrete Gradient Method for Oseen-Frank Gradient Flows

Quantum simulation of discrete linear dynamical systems and simple iterative methods in linear algebra

Quantum simulation is capable of simulating certain dynamical systems in continuous time—Schrödinger’s equations being the most direct and well known—more efficiently than classical simulation. Any linear dynamical system can in fact be transformed into a system of Schrödinger’s equations via a method called Schrödingerisation (Jin et al. 2022. Quantum simulation of partial differential equations via Schrödingerisation. ( https://arxiv.org/abs/2212.13969 ) and Jin et al. 2023. Phys. Rev. A 108 , 032603. ( doi:10.1103/PhysRevA.108.032603 )). We show how Schrödingerisation allows quantum simulation to be directly used for the simulation of continuous-time versions of general (explicit) iterative schemes or discrete linear dynamical systems. In particular, we use this new method to solve linear systems of equations and for estimating the maximum eigenvector and eigenvalue of a matrix, respectively. This method is applicable using either discrete-variable quantum systems or on hybrid continuous-variable and discrete-variable quantum systems. This framework provides an interesting alternative to solve linear algebra problems using quantum simulation.

A field inversion and symbolic regression enhanced Spalart–Allmaras model for airfoil stall prediction

A data-driven turbulence modeling method based on symbolic regression (SR) is proposed in this paper to enhance the prediction accuracy of the Spalart–Allmaras (SA) model for airfoil stall. Unlike traditional methods that rely on neural networks and lack physical interpretability, this paper utilizes SR to establish an analytic expression mapping local flow field variables to the SA model correction factor β. The training data are obtained through field inversion with the discrete adjoint method in the flow field of the S809 airfoil. Additionally, a relearning approach proposed in this paper is applied to the SR process to address the issue arising from the multi-solution nature of field inversion. The SA model embedded with β, referred to as the SA-SR model, can be integrated into computational fluid dynamics solvers with negligible computational cost. The generalization performance of the SA-SR model is tested under various conditions and airfoil types. The results indicate that the new model improves the predictive capability for airfoil stall without compromising the performance of the baseline SA model for attached flows.

The utilization of the discrete event simulation method in scheduling repetitive construction

The utilization of the discrete event simulation method in scheduling repetitive construction

Comparison of heat transfer characteristics of a heat exchanger with straight helical tube and a heat exchanger with coiled flow reverser

Vertical shell and coiled tube heat exchangers are an important aspect of numerous industrial processes and are renowned for their high heat transfer efficiency, versatility, compact design, easy manufacturability and maintenance. As a novel technique, a coiled flow reverser was used to enhance the thermal performance of shell and coiled tube heat exchangers. Therefore, the impact of employing a coiled flow reverser on the heat transfer characteristics and pressure drop in a counter-current shell and tube heat exchanger is investigated using both numerical and experimental approaches. The results are compared to a straight helical coiled tube heat exchanger. The numerical simulations were carried out by employing computational fluid dynamics (CFD) techniques, and the SIMPLE scheme was implemented through the second-order upwind discretization method. The laboratory device used for experimentation features a shell with an inner diameter of 15 cm, and the coil pitch and coil diameter of both coiled tubes are set at 3 cm and 12.1 cm, respectively. Results reveal that the overall heat transfer coefficient of the heat exchanger with a coiled flow reverser is 39.5 % to 49.9 % higher than that of the heat exchanger with a straight helical tube. The effectiveness ratio of the heat exchanger with coiled flow reverser to the straight helical tube heat exchanger ranges from 1.254 to 1.379. Optimal results are observed at a cold flow rate of 6 LPM and a hot flow rate of 1 LPM for both heat exchangers. Additionally, the coiled flow reverser enhances heat transfer performance at the cost of increased pressure drop. Numerical results are in good agreement with experimental findings, establishing the reliability of the study.

Reduced-Order Extended Observer and Correction Discretization Method for High-Speed Permanent Magnet Synchronous Motor With LCL Filter

Reduced-Order Extended Observer and Correction Discretization Method for High-Speed Permanent Magnet Synchronous Motor With LCL Filter

Inverse-Free DZNN Models for Solving Time-Dependent Linear System via High-Precision Linear Six-Step Method.

Time-dependent linear system (TDLS) is usually encountered in scientific research, which is the mathematical formulation of many practical applications. Different from conventional inverse-need models, by utilizing zeroing neural network (ZNN) method twice, an inverse-free continuous ZNN (CZNN) model is developed for solving TDLS. For conveniently practical use, a discrete model is naturally desired. Superior to conventional discretization methods, a general linear six-step (LSS) method with the seventh-order precision and five variable parameters is proposed for the first time. Constraints about five variable parameters are theoretically analyzed to guarantee the efficacy of the general LSS method. Within constraints, 12 specific LSS methods are further developed. Aided with the general LSS method, an inverse-free discrete ZNN (DZNN) is proposed and termed DZNN-LSS model, and its precision is greatly improved compared with conventional discrete models. For comparison, three conventional discretization methods are also utilized to generate DZNN models. Detailed theoretical analyses are provided to prove the efficacy of relevant models. In addition, a specific TDLS example is considered to show the effectiveness and superiority of the DZNN-LSS model. More than that, applications to manipulator control and sound source localization are conducted to illustrate the applicability of the DZNN-LSS model.

Dynamics analysis and chatter control of a polishing and milling nonideal flexible manipulator

Some robots are designed to be lightweight and flexible, enabling them to access small and challenging paths in various applications. These features enable robots to collaborate with humans in performing specific production tasks. However,the movement of the flexible manipulator can become self-excited when handling a cutting tool on a workpiece, which can lead to a control problem. This article presents a control solution for lightweight robotic manipulators with rotating tools, such as polishing and milling, using smart actuators. The control discretization method is also introduced to facilitate integration into digital controllers. The paper starts by describing the governing equations of the non-ideal flexible manipulator for polishing and milling and analysing its dynamic behavior. Subsequently, models for the controllers of the DC motor-only actuators and the hybrid (shape memory alloy and DC motors) actuators were formulated using shape memory alloy and a suboptimal control scheme known as the discrete state-dependent Riccati equation. The coupled system was analysed dynamically, and it was observed that it exhibits chaotic behavior. The results suggest that incorporating smart actuators into the motor group of the system could decrease the positioning error of the manipulator and significantly reduce the oscillation of the robot’s end-effector.

Estimation of Thermal Properties of Solid–Liquid Phase Change Material Using Fuzzy Inference Methods

The accurate measurement of thermal properties in phase change materials holds significant importance for engineering applications. This research introduces fuzzy inference methods to estimate the thermal properties of phase change materials. The solution to the coupled heat transfer involving radiation and conduction in material is achieved through a hybrid approach, which combines the finite volume method with the discrete ordinate method. The estimation process is structured as an inverse problem, where the temperature on the material surface acts as the measurement signal for conducting the inverse analysis. Both the fuzzy inference method and the decentralized fuzzy inference method are utilized to address the inverse heat transfer problem. This enables the determination of latent heat and thermal conductivities for both solid and liquid regions within the phase change material. Retrieval results demonstrate that the thermal properties of phase change materials can be accurately estimated using fuzzy inference methods. While both two fuzzy inference methods perform similarly in estimating a single parameter, the fuzzy inference method has limitations in multiparameter estimation tasks. Conversely, the decentralized fuzzy inference method yields accurate results in simultaneous estimation problems. Furthermore, this method proves robust in estimating the thermal properties of phase change materials, even in the presence of noisy data.

Robust digital-twin airspace discretization and trajectory optimization for autonomous unmanned aerial vehicles

The infiltration of heterogenous fleets of autonomous Unmanned Aerial Vehicles (UAVs) in smart cities is leading to the consumerization of city air space which includes infrastructure creation of roads, traffic design, capacity estimation, and trajectory optimization. This study proposes a novel autonomous Advanced Aerial Mobility (AAM) logistical system for high density city centers. First, we propose a real-time 3D geospatial mining framework for LiDAR data to create a dynamically updated digital twin model. This enables the identification of viable airspace volumes in densely populated 3D environments based on the airspace policy/regulations. Second, we propose a robust city airspace dynamic 4D discretization method (Skyroutes) for autonomous UAVs to incorporate the underlying real-time constraints coupled with externalities, legal, and optimal UAV operation based on kinematics. An hourly trip generation model was applied to create 1138 trips in two scenarios comparing the cartesian discretization to our proposed algorithm. The results show that the AAM enables a precise airspace capacity/cost estimation, due to its detailed 3D generation capabilities. The AAM increased the airspace capacity by up to 10%, the generated UAV trajectories are 50% more energy efficient, and significantly safer.

Simulation of coupled elasticity problem with pressure equation: hydroelastic equation

PurposeThe main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.Design/methodology/approach First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.FindingsNumerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.Originality/valueThis is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

Coupled vibration analysis of the spacecraft with the flexible shaft and solar panels assembly

Dynamical characteristics of the spacecraft with flexible components are investigated in this paper. The component consists of a flexible shaft and solar panels, where solar panels are completely fixed on the shaft. A novel power series multiplier polynomial method is proposed to describe the connecting condition between the shaft and solar panels. The deformation and force matching conditions between the one dimensional shaft and two dimensional panels are established, which is the major contribution about the multidimensional combined structure of this study. According to constraining boundary conditions of the shaft and solar panels, the bending-torsion coupled deformation of the component is expressed. The attitude maneuvering of the central platform and the elastic vibration of the component cause the rigid-flexible coupled effect. Then the discrete dynamical equation of the spacecraft is obtained by using the rigid-flexible coupled modal shapes. The correctness of the proposed method is verified by comparing the natural characteristics with that of the finite element model. The relative error of the higher order frequency is not more than 4%. Numerical results show that the degree of the bending angle polynomial should be at least quadratic to ensure the accuracy of the calculation. The proposed power series multiplier polynomial breaks through the barriers that the traditional Lagrange multiplier cannot express the continuous constraint. It is an effective method to fit the constrained boundary by the polynomial rather than by the discrete multipoint method.

Probing the Gamma-Ray Emission Region of Five TeV Flat Spectrum Radio Quasars

The location of γ-ray emission of blazars remains a contested topic, inspiring the development of numerous investigative techniques to address this issue. In this work, we analyzed Fermi γ-ray light curves in the GeV and MeV bands, employing the discrete cross-correlation function method to discern time lags between the two bands. For 4C +21.35, Ton 599, B2 1420+32, and PKS 1510-089, we identified a time lag spanning several days, while for PKS 1441+25, the time lag was not statistically found. The results imply that the soft photons necessary for inverse Compton scattering predominantly originate from the dusty torus in the first four sources, whereas for PKS 1441+25, they seem to be sourced mainly from the broad-line region. Further analysis of the opacity (τ γ γ ) and the GeV spectra study supports the conclusion that the location of the dissipation region must be beyond the BLR to avoid significant absorption. Notably, for PKS 1441+25, the emission region is also posited to lie outside yet proximate to the BLR. The parameters of describing the emission region were obtained by fitting broadband spectral energy distribution with contemporaneous observation data. Our findings suggest that for the five TeV FSRQs, during TeV flaring events, the jet appears to maintain an equilibrium between the energy density of the magnetic field and that of the particles for all investigated sources, with the exceptions of 4C +21.35 and PKS 1441+25. In terms of the overall jet power, particle energy is the dominant contributor, and the observed blazar radiation cannot be solely attributed to the magnetic field, except in the case of 4C +21.35. Consequently, magnetic reconnection is unlikely to be the primary mechanism behind particle acceleration in these systems.

Extended finite element method with cell-based smoothing for modeling frictional contact crack-induced acoustic nonlinearity involving distorted mesh

This study proposes a 2D cell-based smoothed extended finite element method (CS-XFEM) for accurate and efficient simulation of nonlinear ultrasonic wave propagation in solid structures, specifically addressing the effects of frictional contact in cracks. Traditional mesh discretization methods for cracks often suffer from mesh distortion and computational instability owing to their high aspect ratios. To overcome this, CS-XFEM integrates a cell-based smoothing technique into XFEM to model the frictional contact of a crack. A comprehensive numerical example demonstrates the advantages of CS-XFEM. The results show that CS-XFEM exhibits a higher convergence rate and enables a larger critical time increment than XFEM. Specifically, the critical time increment of CS-XFEM was found to be twice that of XFEM, leading to a 50% reduction in the total computational time. These findings confirm that CS-XFEM is an efficient, accurate, and robust numerical method for studying the acoustic nonlinearity induced by crack-induced frictional contact.

Impact of different numerical approaches on the magnetocaloric effect modeling

As an ecological alternative to the conventional refrigeration technology, magnetocaloric refrigeration is still facing scientific and technological challenges hindering their application. Magnetocaloric devices rely on the magnetocaloric effect, where temperature variations result from magnetic field changes. The correct implementation of the magnetocaloric effect in numerical models is crucial before prototyping the related solutions. Here, we present a comparison between the three most used numerical methods to simulate the magnetocaloric effect: continuous temperature change, discrete temperature change step and heat source obtained from adiabatic temperature. By varying the time and space steps, it was observed that the continuous temperature change method is the most appropriate for small time steps, but has the largest computational cost. The discrete method can only be applied to small time steps, but is the fastest method. Finally, the adiabatic temperature change power source method can be applied in the entire range and is the one that presents the best results for larger time steps.

Analysis and Modeling of Queuing System Simulation in Payment Process at Minimarket (Case Study of Minimarket X Yogyakarta)

The purpose of this study is to analyze the queuing system at Minimarket X so as to minimize the number of queues. The method in this writing uses a discrete event simulation method with the help of Anylogic software. Simulation is used as a method to analyze problems that exist in the payment system at the Minimarket, while Anylogic software as an illustration that can facilitate the queuing system so that it is clearer and easier to understand. Observations conducted for 20 days show that the queuing system applied by Minimarket X is Multi Channel – Single Phase and First Come First Served (FCFS) or First In First Out (FIFO) with 4 cashier officers (M/M/4). Each cashier is tasked with serving customers and processing the products purchased by scanning them and packing each grocery. The average number of Minimarket visitors every day is 22 customers per hour with a queue length of 10 to 11 customers per hour, and the average waiting time for customers to be served is around 27 to 30 minutes. The simulation used through anylogic software obtained error results within 67 minutes 43 seconds.

A new kinematic model for revolute clearance joints with noncircular bushing and pin in planar multibody systems

Due to manufacturing error and wear, the profiles of bushing and pin in revolute clearance joints are usually non-circular, while traditional kinematic models for revolute clearance joints generally assume a circular profile for the pin or bushing, which may not be able to accurately capture the kinematic behavior of real revolute joints in multibody systems. To this end, a curvature center method was proposed for the kinematic modeling of revolute clearance joints with noncircular bushing and pin. In the proposed method, both the noncircular bushing and pin were discretized, and the curvature center and radius corresponding to each discrete point were then obtained to detect the point of contact between the bushing and pin based on their kinematic constraints. The number of discrete points was determined based on both accuracy and efficiency of the proposed curvature center method, which was then applied to model a revolute clearance joint between the slider and the link of a slider-crank mechanism. Via comparison to the traditional geometric center method and discrete point method, it was validated that the proposed method can achieve a comparable accuracy in modeling revolute joints with a circular pin. The unique advantage of the proposed curvature center method is its ability in modeling revolute joints with noncircular bushing and noncircular pin, which was also demonstrated via a series of simulations. Simulation results show a significant influence of noncircular bushing and noncircular pin on the dynamic response of the multibody system. This implies the promising value of the proposed method in studying the effect of manufacturing error or wear of both bushing and pin on the system’s dynamic performance and/or service life.

Multiphase Numerical CFD Simulation of the Hydrothermal Liquefaction Process (HTL) of Sewage Sludge in a Tubular Reactor

This article presents multiphase numerical computational fluid dynamics (CFD) for simulating hydrothermal liquefaction of sewage sludge in a continuous plug-flow reactor. The discrete particle method (DPM) was used to analyze the solid particles’ interaction in liquid–solid high shear flows to investigate coupling computational fluid dynamics (CFD). Increasing solid particles’ interactions were observed with the increasing liquid velocity. The study examined the influence of parameters such as flow rate, temperature, and residence time on the efficiency of bio-oil production. An increase in temperature from 500 to 800 K caused an increase in the amount of biocrude oil produced from 12.4 to 32.9% within 60 min. In turn, an increase in the flow rate of the suspension from 10 to 60 mL/min caused a decrease in the amount of biocrude oil produced from 38.9 to 12.9%. This study offers insights into optimizing the flow channel of tubular reactors to enhance the HTL conversion efficiency of sewage sludge into biocrude oil. A parametric study was performed to investigate the effect of the slurry flow rate, temperature, and the external heat transfer coefficient on the biocrude oil production performance. The simulation data will be used in the future to design and scale up a large-scale HTL reactor.