Two-dimensional Model Research Articles

Orientation-position coupling in the two-dimensional XY model

Orientation-position coupling in the two-dimensional XY model

Simulation and Experimental Investigation of the Effect of Pore Shape on Heat Transfer Behavior of Phase Change Materials in Porous Metal Structures

With the gradual increase in energy demand in global industrialization, the energy crisis has become an urgent problem. Due to high heat storage density, small volume change, and nearly constant transition temperature, phase change materials (PCMs) provide a promising method to store thermal energy. In this work, we designed and fabricated three kinds of porous metal structures with hexagonal, rectangular, and circular pores and explored the phase change process of PCMs within them. A two-dimensional numerical model was established to investigate the heat transfer process of PCMs within different shapes of porous metal structures and analyze the influence of heat source location on the thermal performance of the thermal storage units. Visualization experiments were also carried out to reveal the melting process of PCMs within different porous metal structures by a digital camera. The results show that paraffin in a porous metal structure with hexagonal pores has the fastest melting rate, while that in a porous metal structure with circular pores has the slowest melting rate. Under the bottom heating mode, the melting time of the paraffin in porous metal structures with hexagonal pores is shortened by 18.6% compared to that in porous metal structures with circular pores. Under the left heating mode, the corresponding melting time is shortened by 16.7%. These findings in this work will offer an effective method to design and optimize the structure of porous metal and improve the thermal properties of PCMs.

Selective anticancer effects of plasma-activated saline in 3D tumor model co-culturing of normal cells and cancer cells

Compared with conventional two-dimensional (2D) cell culture model, the 3D tumor model constructed in vitro is better representative of the tumor microenvironment in vivo. Here, we proposed the utilization of 3D tumor model of co-cultured cancer cells and normal cells to evaluate the selective anticancer effects of cold atmospheric plasma-activated saline (PAS), and expected to provide more precise information about PAS-tumor interactions. By cell sorting, we clarified that A375 melanoma cells and HaCaT normal skin cells purified from the 3D multicellular tumor model differ in sensitivity and responsiveness to PAS compared to the 2D culture model. And during the optimization of PAS treatment parameters, we further found that A375 cells were almost completely killed while HaCaT cells were still present in large numbers after 5 d of certain PAS treatment. Our experiment innovatively carries out the selective study of plasma technology in 3D co-culture system and provides a theoretical basis for further clinical and practical applications of PAS.

Brush Seal Performance with Ideal Gas Working Fluid under Static Rotor Condition

The study investigated how variations in pressure ratio affect the leakage flow of a brush seal for both contact and clearance structures, in which the clearance is measured as the distance between the bristles tip and the rotor surface. This investigation utilized the Reynolds-Averaged-Navier–Stokes (RANS) equations alongside a two-dimensional axisymmetric anisotropic porous medium model. To verify the model’s accuracy and dependability, the obtained results were compared with previous numerical results and experimental observations, showing a satisfactory level of agreement. The results indicate that the predominant pressure drop occurs downstream of the bristle pack with the clearance model exhibiting a higher leakage rate compared to the contact model. Leakage increases proportionally with the pressure ratio, while axial velocity gradually rises and radial velocity experiences a significant increase. In conclusion, the leakage in the brush seal contact structure is significantly lower than in the clearance structure, resulting in the best performance.

Applicability of Measurement-based Quantum Computation towards Physically-driven Variational Quantum Eigensolver

Abstract Variational quantum algorithms are considered one of the most promising methods for obtaining near-term quantum advantages; however, most of these algorithms are only expressed in the conventional quantum circuit scheme. The roadblock to developing quantum algorithms with the measurement-based quantum computation (MBQC) scheme is resource cost. Recently, we discovered that the realization of multi-qubit rotation operations only requires a constant number of single-qubit measurements with the MBQC scheme, providing a potential advantage in terms of resource cost. The structure of the Hamiltonian variational ansatz (HVA) aligns well with this characteristic. Thus, we propose an efficient measurement-based quantum algorithm for quantum many-body system simulation tasks, called measurement-based Hamiltonian variational ansatz (MBHVA). We then demonstrate its effectiveness, efficiency, and advantages with the two-dimensional Heisenberg model and the Fermi-Hubbard chain. Numerical experiments show that MBHVA can have similar performance as circuit-based ansatz, and is expected to reduce operation counts during execution compared to quantum circuits, bringing the advantage of running time. We conclude that the MBQC scheme is potentially feasible for achieving near-term quantum advantages in the noisy intermediate-scale quantum (NISQ) era, especially in the presence of large multi-qubit rotation operations.

Quasinormal modes of analog rotating black holes in a two-dimensional photon-fluid model

Quasinormal modes of analog rotating black holes in a two-dimensional photon-fluid model

Investigation on the explosion of ammonia/hydrogen/air in a closed duct by experiments and numerical simulations

In this paper, the evolution of an explosion of hydrogen/ammonia/air flames is studied by experiments and numerical simulations. A two-dimensional (2D) model is utilized, employing the detached eddy simulation (DES) method along with the thickened flame model. Stoichiometric hydrogen/ammonia/air with ∂ = 25% and 50% (ammonia blending ratio in fuel) is considered. The results indicate that the numerical simulations accurately replicate the experimental findings. Both the experimental observations and numerical simulations reveal the formation process of the “tulip” flame. With increasing ∂, both the flame front speed and overpressure decrease. A change in ∂ can alter the contribution of intermediate reactions, thereby influencing the explosion behavior.

Metaheuristic algorithms for calibration of two-dimensional wildfire spread prediction model

Wildfires are complex phenomena with harmful consequences, ranging from environmental and property destruction to loss of human lives. In this sense, predicting wildfire behaviour is essential to mitigate its impacts and consequences. The Rothermel model is the most used fire rate of spread prediction model. However, input parameter uncertainty is a significant source of prediction error. In this paper, we propose the calibration of the input parameters of the fire propagation model by metaheuristic algorithms under a two-stage framework. The fire spread model consists on the Rothermel model in a two-dimensional approach for surface fires. The proposed calibration is performed in two stages iteratively repeated over time: (i) the calibration of the fire spread model’s input parameters and (ii) the wildfire spread prediction using the calibrated input parameters. The calibration was performed by the genetic algorithm, differential evolution, and simulated annealing, which calibrates the surface-area-to-volume ratio, fuel bed depth, live fuel moisture and dead fuel moisture. The symmetric difference between the real and predicted fire map shapes was defined as the fitness function of all three metaheuristic algorithms. For validation, simulations were done on two prescribed fires. The results for the real and estimated fire behaviour were then compared and revealed that all the tested metaheuristic algorithms produce a better fit to the real fire’s perimeter when compared to the uncalibrated Rothermel model. From the results, differential evolution provided the majority of best results when compared to genetic algorithm and simulated annealing algorithms in each scenario.

Analytical model for the mitigation of VOC vapor with horizontal permeable reactive barrier in the contaminated site considering non-uniform source

Volatile organic compounds (VOCs) contamination at the groundwater may cause vapor intrusion and pose significant threats to human health. As a novel low-carbon mitigation technology, a horizontal permeable reactive barrier (HPRB) is proposed to remove the VOC vapor in the vadose zone and mitigate the vapor intrusion risk. To estimate the performance of HPRB in the contaminated site with a non-uniform source, a transient two-dimensional analytical model is developed in this study to simulate the VOC vapor migration and oxidation processes in the layered soil. The analytical model is verified against the experimental results and numerical simulation first and the parameter study is then conducted. The HPRB has good performance for the contaminated sites involving factors including deep source and local soil with low effective diffusivity. To consider the vertical heterogeneity of the local soil, the traditional equivalent homogeneity method has limitations in considering the horizontal migration of VOC vapor and is not suitable for the two-dimensional model. On the contrary, the artificial layered method based on the proposed analytical model has better accuracy and is recommended to be adopted in practice. Leading to the exponential decrease in the VOC vapor concentration at the ground surface, increasing the thickness of HPRB is an effective measure to enhance the performance of HPRB. The fitting exponential function can be applied to determine the minimum design value of the thickness of HPRB in practice.

Exploring historical changes in mountain river hydrodynamics induced by human impact

During the 20th-century many mountain rivers in Europe were subjected to intensive human impacts which substantially modified their channel morphology. How these changes affected river hydrodynamics and response to floods remains uncertain. In this work, we perform hydraulic modelling using data from archival aerial photos to explore relations between hydraulic parameters of floods and human-induced channel incision occurring on the Czarny Dunajec River (Polish Carpathians) between 1964 and 2012. Data on vertical position of the channel used for two-dimensional modelling of flood flows were extracted (as Digital Elevation Models DEMs) from archival aerial photos from 1964 and 1983 and ALS (Airborne Laser Skanning)-derived DEM from 2012. Water depth, flow velocity, bed shear stress, and sediment critical diameter were modelled for four flood scenarios (2-year, 5-year, 20-year, and 50-year floods) as well as the extent of flooded area and additionally the grain size of channel sediment was calculated. The values of water depth, flow velocity, bed shear stress and sediment critical diameter increased significantly between 1964 and 1983, especially for 20-year and 50-year floods. Only the flow velocity within the floodplain zone did not increase for the two largest flood scenarios due to the expansion of riparian forest in the second half of the twentieth century. The increase in flow rate was accompanied by a progressive reduction of the extent of flooded area, especially between 1964 and 1983, as well as by increase in mean grain size of channel sediment. Between 1983 and 2012 changes in hydraulic parameters were less pronounced, and coarser and well packed channel sediment dominated on the river bed. Our work demonstrates that reconstruction of past river hydrodynamics, rather than river state at time horizons, can give essential insights into functioning of the river channel and floodplain during the intensification of human impacts after 1950s.

Dimensional effects in analysis of laser-induced-desorption diagnostics data

The accurate assessment of the local tritium concentration in the tokamak first wall by means of the laser-induced desorption (LID) diagnostics is sought as one the key solutions to monitoring the local radioactive tritium content in the first wall of the fusion reactor ITER. Numerical models of gas desorption from solids used for LID simulation are usually closed with the one-dimensional transport models. In this study, the temperature and particle dynamics in the target irradiated by a short laser pulse during LID are analyzed by means of the two-dimensional model to assess the validity of using one-dimensional approximation for recovering the diagnostics signal. The quantitative estimates for the parameters governing the heat and particle transfer are presented. The analytical expressions for the sample spatiotemporal temperature profiles driven by the target irradiation with a Gaussian laser beam with the trapezoid temporal shape are derived. The obtained relations are used to simulate tritium desorption from a tungsten sample driven by pulsed heating. It is shown that depending on the ratio between the laser spot radius and the heat diffusion length, the one-dimensional approach can noticeably overestimate the sample temperature in the limit of small laser spot radius (estimated for tungsten as ∼0.5–1.0 mm), resulting in more than 100% larger amounts of tritium desorbed from the target, compared to the two-dimensional approximation. In the limit of large laser spot radius (≥1.5 mm), both approaches yield comparable amounts of tritium desorbed from the sample.

Promoting precision surface irrigation through hydrodynamic modelling and microtopographic survey

Precision irrigation aims to deliver water and nutrients to crops at exactly the right time, in the right place and in the right amount. While surface irrigation is often perceived as less precise, accurate water distribution, wise use of resources and high efficiency can still be achieved with careful land preparation, astute irrigation management and rigorous performance monitoring. In this study, we advocate the innovative concept of Precision Surface Irrigation, centered around three key design and operational principles that are: well-organized field geometry (and microtopography), precise control of hydraulic-hydrological variables and, finally, regular performance evaluation. These pillars are then integrated into a simulation environment capable of capturing the intricacies of surface irrigation dynamics. Conducted in northern Italy's Padana plain, the study contrasts one-dimensional (WinSRFR) and two-dimensional (IrriSurf2D) irrigation dynamic modeling. Data collection includes boundary geometries, inflow rates, intervention durations, and microtopography, facilitating spatial performance assessment from both the models. The results show that the versatility of the two-dimensional modelling approach was able to reproduce well the observed water depths and the phases of water advance and depletion both in time and space within the studied border irrigation strips, even in complex situations where the strips were hydraulically connected. The RMSE between observed and simulated maximum water depth and waterfront advance time was less than 2.1 cm and 1.9 min, respectively. The two-dimensional approach was also able to detect the cross variability of irrigation dynamics, and to provide a spatial assessment of irrigation performance at high resolution. In conclusion, while the one-dimensional hydrodynamic approach to describing the hydraulic behavior of surface irrigation and field-scale irrigation performance remains valid, the two-dimensional approach provides, in our case study and reasonably elsewhere, a valid simulation environment for spatially characterizing irrigation dynamics in the context of Precision Surface Irrigation.

Determining the Size of Representative Volume Elements for a Two-Dimensional Random Aggregate Numerical Model of Asphalt Mortar without Damage

The size of the representative volume element (RVE) for the two-dimensional (2D) random aggregate numerical model of asphalt mortar in a non-destructive state, which directly affects the time required to simulate the linear viscoelastic behavior from asphalt mastic to asphalt mortar. However, in the existing literature, limited research has been conducted on the size determination of the numerical model RVE for asphalt mortar. To provide a recommended size for the typical 2D random aggregate numerical model RVE of asphalt mortar in a nondestructive state, this paper first applies the virtual specimen manufacturing method of asphalt concrete 2D random aggregate to asphalt mortar. Then, it generates numerical model RVEs of asphalt mortar with different maximum particle sizes, after which geometric and numerical analyses are conducted on these models. Finally, based on the geometric and numerical analysis results, the recommended minimum sizes of RVE for the 2D asphalt mortar numerical model are provided.

Determining the Size of Representative Volume Elements for a Two-Dimensional Random Aggregate Numerical Model of Asphalt Mortar without Damage

The size of the representative volume element (RVE) for the two-dimensional (2D) random aggregate numerical model of asphalt mortar in a non-destructive state, which directly affects the time required to simulate the linear viscoelastic behavior from asphalt mastic to asphalt mortar. However, in the existing literature, limited research has been conducted on the size determination of the numerical model RVE for asphalt mortar. To provide a recommended size for the typical 2D random aggregate numerical model RVE of asphalt mortar in a nondestructive state, this paper first applies the virtual specimen manufacturing method of asphalt concrete 2D random aggregate to asphalt mortar. Then, it generates numerical model RVEs of asphalt mortar with different maximum particle sizes, after which geometric and numerical analyses are conducted on these models. Finally, based on the geometric and numerical analysis results, the recommended minimum sizes of RVE for the 2D asphalt mortar numerical model are provided.

Competition between d-wave superconductivity and magnetism in uniaxially strained Sr2RuO4

The pairing symmetry of Sr2RuO4 is a long-standing fundamental question in the physics of superconducting materials with strong electronic correlations. We use the functional renormalization group to investigate the behavior of superconductivity under uniaxial strain in a two-dimensional realistic model of Sr2RuO4 obtained with density functional theory and incorporating the effect of spin-orbit coupling. We find a dominant dx2−y2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${d}_{{{{{\\rm{x}}}}}^{2}-{{{{\\rm{y}}}}}^{2}}$$\\end{document} superconductor mostly hosted by the dxy-orbital, with no other closely competing superconducting state. Within this framework, we reproduce the experimentally observed enhancement of the critical temperature under strain and propose a simple mechanism driven by the density of states to explain our findings. We also investigate the competition between superconductivity and spin-density wave ordering as a function of interaction strength. By comparing theory and experiment, we discuss constraints on a possible degenerate partner of the dx2−y2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${d}_{{{{{\\rm{x}}}}}^{2}-{{{{\\rm{y}}}}}^{2}}$$\\end{document} superconducting state.

Competition between d-wave superconductivity and magnetism in uniaxially strained Sr2RuO4

The pairing symmetry of Sr2RuO4 is a long-standing fundamental question in the physics of superconducting materials with strong electronic correlations. We use the functional renormalization group to investigate the behavior of superconductivity under uniaxial strain in a two-dimensional realistic model of Sr2RuO4 obtained with density functional theory and incorporating the effect of spin-orbit coupling. We find a dominant dx2−y2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${d}_{{{{{\\rm{x}}}}}^{2}-{{{{\\rm{y}}}}}^{2}}$$\\end{document} superconductor mostly hosted by the dxy-orbital, with no other closely competing superconducting state. Within this framework, we reproduce the experimentally observed enhancement of the critical temperature under strain and propose a simple mechanism driven by the density of states to explain our findings. We also investigate the competition between superconductivity and spin-density wave ordering as a function of interaction strength. By comparing theory and experiment, we discuss constraints on a possible degenerate partner of the dx2−y2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${d}_{{{{{\\rm{x}}}}}^{2}-{{{{\\rm{y}}}}}^{2}}$$\\end{document} superconducting state.

Numerical simulation of multiple hydraulic fracture propagation in heterogeneous coal reservoirs based on combined finite-discrete element method

Multi-stage fracturing in Horizontal well increases the permeability of coalbed methane by generating multiple fractures. However, the heterogeneity of coal reservoirs is a crucial factor that cannot be ignored in the study of multiple hydraulic fracture propagation. Therefore, we established a two-dimensional model for multiple hydraulic fracture propagation based on the combined finite-discrete element method (FDEM) and assigned a Weibull distribution function to the heterogeneity of the physical parameters of the cohesive elements in the model. The objective was to simulate and study the fracture propagation law of multi-cluster fracturing in horizontal wells in heterogeneous coal reservoirs. The research results indicated that: 1) as the heterogeneity of the coal reservoir weakened, the deflection angle of the main fracture increased. More secondary fractures were generated in the coal reservoir, leading to significant discontinuity. 2) The fracturing disturbance area was always concentrated at the tip of the main fracture, with a double wing shape. However, the fracturing disturbance areas at the tips of multiple main fractures could easily converge, with a square shape; 3) It is recommended to use a moderate injection rate and increase the perforation spacing appropriately when hydraulic fracturing is carried out in coal reservoirs with a heterogeneity coefficient m=5.

Emergence and evolution of electronic modes with temperature in spin-gapped Mott and Kondo insulators

Electronic modes emerge within the band gap at nonzero temperature in strongly correlated insulators such as Mott and Kondo insulators, exhibiting momentum-shifted magnetic dispersion relations from the band edges. As the temperature increases, the emergent modes gain considerable spectral weights and form robust bands that differ from the zero-temperature bands. Here, the origin of the emergent modes, their relation to doping-induced modes, and how their spectral weights increase with temperature are clarified in the ladder and bilayer Hubbard models and one- and two-dimensional Kondo lattice models using effective theory for weak inter-unit-cell hopping and numerical calculations. The results indicate that the temperature-driven change in the band structure reflecting spin excitation, including the change in the number of bands, can be observed in various strongly correlated insulators even with a spin gap, provided that the temperature increases up to about spin-excitation energies. The controlled analyses developed in this study significantly contribute to the fundamental understanding of the band structure of strongly correlated insulators at nonzero temperature. Published by the American Physical Society 2024

Numerical simulation of dynamics behavior of pulsed-DC helium plasma jet confined by parallel magnetic field at atmospheric pressure

A two-dimensional axisymmetric fluid model is used to simulate the dynamics behavior of an atmospheric-pressure helium plasma jet in the presence of a parallel magnetic field. The plasma jet is generated in a coaxial dielectric barrier discharge (DBD) driven by pulsed direct-current voltage. Comparative analysis of the plasma jet with and without the parallel magnetic field indicates that a slightly thinner plasma sheath inside the tube is present with the parallel magnetic field as a result of the decreased accumulated electrons on the inner surface of dielectric tube. After the streamer propagates outside the tube, a little more concentrated electron distribution in the annular wall is observed by applying the magnetic field because of the reduced electron diffusion in the radial direction and the confinement effect of the magnetic field on the electrons in the avalanche heads. The tiny reduction in the length of plasma jet is attributed to the E × B drift of charged particles. These results demonstrate that the parallel magnetic field has no apparent effect on the propagation of the plasma jet, and it contributes little to the performance improvement of the coaxial DBD, which agrees well with the previous experimental observations. This little impact of the parallel magnetic field on the coaxial DBD plasma jet may result from negligible contribution of the memory effect to the sole discharge pulse as well as the weak confinement effect of the applied magnetic field on the surface electrons that moves along the magnetic field lines under electrostatic repulsion. Published by the American Physical Society 2024

The molecular binding sequence transformation of soil organic matter and biochar dissolved black carbon antagonizes the transport of 2,4,6-trichlorophenol

Dissolved organic matter (DOM) and dissolved black carbon (DBC) are significant environmental factors that influence the transport of organic pollutants. However, the mechanisms by which their molecular diversity affects pollutant transport remain unclear. This study elucidates the molecular binding sequence and adsorption sites through which DOM/DBC compounds antagonize the transport of 2,4,6-trichlorophenol (TCP) using column experiments and modelling. DBC exhibits a high TCP adsorption rate (kn = 5.32 × 10−22 mol1–n∙Ln–1∙min−1) and conditional stability constant (logK = 5.19–5.74), indicating a strong binding affinity and antagonistic effect on TCP. This is attributed to the high relative content of lipid/protein compounds in DBC (25.65 % and 30.28 %, respectively). Moreover, the small molecule lipid compounds showed stronger TCP adsorption energy (Ead = −0.0071 eV/−0.0093 eV) in DOM/DBC, combined with two-dimensional correlation spectroscopy model found that DOM/DBC antagonized TCP transport in the environment through binding sequences that transformed from lipid/protein small molecule compounds to lignin/tannin compounds. This study used a multifaceted approach to comprehensively assess the impact of DOM/DBC on TCP transport. It reveals that the molecular diversity of DOM/DBC is a critical factor affecting pollutant transport, providing important insights into the environmental trend and ecological effects of pollutants.