COMSOL Multiphysics Software Research Articles

Design and Development of Mathematical and Thermal Load Modelling for Induction Heating Systems

Complex systems can be modelled and their performance can be anticipated with the help of current computers and software tools that are applicable in real-world situations. The IH system is capable of being modelled and put through analysis in two distinct domains. The modelling is compatible with the research and applications in their entirety, including the particular control circuits and converter systems under consideration. The use of electrical modelling makes it possible to analyse the features of the IH load and to examine the variations in the parameters in order to determine the impact of load variations. In order to gain an understanding of the load's flow and thermal distribution, thermal modelling is a crucial component of this study. For the purpose of analysing the performance of the coil in terms of flux and temperature from its comparable model, the COMSOL Multiphysics software that is utilised is utilised.

Оценка возможности использования пакета COMSOL для моделирования теплогидравлических процессов в реакторах ВВЭР

Thermophysical modeling software packages are widely used for modeling various thermohydraulic processes including modeling nuclear reactor fuel element. Verification of such calculations is possible only if we compare the results obtained in different packages (cross-verification) due to the extreme complexity or impossibility to carry out analytical calculations or experiments. It seems relevant to carry out cross-verification of calculations of WWER reactor fuel elements in COMSOL Multiphysics software package using another popular ANSYS software package. Models of fuel elements for the analysis of thermohydraulic processes are developed using thermophysical modeling packages ANSYS Fluent and COMSOL Multiphysics for cross-verification of calculation results in the latter. Modeling of temperature fields is carried out based on the heat conduction equation with an internal heat source. The authors have considered several formulations of problems: with fuel element cladding and without it and the central hole in the fuel. Models of fuel elements of the WWER-1000 reactor have been designed for three formulations of problems with different designs in ANSYS Fluent and COMSOL Multiphysics software packages. Temperature fields have been determined using the finite element method. The results of the calculations of the fuel core maximum temperature are presented and compared for two software packages. The maximum relative deviation of the temperatures obtained is 0,97 %. A comparison of the temperature calculations results obtained using models developed in two packages shows that the relative deviation of the values corresponds to the error of the measuring instruments. It confirms that COMSOL Multiphysics can be successfully used to simulate processes in WWER reactors along with ANSYS Fluent.

Aluminium Scandium Nitride Thin Film BAW Resonators (FBAR) for 5G Applications

With the growing impact of 5G communication system, Radio-Frequency Micro-Electro-Mechanical-Systems (RF MEMS) technologies are becoming increasingly important in our daily lives. Development of filters in the 5G frequency range is being carried out tremendously. The rapid rise in smartphone usage, heightened data demands, congestion in the frequency spectrum, and advancements in technology have posed numerous challenges for the development of RF filters. The ultimate goal is to simplify the RF Front End Module used in smart phones, thereby achieving better selectivity. In this paper, a high frequency - high Q thin film bulk acoustic wave resonator of piezoelectric material Aluminium scandium nitride (AlScN)based FBAR is designed for 5G Application. Finite element analysis conducted using the COMSOL Multiphysics software assess the characteristics of the constructed FBAR. The outcomes indicate that the device resonates effectively at 5 GHz, exhibiting commendable return loss and a Q-factor of 1490. A high-resonance-frequency FBAR device offers advantages for 5G mobile technology.

Harnessing CFD simulations and CAD design to optimize hydropower efficiency of a propeller turbine in the Maliringan River: acase study

This research focuses on improving hydropower production by designing a turbine for the MaliringanRiver in Kalimantan Selatan. The rotational speed and torque will be the focus of maximizing power generation. This research used Computational Fluid Dynamics (CFD) to understand how the fluid flows and how efficiently the turbine works. The simulations helped us see how the fluid flowed, the pressure differences, and the speed of the water inside the turbine. The software COMSOL Multiphysics imitates how the fluid behaves and interacts in real-life situations. The utilization of SOLIDWORKS played a crucial role in the turbine's design process, facilitating an accurate representation of the turbine's geometry and the subsequent fabrication of a prototype propeller turbine, featuring an outer diameter measuring 0.27 meters and an inner diameter measuring 0.113 meters. The methodology resulted in a power efficiency of 76.45%, showcasing the possibility of significant enhancements in the efficiency of hydropower generation. The broader ramifications of this study emphasize the feasibility of tailor-made turbines for local hydropower initiatives, thereby supporting Indonesia's renewable energy plan by providing sustainable and efficient energy alternatives. This study emphasizes the collaborative utilization of Computer-Aided Design (CAD) and CFD technologies in the progression of turbine technology, thereby establishing a basis for future investigations in hydropower optimization.

Efficiency of multi-armature linear pulse electromechanical power and speed converters

Introduction. High-speed linear pulse electromechanical converters (LPEC) provide acceleration of the executive element in a short active section to high speed with significant displacement, while power-purpose LPECs create powerful power impulses of the executive element on the object of influence with minor movements. One of the areas of improvement of LPEC is the creation of multi-armature structures. Methodology. To analyze the electromechanical characteristics and indicators of LPEC, a mathematical model was used, which takes into account the interconnected electrical, magnetic, mechanical and thermal processes that occur when connected to a pulse energy source with a capacitive energy storage. The main results of the calculations were performed in the COMSOL Multiphysics software environment and confirmed by experimental studies in laboratory conditions. Results. The features of the electromechanical processes of multi-armature LPECs are established and their indicators are determined. With the help of efficiency criteria, which take into account electrical, power, speed and magnetic indicators in a relative form with different options for their evaluation strategy, it was established that multi-armature LPECs for power purposes have increased efficiency, and for high-speed LPECs the use of multi-armature configurations is impractical. The conducted experimental studies confirm the reliability of the calculated results. Originality. It has been established that almost all multi-armature LPECs for power purposes have higher efficiency compared to a converter with one armature, and for high-speed LPECs it is advisable to use traditional LPECs with one armature. Practical value. On the basis of multi-armature LPECs, models of an electromagnetic UAV catapult, a magnetic pulse press for ceramic powder materials, an electromechanical device for dumping ice and snow deposits from a power line wire, a device for destroying information on a solid-state digital SSD drive have been developed and tested. References 20, tables 4, figures 8.

Numerical Multifield Coupling Model of Stress Evolution and Gas Migration: Application of Disaster Prediction and Mining Sustainability Development

At present, coal mining is gradually shifting towards deep areas, and coal mines under deep mining conditions are more prone to coal and gas outburst accidents. In this research, we aim to explain the causes and mechanisms of dynamic disasters, which are caused by the combined action of static load, gas, and dynamic load on tectonic regions in complex stress field environments. Through numerical simulation using COMSOL Multiphysics software, based on the geological conditions of a mine in Jilin Province, it was found that faults lead to abnormal stress in tectonic regions. The combined action of dynamic and static loads results in excessive stress, causing the fragmentation and displacement of the coal body, leading to coal mine disasters, thus disrupting sustainability. Additionally, the coal matrix gas entering fractures raises the gas pressure and leads to the accumulation of methane near earthquake sources. Dynamic loads accelerate gas desorption in coal and increase porosity and permeability, facilitating rapid gas migration. This influx of gas into the roadways exceeds safety limits. Then, based on these findings and on-site conditions, a set of sustainable measures for coal mines has been proposed. This research offers theoretical guidance for enhancing safety, stability, and sustainability in coal mining processes.

A generalized model for tribovoltaic nanogenerator

Converting mechanical energy into direct-current electric power based on the tribovoltaic effect is a typical characteristic of tribovoltaic nanogenerators (TVNGs). Although this newly discovered physics effect has been devoted to numerous research studies recently, a generalized theoretical model is still missing, thus unable to comprehensively elaborate the working principles of TVNG. Unlike previous qualitative explanations restricted to the conventional diffusion-drift theory, a new theoretical model is proposed according to classical semiconductor physics. Using the model, the governing equation of a TVNG is derived for the first time, which provides possibilities for revealing the variations of basic physical variables whether within the device or in an external circuit. The direct-current output is suggested to be the coupling of the tribovoltaic effect and contact electrification; in detail, it directly results from the movement and realignment of quasi-Fermi levels for excess carriers that are contiguous to the junction/contacting interface under non-equilibrium conditions. Moreover, an equivalent circuit model is established, equivalent to a constant current source parallel to a p–n junction diode according to the lumped parameter circuit theory. Notably, a new term, mechano-induced electric field EM, is defined and introduced to describe the impact of triboelectric charges at interfaces. Furthermore, using the COMSOL Multiphysics software, a dynamic simulation model for TVNGs is proposed, allowing the simulation and calculation of various TVNGs with different geometric constructions and charge distributions.

Моделирование и анализ распределений температуры и термомеханических напряжений в многокристальном электронном модуле

A mathematical thermomechanical model of a multi-chip electronic module (MEM) containing three silicon dies of high-power transistors fixed with a conductive adhesive on a copper plate placed on a radiator is considered in the form of a system of equations of thermal conductivity and thermoelasticity with specified boundary conditions. As a result of computational studies of the model in the COMSOL Multiphysics software environment, distributions of temperature and thermomechanical stresses in the MEM structural elements were obtained depending on the heating power, the thickness of the adhesive and the size of the model defect in the contact connection of one of the MEM crystals with a copper plate in the form of voids in the adhesive layer. It is shown that voids in the adhesive layer lead to a sharply non-uniform temperature distribution over the crystal area and thermomechanical shear stresses in the contact layer that exceed the maximum permissible values for the adhesive. It has been established that thermomechanical stresses decrease with increasing thickness of the adhesive layer and increase with increasing size of the defect(voids) in this layer. The simulation results are in good agreement with the results of measuring the thermal impedance of power transistor crystals using the modulation method, which indicates the correctness and adequacy of the developed model.

Three-dimensional simulation of heat and moisture transfer in woven fabric structures

Fabric structure parameters have a significant impact on the comfort of heat and moisture transfer in garments. Previous numerical simulations required extensive mathematical calculations and mostly investigated one- or two-dimensional models of fiber assembly without considering the weave structure, which is a key parameter, that significantly influences the porosity of the woven structure and consequently its heat and moisture management. While the finite element method supports the simulation of the optimal shape and material properties with better visibility, previous finite element models focused on heat transfer and neglected water vapor transfer in fabrics. In this article, the finite element simulation of heat and moisture transfer in woven fabrics is established based on the testing principle of thermal resistance and moisture resistance tester using COMSOL Multiphysics software. In this simulation, three-dimensional parametric geometrical models of the fabric are created using curve interpolation methods by acquiring the control point coordinates of different weaves (plain, 2/2 balanced twill, and 4/1 unbalanced twill weaves). Heat and moisture transfer properties of fabric models in the horizontal and vertical directions were analyzed, including the heat flux, moisture resistance, water vapor permeability, and water vapor concentration. The article also deals with the effects of weave structure and fabric cover in a range of 72.1–85.1% on the fabric heat flux and water vapor concentration. Comparison between model and experimental results revealed that the three-dimensional simulation can accurately predict the impact of weave pattern and fabric cover on the fabric heat and moisture transfer performance. In addition, this model can be utilized to study the distribution of heat and water vapor within fabrics, providing a theoretical foundation for optimizing heat and moisture comfort in woven fabrics.

Analysis of the influence of local electrodeposition conditions on the accuracy of electrochemical 3D-printing

Problem. Additive manufacturing of metal parts by local electrodeposition is a new promising area with several unique features. The use of electrolysis makes it possible to manufacture metal parts at room temperature with high accuracy and low energy consumption. The accuracy and speed of electrochemical 3D printing are inversely related, and therefore it is important to find optimal electrodeposition conditions at maximum speed without degrading accuracy. The aim of the study. The purpose of this work was to analyse the influence of electrodeposition parameters (electrode distance, electrical conductivity and electrolyte polarization) on the accuracy of metal deposit formation in a computer model and to conduct experimental verification of optimal conditions for electroforming a real object from copper sulphate electrolyte. Methodology of implementation. To achieve the goals of the study, computer simulation of the electrochemical deposition process in the COMSOL Multiphysics software and electrochemical measurements of the characteristics of copper sulphate electrolytes were used. Research results. Computer simulation established optimal conditions: the distance between the edge of the capillary and the surface on which deposition occurs is not more than 0.5 mm, the electrolyte in which the inverse slope of the cathode polarization curve is not lower than 2000 mA/(V×cm2) and the electrical conductivity is not higher than 0.02 S/сm. The influence of the electrolyte composition on the inverse polarization of the cathode process and its electrical conductivity was investigated. The optimal composition of the electrolyte was selected, containing 200 g/L CuSO4, 60 g/L H2SO4, 0.2 g/L KCl and RUBIN T-200 additive. In the selected electrolyte, the value of the inverse slope of the cathode polarization curve is 2120 mA/(V×cm2). Verification of the process of local electrodeposition in the selected electrolyte during the electroformation of a cylindrical object with a diameter of 4 mm and a height of 100 μm was carried out. It was determined that no more than 5 % of metal was deposited outside the capillary of the working electrode. Conclusions. The results of the work can be used to create an electrochemical 3D printing system. Further research should be aimed at approbation of the obtained parameters of local electrodeposition in an installation that simulates the operation of a 3D printer and establish the accuracy and speed of printing of a three-dimensional object. Key words: additive manufacturing, electrical conductivity, electrolyte composition, polarization, local electrodeposition.

One-click startup operation of a micro-tubular solid oxide fuel cells using propane/air/steam

The design and performance of a tubular solid oxide fuel cell (SOFC) integrated with a one-click startup thermal self-sustaining reformer is presented, allowing the tubular SOFC to operate without any external heat source by using propane/air/steam. Double conductive Ni-pads (DCNPs) are utilized as anode current collector in μT-SOFC. The reformer consists of a combustion chamber for initial heating and a reforming chamber for syngas production for the SOFC. The thermal energy generated by complete propane combustion is adequate to achieve thermally self-sustaining operation and efficiently raise the cell temperature above 700 °C. Operating at 700 °C with pure H2, a cell with a length of 12 cm achieves a maximum power density of 138 mW cm−2, which decreases by 13 %–120 mW cm−2 when using syngas produced by steam-propane reforming. The syngas composition from steam-propane reforming reaches 7.6 % CO and 61.7 % H2. COMSOL Multiphysics software is used to simulate the multi-physical field of a single cell exposed to various syngas compositions, showing that syngas produced by steam-propane reforming exhibits excellent electrochemical performance in tubular SOFCs compared to catalyst partial oxidation (CPOX).

A smart structural optimization method of magnetorheological damper for ultra-precision machine tool

To address the problem of multi-source vibration in ultra-precision machine tools, a vibration reduction stand was designed by replacing passive damping components with magnetorheological dampers (MRDs). In this work, the structural parameters of MRDs were optimized using an improved pelican optimization algorithm (IPOA) to realize the maximum capability in reducing vibration. Firstly, the working principle of MRDs was explained, and the mathematical models of MRDs were established. Then, an IPOA based on singer chaotic mapping, nonlinear inertia weight factor, and Cauchy mutation strategy was proposed to enhance the global search capability and convergence efficiency of the algorithm. Subsequently, the IPOA was applied to optimize key structural parameters of MRDs, including output damping force, controllable damping range, response time, and power consumption. Finally, COMSOL Multiphysics software was used to verify the effectiveness of the proposed algorithm by comparing the magnetic induction intensity distribution of MRDs before and after optimization, as well as the variation of the four performance indexes under the different applied currents. After being optimized using the proposed IPOA, the MRDs can deliver a larger maximum damping force and a wider damping controllable range, with less power consumption and quick response, which could meet the requirement for vibration suppression of ultra-precision machine tools.

Computational and experimental microfluidics: Total analysis system for mixing, sorting, and concentrating particles and cells.

Body fluids can potentially indicate the presence of non-small cancer cells. Studying these fluids is an emerging field that could be crucial for cancer detection and monitoring treatment effectiveness. Meanwhile, the examination of fluids on a microscopic level is part of the field of microfluidics. This study focuses on the development of a total analysis system that consists of various interconnected structures that are designed to mix, classify, concentrate, and isolate particles in fluids that mimic the behavior of cancer and normal cells. Using the COMSOL Multiphysics software, the device's performance was optimized to use a pressure input of 35 kPa for water or serum and 29.4 kPa for a mixture of liquid and serum samples, which are the optimal pressure inputs. The numerical models were validated by experiments using two types of polystyrene particles, with diameters of 5 and 20 μm. Moreover, the developed system was applied to monitor the behavior of red blood cells. The microfluidic chip is capable of addressing several challenges through visual detections, including mixing tests of two fluids with similar densities, proper particle size classification using Dean flow fractionation, and single-step recovery of large, labeled particles. Finally, the collected particles were examined using an environmental scanning electron microscope to determine their size, and the results demonstrated that successful size separation was achieved, with particles around 20 μm completely separated from the smaller ones.

Trajectory analysis of circulating tumor cells through contorted deterministic lateral displacement array for unruptured trapping: a simulation study

During the multistage process of cancer spread, fragments of cells branch out from the primitive tumor tissue and scatter throughout the body by blood flow stream or lymph, which are termed as circulating tumor cells (CTCs), regarded as an important biomarker for early cancer diagnosis. Deterministic lateral displacement (DLD) is one of the most efficient passive type particle isolation techniques which functions on asymmetrical laminar flow diversion around the placed microposts forming an array, based on the critical dimeter of the device. In this simulation study, complete trajectory analysis of CTCs and WBCs has been performed using COMSOL Multiphysics software with respect to the angular orientation of sample inlet channel. A simulated model of an asymmetric DLD array-based microfluidic device with latest empirical expression for critical diameter has been presented with optimized sample inlet orientation so as to maintain perfect separation efficiency of CTCs along with arrangement for unruptured trapping of CTCs and WBCs. The approximate CTC isolation efficiency was obtained above 95% for sample inlet inclinations up to 0.5°, and further tilt showed reduction in efficiency. Migration angle of CTC trajectory for the modeled device was obtained around 15° for both the sample inlet orientations leading to least intermingling of CTCs and WBCs. A smaller number of CTC–WBC interaction points was observed for clockwise sample inlet tilt compared to the anticlockwise displacement reducing possibilities of cell rupture. The shear rate around central microposts was found to be higher as well in clockwise sample inlet orientation. This simulation study establishes that sample inlet clockwise tilting of maximum around 0.5°–0.7° leads to better CTC isolation and provides better chances of unruptured trapping of cells at the designated outlets. The results of this study provide an approach toward further optimization of DLD devices' functioning and, thus, could help fabricate better DLD-based microfluidic devices for efficient trapping of CTCs.

Deep learning-assisted dual-modal tomography for phase flow rate estimation in two-phase oil-water flow systems

Accurately estimating phase flow rates in multiphase systems is crucial for many industries, where precise measurements are essential for operational efficiency and safety. Addressing this issue, this paper introduces an approach that employs deep learning-assisted dual-modal electromagnetic flow tomography (EMFT) and electrical tomography (ET) to predict both oil and water flow rates in two-phase oil-water flows. To facilitate the generation of the data, we first simulate diverse flow conditions using COMSOL Multiphysics software and the convection–diffusion equation, aiming to create a realistic representation of two-phase oil-water flows. The dual-modal system measurement data, generated from these simulations and simulated by using a dense finite element mesh, provide reliable inputs for the deep learning model. Moreover, this study also integrates experimental data into both the training and testing phases, improving the ability of the proposed approach to estimate flow rates accurately in practical investigations. The results from laboratory experiments demonstrate the potential of the deep learning-assisted dual-modal ET and EMFT approach in effectively resolving the challenges of estimating flow rates in two-phase oil-water flow systems. By combining the deep learning capabilities with dual-modal tomography, this study offers valuable insights for future applications and represents a significant step forward in the field of multiphase flow rate estimation.

A simulation study on the contact‐separation triboelectric nano‐generator for magnetic energy harvester

AbstractTransmission lines and outdoor substations are replete with various forms of micro energy such as wind energy, solar energy, and electromagnetic energy. There exists micro energy in the forms of mechanical vibrations, temperature differences, and humidity in power distribution equipment. Innovative sensor or monitoring methods are needed to maintain the stability and digitisation of the grid. Unfortunately, there is limited study on the power supply of these sensor systems. Triboelectric nanogenerators, which are environmentally friendly and use simple materials, show excellent performance in environmental nano‐energy collection and self‐powered online monitoring. Therefore, environmental energy collection systems based on triboelectric nanogenerators are one of the selected methods to convert magnetic energy in the magnetic field into electrical energy. A model structure was designed using the contact separation mode, which is one of the four working modes of triboelectric nanogenerators, based on the strength of the magnetic field in the environment. This structure mainly consists of the friction layer, electrodes for current conduction, and connected loads. The research includes a comparison of four inherent electrical outputs of the triboelectric nanogenerator: open‐circuit voltage, short‐circuit current, capacitance, and power. COMSOL Multiphysics software was used for all modelling and simulation of the TENG. This software was used for the design, material selection, and static study of the TENG. When the relative dielectric constant was fixed, the output voltage reached and the energy reached . Overall, ideal reference can be provided for researchers studying power supply issues for sensors in complex magnetic field situations and help them design high‐performance TENGs.

Comparison of Single-Phase Mathematical Models for Solid-State Packed Beds for Thermal Energy Storage

This article presents an analytical solution for the evaluation of the thermal performance of packed bed sensible heat storage. The numerical model developed was tested for four different solid storage mediums. The thermal energy equation is solved numerically by deploying the finite difference method. The presented analytical solution is based on a novel mathematical approach. The numerical model was validated using the computer simulation package Comsol Multiphysics v5.3. Our numerical model results are in good agreement with the published experimental data, with an overall difference of ~10%. Hence, the numerical model is an efficient way of evaluating the thermal performance of packed bed thermal energy storage systems compared to other numerical strategies or computer simulation techniques. This proves that the novel analytical model has shown to be a reliable and broadly accurate approach to acquire the thermal performance of sensible heat storage.

Influence of Magnetic Field Configuration in Helicon-Negative Ion Source on Increase of the Efficiency of NNBI Heating System

The influence of the magnetic field configuration on the performance of a helicon-based negative ion source is investigated with simulation experiments. Using COMSOL Multiphysics software, a three-dimensional simulation model for a negative ion source, based on a helicon plasma source, is presented in two magnetic field configurations: uniform and nonuniform configurations. The helicon plasma source employed a Nagoya-type antenna to apply radio-frequency (RF) power at a frequency of 13.56 MHz. The injected gas is hydrogen with a flow of 10 standard cubic centimeters per minute. Using a three-dimensional model, helicon wave propagation in the presence of a magnetic filter and the energy absorption mechanism in the helicon system are investigated. In this context, in the presence of the two magnetic field configurations, the influence of the important parameters’ working pressure and RF power on the optimization of negative ion production under volume mode is studied. Six electromagnetic coils at the same current are used for producing the magnetic field in both cases of uniform and nonuniform configurations. The variation of the electron density and electron temperature, in both regions of driver and expansion, are calculated and represented with respect to the different power and the gas pressure. The simulation results of the negative ion density in the expansion region for the uniform and nonuniform magnetic field configurations are compared. The results indicate that at the same applied current of coils, the negative ion density in the presence of the nonuniform magnetic field is about 1.75 times higher than the negative ion density of the uniform case. Moreover, the results show that the negative ion density is decreased by decreasing the magnetic field of the driver region in the nonuniform cases.

Experimental and theoretical study on the influence of bundle number on corona loss characteristics of ultra high voltage corona cage positive/negative DC conductors at high altitude

The corona loss generated by the conductor is one of the important factors in the design of ultra-high voltage (UHV) direct current transmission lines. This article is based on the high-altitude (2200 m) ultra-high voltage corona cage experimental platform built in Xining City, Qinghai Province. For the LGJ400-35 conductor used in the typical ± 800 kV DC transmission lines from Qinghai to Henan Province, positive and negative polarity corona experiments were conducted on LGJ400-35 conductors with 4, 6, 8, and 10 bundle under surface field strength conditions of 14–22 kV/cm. For the first time, corona loss data of bundle conductors in positive/negative polarity ultra-high voltage DC transmission lines were obtained, Fit the prediction formula for corona loss of transmission lines under different bundle numbers and surface electric field strengths of wires. And build a model in COMSOL Multiphysics software for simulation analysis, fitting the correction factor for the impact of wind speed on wire corona loss. The research results in the article can provide data reference and support for the construction of ultra-high voltage direct current transmission lines in high-altitude areas.

Phase Separation Study on the Storage of Technically Specification Natural Rubber Modified Bitumen

Overloading and climate change are often problems in pavement structures. For this reason, hard asphalt binders have high softening points, are elastic, and have good adhesion, which is needed to improve pavement performance. Asphalt binder performance can be enhanced by adding additives such as natural rubber or natural-rubber-modified asphalt. However, natural-rubber-modified asphalt shows poor storage stability problems. This is due to differences in density and viscosity between the constituent components of natural-rubber-modified asphalt. This study examines the phase separation mechanism in technically specified natural rubber (TSNR) modified asphalt. Prediction of the optimum storage length of modified asphalt before phase separation occurs, using a combined incompressible Navier–Stokes and phase field model and carried out with COMSOL Multiphysics software version 5.5. Experimental validation was conducted at TSNR levels of 8, 10, and 12% at 160 °C for 48 h, with and without sulfur. The simulation showed that the asphalt modified with TSNR experienced phase separation after 12 h of storage at 160 °C under conditions without stirring. This aligns with the experimental results, which showed phase separation at 160 °C after 48 h. Adding sulfur additives did not have much effect on improving storage stability. The combined incompressible Navier–Stokes and phase field model accurately describes the phase separation in TSNR-modified asphalt. The results of this research recommend that the industry store natural-rubber-modified asphalt in a constantly stirred condition to prevent phase separation of modified asphalt. In addition, the results of this research help the industry predict or increase the homogeneity of polymer-modified asphalt production and save time and costs.