Two-dimensional Numerical Model Research Articles

Mechanism study on the synergistic coupling of the double-dish solar latent heat storage system to enhance heat transfer

With the rapid advance of the global economy, latent heat storage (LHS) is critical to solar thermal utilization. In this study, a double-dish solar Stirling LHS power generation system was designed. The heat transfer performance of the thermal storage system was improved by using the gradient tree-shaped fins and adding graphene nanoparticles to paraffin. Two-dimensional numerical models of three thermal storage systems were developed to study the effects of natural convection and graphene nanoparticle concentration on the melting characteristics of paraffin. From the perspective of field synergy, the coupling effect of fin structure and graphene nanoparticle concentration on the heat transfer process was analyzed. The results show that natural convection exerts an important role in the heat transfer characteristics of paraffin melting. Compared with the six-longitudinal and snowflake fin structures, the gradient tree-shaped fin structure shows a 52.13% and 27.73% reduction in the complete melting time and a 52.16% and 9.53% increase in thermal storage efficiency. Little variation is observed in the liquid phase rate and the average temperature of the composite phase change materials (PCMs) with different concentrations of graphene nanoparticles. The total thermal storage capacity of composite PCMs is higher than that of pure paraffin. The synergistic coupling of gradient tree-shaped fin arrangement and graphene nanoparticles added in paraffin can enhance heat transfer.

Volumic Eddy-Current Losses in Conductive Massive Parts with Experimental Validations

In this paper, the analytical determination of volumic eddy-current losses in rectangular-shaped conductive massive parts is presented with experimental validations. Eddy currents, as well as the resulting volumic losses, are generated by a sinusoidal spatially uniform applied magnetic field. A U-shaped electromagnetic device with a flat mobile armature (or adjustable air gap) is used to measure the eddy-current losses. The experimental device, its instrumentation, and the conductive massive parts are presented in detail in the paper. Thereafter, the magnetic field distribution applied on the conductive massive parts, which is the mean input data for the eddy-current loss model, is studied. A two-dimensional (2D) numerical model, under the FEMM software, for the magnetic field calculation was also developed. A comparative analysis between the experimental measurements and the numerical results allowed the distribution of the applied magnetic field to be accurately validated. In the final phase, the objective was to estimate the volumic eddy-current losses in rectangular-shaped aluminium conductive massive parts generated by the variation of this applied magnetic field. An analytical model, based on the Maxwell–Fourier method, for the accurate prediction of eddy-current losses has also been developed. An electromagnet with and without the conductive massive parts is characterized in terms of power consumption. By using the power conservation method (i.e., Boucherot’s theorem), the eddy-current losses could be quantified experimentally. The influence of segmentation is also studied. The analytical results are compared to the experimental test results.

Numerical Analysis on Spanwise Correlation of Vortex-Induced Force of Split Double-Box Beam

The vortex-induced force of bridges is not synchronized along the span direction. The spatial correlation of vortex-induced force becomes obvious when the span of the bridge structure gradually increases, and the mass gradually decreases. Therefore, it is important to study the spanwise correlation of vortex-induced forces of long-span bridges to ensure the accuracy of the prediction of vortex-induced vibration response of bridges. The study established and improved the theory and experimental research methods for vortex-induced vibration analysis of large-span bridges by discussing the spanwise correlation during vortex-induced vibration of the split double-box beam. Taking Xihoumen Bridge as the object, two-dimensional (2D) and three-dimensional (3D) numerical models of the scaled-down sections were designed and established based on the ANSYS Fluent platform and the RANS SST turbulence model. Based on the Newmark-β algorithm, a User Defined Function (UDF) was written for vortex vibration calculation, the three-dimensional bypass of the split double-box beam in static conditions and vortex-induced vibration of the split double-box beam of which were calculated, and the spanwise correlations of the aerodynamic coefficients, the surface pressure coefficients, and the wake wind velocity of the main girder section were analyzed for the static and vibration states of the bridge. The results show that the self-excited force of the split double-box beam in vortex-induced vibration improved its spanwise correlation. Compared with those in a static state, the spanwise correlation of the lift coefficient and torque coefficient increased by 55% and 87%, respectively, and the resistance coefficient increased by more than 10 times. The correlation of the pressure coefficients increased by 153%. The correlation of wake wind velocity increased by 37% in the along-wind direction and that in the across-wind and vertical-wind direction increased by more than 10 times. The accuracy of the numerical simulation results was verified by comparing the pressure distribution and pressure spanwise correlation of the main beam with field-measured data.

Dynamic Evolution Modeling of a Lake-Terminating Glacier in the Western Himalayas Using a Two-Dimensional Higher-Order Flowline Model

To better understand the future evolution of Jiemayangzong Glacier (JMYZG), the headstream of the Yarlung Zangbo River, we simulated its future ice thickness evolution using a two-dimensional higher-order numerical flowline model. Due to the sparsity of in situ observational data, we used a combination of field observations and inversion models of velocity and ice thickness to initialize the model parameters. We validated the parametrizations of the calving scheme by comparing the modeled and observed glacier terminus retreats. To estimate the response of JMYZG to climate change, the ice flow model was forced with different climate scenarios. We found that the JMYZG will retreat under different climate scenarios. By 2100, the volume loss of JMYZG will be approximately 34%, 67% and 81% under SSP1-2.6, SSP3-7.0 and SSP5-8.5, respectively.

Numerical Simulation of Wave–Current Force Characteristics of Horizontal Floating Cylinder in Heave Motion

This paper presents the characteristics of the heave motion responses and hydrodynamic forces of a horizontal floating circular cylinder during a wave–current interaction. A two-dimensional numerical model based on the CFD (Computational Fluid Dynamics) method for modeling wave flow is validated and verified. The hydrodynamic characteristics of the horizontal floating cylinder during heave motion were calculated and analyzed under the conditions of different k values (stiffness of spring), wave amplitudes, submerged depths, and flow rates. The results show that, with the increase in the k value, the vibration amplitude of the cylinder first increases and then decreases. The vibration amplitude peak is achieved, the vibration frequency is consistent with the wave frequency, and a resonant motion takes place. When the wave amplitude and flow rate are fixed, the maximum vibration amplitude decreases as the wave period increases. When the cylinder is half-submerged, the effect of the current on motion is significant; the vibration amplitude is less than the wave amplitude. When a quarter is submerged, the vibration amplitude is larger than that of the half-submerged cylinder at each k. The maximum amplitude is greater than the wave amplitude, and the vibration amplitude reaches the minimum at the moderate flow rate for each k.

Thermal control of lithium-ion batteries for electric cars by metal foam partially filled with Phase Change Material

A viable solution to environmental problems related to pollutant emissions from internal combustion engine vehicles is electric cars. The main problem with electric vehicles, which also limits their widespread use, is due to lithium batteries and in particular their thermal control. Thermal control of lithium batteries is an extremely topical issue since it is not only related to the performance of the vehicles on which they are mounted, but also to the safety of the vehicles and drivers. For this purpose, both factors internal to the battery itself and external causes such as environmental conditions must be considered.In this work, a thermal control system based on metal foams partially filled by the phase change material (PCM) was numerically investigated as a possible application for cooling lithium batteries. For this purpose, a two-dimensional numerical model was constructed in which the PCM partially fills the metal foam. The finite volume method used through the ANSYS Fluent code was used for the governing equations, which are written assuming local thermal equilibrium. Several cases were simulated for different values of the external convective heat transfer coefficient. The results obtained, are reported in terms of temperatures and liquid fractions. Comparisons with partially filled and totally filled cases by PCM in the thermal control system are reported to show the advantage obtained, in terms of thermal control, with a composite system of metal foam and PCM.

CFD modelling of spatiotemporal evolution of detailed tar species in a downdraft gasifier

The rising demand for renewable energy around the world has sparked interest in biomass gasification. However, the technology greatly suffers because of tar species produced during the gasification process, which limits direct use of the produced gas. To address this issue, the paper presents a novel piece of work that focuses on the formation and evolution of tar species consisting of benzene, naphthalene, toluene, and phenol. A two-dimensional numerical model for a downdraft biomass gasifier is developed with a total of 20 thermochemical kinetic reactions to investigate the formation of tar species in the gasifier with the effect of residence time. The model's predictions are validated with the experimental and kinetic data and found to be in good agreement. Besides, the model's ability to simulate the producer gas production from a downdraft gasifier is examined. Reaction rates for volatiles decomposition, combustion, and gasification reactions under different working conditions are investigated. Overall, benzene has the highest concentration of the selected tar species, followed by naphthalene, and with relatively modest amounts of phenol and toluene.

Accuracy of Seismic Response Evaluation of Two-Dimensional Analysis Model with Rigid Joints for RC Frame Buildings.

Three- or two-dimensional (2D) numerical models are used for the evaluation of the seismic performance of reinforced concrete (RC) buildings. This study examines a 2D numerical model for a specimen used in a full-scale four-story RC shaking-table test and evaluates the accuracy of the seismic response of the 2D numerical model, which is composed of a square fiber section model for the columns, a T-shape fiber section model for the beam and slab, and a rigid joint model for the beam-column joint. A parametric analysis of the effective slab width is performed to analyze its effects on the modal shape and natural period. The results suggest that the primary natural period of the considered model is almost identical to that associated with the experimental results. The applicability of the 2D numerical model for estimating the seismic response of the structure is established. By comparing the results of the seismic analysis and the experiment in the 50% amplitude of the JMA-Kobe wave, which corresponds to slightly exceeding VII on the modified Mercalli intensity scale, the root-mean-square percentage error of the 2D numerical model is 1.03% for the floor acceleration and 4.7% for the inter-story drift. Thus, the analytical model used in this study has sufficient accuracy in evaluating the seismic performance of buildings constructed in regions with a maximum seismic intensity of VII.

Soil-structure interaction of integral abutments

Integral abutment bridges (IABs) are pioneering structures that offer certain advantages over conventional bridges. The main difference is the ability of IABs to eliminate the necessity of utilizing bearings and expansion joints. This results in more economical construction and enhanced structural performance. However, the intricate soil-structure interactions of IABs under the influence of seasonal changes are not yet fully understood. This paper aims to improve the understanding of these innovative structures by utilizing comprehensive numerical models to investigate their geostructural performance under various conditions. The case study selected for this research was the Middlesex bridge in Vermont, USA. A full-scale monitoring program at this bridge measured and recorded the bridge acting pressures, deformations, and internal forces, making it possible to gain insight into the thermal response of such structures. For the numerical research, a comprehensive two-dimensional numerical model, developed by using the PLAXIS software, was verified against the available field measurements. Subsequently, parametric studies were conducted to investigate earth pressure and pile bending moment variations influenced by different model conditions. Key model parameters were varied, including the constitutive soil model, thermal loading, backfill stiffness, pile size and orientation, and span length. It was found that utilizing a linear constitutive soil model could lead to significant inaccuracies in the results. Increasing the backfill stiffness was seen to increase passive earth pressures and decrease pile bending moments. With smaller pile sections, oriented for weak-axis bending, pile bending moments decreased substantially and lateral earth pressures increased. Because dead load magnitudes were closely linked to span lengths, pile bending moments and backfill pressures increased when fewer and longer bridge spans were used.

Modeling the hydrodynamics of a tidal inlet during bay-side storms

Bay-side storms, defined here as storms with tracks on the landward side of barrier islands, may disturb the hydrodynamics of inner bays to a larger extent than on the ocean side. These storms are common in large-scale O(>100,000 m) estuarine systems and have the potential to modify the circulation in bays and within tidal inlets. Here, we provide an overview of the hydrodynamic response of a tidal inlet under forcings caused by bay-side storms and explore the role of waves in modulating the release of storm surge from the back-barrier regions into the ocean. A two-dimensional horizontal numerical model including wave-current interactions is calibrated and validated against field observations of water levels and depth-averaged velocities at Oregon Inlet, NC. The model is then used to investigate the effect of synthetic bay-side storms with varying wave conditions and water levels based on those generated by Hurricane Irene (2011), which is the strongest bay-side storm to hit the Outer Banks of North Carolina in the last two decades. Effect of timing of the peak storm surge during the ebb and flood phases of the tide is also explored.Results from synthetic storms indicate that, during bay-side storms, the water level gradient along the inlet favors ebbing flows regardless of the timing of the storms relative to tidal phase. These results suggest that waves might be responsible for any influx of volume to the bay during high bay-side surge events. Wave blocking effects were found to be stronger along the ebb shoal and only reached the flood delta when bay water levels were nearly the same as the ocean water levels. Reduction of currents by waves in the inlet have the potential to extend the duration of the inundation period in the back barrier region. Bay-side storms also caused flux enhancement over inlet shoals and channels in the flood delta, which could have implications for circulation patterns as well as the morphodynamics of the system.

Effects of Different Heat Transfer Conditions on the Hydrogen Desorption Performance of a Metal Hydride Hydrogen Storage Tank

To investigate the influence of thermal effects on the hydrogen desorption performance of the metal hydride hydrogen storage system, a two-dimensional numerical model was established based on a small metal hydride hydrogen storage tank, and its accuracy was verified by the temperature variations in the reaction zone of the hydrogen storage tank during hydrogen desorption. In addition, the influence of the heat transfer medium on the heat and mass transfer performance of the hydrogen desorption reaction was analyzed. An external heat transfer bath was added to simulate the thermal effect of the model during the hydrogen desorption reaction. The temperature and type of heat transfer medium in the heat transfer bath were modified, and the temperature and reaction fraction variations in each zone of the hydrogen storage model were analyzed. The results showed that under heat transfer water flow, the reaction rate in the center region of the hydrogen storage tank was gradually lower than that in the wall region. The higher the temperature of water flow, the shorter the total time required for the hydrogen desorption reaction and the shortening amplitude is reduced. The variations in the temperature and hydrogen storage capacity during hydrogen desorption were similar, with water and oil as the heat transfer medium, under the same flow rate and heat transfer temperature, however, the heat transfer time and hydrogen desorption time of water were about 10% and 5% shorter than that of oil, respectively. When the air was used as the heat transfer medium, the heat transfer rate of the air convection in the channel was lower than the heat transfer rate of the tank wall, reducing the temperature difference between the air and alloy on both sides of the wall, decreasing heat transfer efficiency, and significantly prolonging the time required for hydrogen desorption.

Numerical investigation of solute transport into deformable marine sediments driven by ocean waves

Most existing studies do not consider the effect that the deformability of submarine sediments may have on solute transport under the action of ocean waves. In this study a two-dimensional numerical model based on Biot’s poro-elastic consolidation theory and an advection–dispersion equation is established to investigate the transport of non-reactive solute into a deformable porous seabed under the influence of ocean waves. The numerical model is then validated against experimental and analytical results in order to demonstrate its accuracy and robustness. The simulation results indicate that ignoring the deformability of seabed soil and associated seepage force due to ocean waves may lead to a significant underestimation of the rate of solute transport in shallow layers of marine sediment, particularly fine sand or silt. The effects of parameters such as the shear modulus, non-dimensional mechanical dispersion coefficient (related to particle size and permeability), dimensionless sediment thickness, wave steepness and the relative water depth have on solute transport into deformable sediments were investigated. The results indicate that the ratio of the longitudinal coefficient of mechanical dispersion between deformable and non-deformable sediments that can represent the enhancement of solute migration due to soil deformation may reach up to 85. This becomes more pronounced as the shear modulus and permeability of sediments decrease and the depth of water increases when the sediment thickness is less than 0.8 times the wavelength. The rate of solute transport in deformable sediments reaches its maximum when the sediment thickness is about 0.2 times the wavelength. In fact, as the relative water depth increases from 0.1 to 0.3, the enhancement of solute migration due to wave-induced soil deformation increases by about 60%-80%. In comparison to linear waves, shallow water waves, e.g., the first-order cnoidal waves, can increase the transport of solute more due to its strong non-linear behaviour, particularly in deformable seabed.

Role of Hydrothermal Circulation along and above Inherited Basement Structures Relating to Unconformity-Related Uranium Mineralization

The reactivation of inherited tectonic structures formed during the Paleoproterozoic Trans-Hudson Orogeny (THO) has played a significant role in generating high-grade unconformity-related uranium deposits in the eastern Athabasca Basin. The role of these tectonic structures is now investigated through a series of two-dimensional hydrothermal numerical models. Two modelling scenarios are considered: (1) models during the THO peak of metamorphism and (2) models with a permeable layer mimicking the presence of the Athabasca Basin, deposited unconformably over the THO basement. In the first scenario, general fluid patterns are strongly affected by the applied permeability configurations. Unidirectional high fluid flow zones (from 10-9 to 10-8 m·s-1) and high thermal gradients (up to 65°C·km-1) can be observed above and within the deep-seated tectonic structures. In the second scenario, well-established fluid convection cells or unidirectional fluid flow zones are observed within the basin layer, with upflow originating from the core of the deep-seated structures, regardless of the applied fluid pressure regime. These results highlight that these deep-seated structures can efficiently transport fluids and heat towards the upper parts of the crust and the basin. In the second scenario, the loci for preore alteration are then evaluated by computing a rock alteration index based on temperature and fluid velocity constraints. These alteration areas reside along and above the deep-seated structures and are potential regions for structural reactivation during mineralization. These results imply that the analysis of the inherited tectonic structures, combined with the alteration regions, can serve as markers for uranium exploration.

Numerical simulation of combustion surface regression based on Butterworth filter in hybrid rocket motor

This paper provides a valuable method to perform transient numerical simulations for hybrid rocket motors (HRMs). A two-dimensional transient numerical simulation model based on Butterworth filter and dynamic mesh is established. The Butterworth filter adopted in this project is a second-order low-pass filter. Then a firing test is arranged to prove the accuracy of the numerical simulation model of the HRMs. This rocket motor adopts 98% hydrogen peroxide (HP) and hydroxyl-terminated polybutadiene (HTPB) fuel containing aluminum as the propellants. The results demonstrate that the Butterworth filter can significantly improve the smoothness of the combustion surface compared with non-filter. And the combustion chamber pressure and thrust of simulation and test results are in good agreement. The average deviation of combustion chamber pressure between the simulation and test is 0.9%, and the average deviation of thrust between the simulation and test is 3.3%, which verifies the feasibility of the numerical simulation model for this HRM. Comparing the regression rate and filtered regression rate along with axial distance over time, it is indicated that the high-frequency oscillations are filtered by the Butterworth filter. In order to compare the port of simulation and test, the combustion chamber is scanned by computed tomography (CT) after the test. The results demonstrate that the simulation and test results are in good agreement, and the average deviation of simulation and test results is 4.5% mainly due to the end-burning of grain and nitrogen blowing. Then, the long-time simulations that adopted the Butterworth filter and non-filter were conducted. The results demonstrate that the use of the Butterworth filter can significantly smooth the combustion surface and the model has the ability to perform long-time simulations. This project can make the simulation for HRMs closer to the test results and provide a valuable guide for HRMs simulation.

Effect of pyrolysis kinetic parameters on the overload ignition of polymer insulated wires in microgravity

The overload ignition of polymer insulated wires in microgravity was investigated as basic information for space fire safety. This paper focuses on the effect of pyrolysis kinetic parameters of the insulation material on the ignition characteristics. A two-dimensional numerical model calculation under varied pyrolysis model was carried out to simulate the overload ignition, and microgravity experiments have been conducted for the model validation. With the help of thermogravimetric analysis, the original pyrolysis model of the polyethylene insulation used for wire ignition experiments was obtained. Three groups of varied pyrolysis models were proposed by changing the activation energy (E), pre-exponential factor (A) from the original model. Numerical simulations were conducted with the varied pyrolysis models. The simulated results employing the original pyrolysis model showed a good agreement with the microgravity experiments for ignition delay time with continuous current and the minimum ignition energy (MIE) with short-term current, respectively. From the simulation results with continuous overloading, the independent variation of A and E in the pyrolysis model has a limited effect on the ignition delay time but a significant effect on ignition modes in microgravity. Such effects are mainly reflected by changing the pyrolysis temperature during overloading, which determines the ease of spontaneous ignition. With a constant pyrolysis temperature, the change of mass loss rate in the pyrolysis model may affect the transition point from the “spontaneous ignition” to the “wire assisted ignition” when the current increases. Finally, for short-term current supply in microgravity, the single change of A or E showed higher MIE than the original pyrolysis model.

Simulation of Inclusion Particle Motion Behavior under Interfacial Tension in Continuous Casting Mold

Inclusions entrapped by the solidifying front during continuous casting adversely affect the properties of the final steel products. In this study, we investigated the effect of the interfacial tension due to surfactant concentration, particularly sulfur, on alumina inclusion motion behavior during molten steel solidification in a continuous casting mold. A two-dimensional numerical model was developed in Ansys Fluent software to simulate the inclusion motion in a continuous casting mold. Further, the impacts of different values of the alumina inclusion diameter, sulfur concentration, and melt temperature were studied to understand the inclusion motion behavior. The inclusion diameter affected the inclusion distribution throughout the domain. The alumina inclusion entrapment percentage varied in the case of sulfur mixing (using an empirical relationship for modeling). It was found that the removal percentage varied according to the sulfur concentration. The addition of sulfur at concentrations from 10 ppm to 70 ppm resulted in a 4% increase in the removal of alumina inclusions (trapped in the solidifying shell), except for the 100-ppm case. Smaller-sized inclusion particles had a 25% higher chance of entrapment at the top level of the mold. Under the effect of a higher surface tension gradient between inclusions and the melt, the predicted findings show that inclusions were vulnerable to engulfment by the solidification front.

Numerical analysis of inner heating tube position for improving solid-phase transition in a shell-and-tube heat accumulator

Latent heat thermal storage (LHTS) system is vital to reduce environment pollution. In the shell-and-tube heat accumulator, the position of the inner heating tube plays a vital role in the thermal storage. To analyze the effect of the inner tube position on the phase transition, a two-dimensional numerical model is developed. The structure has the minimum full melting time of 3480 s when the inner tube is 12 mm (L = 12 mm) from the center. Compared with L = 0 mm, the full melting time at L = 12 mm can be reduced by 13.4%.

Wave-triggered breakup in the marginal ice zone generates lognormal floe size distributions: a simulation study

Abstract. Fragmentation of the sea ice cover by ocean waves is an important mechanism impacting ice evolution. Fractured ice is more sensitive to melt, leading to a local reduction in ice concentration, facilitating wave propagation. A positive feedback loop, accelerating sea ice retreat, is then introduced. Despite recent efforts to incorporate this process and the resulting floe size distribution (FSD) into the sea ice components of global climate models (GCMs), the physics governing ice breakup under wave action remains poorly understood and its parametrisation highly simplified. We propose a two-dimensional numerical model of wave-induced sea ice breakup to estimate the FSD resulting from repeated fracture events. This model, based on linear water wave theory and visco-elastic sea ice rheology, solves for the scattering of an incoming time-harmonic wave by the ice cover and derives the corresponding strain field. Fracture occurs when the strain exceeds an empirical threshold. The geometry is then updated for the next iteration of the breakup procedure. The resulting FSD is analysed for both monochromatic and polychromatic forcings. For the latter results, FSDs obtained for discrete frequencies are combined following a prescribed wave spectrum. We find that under realistic wave forcing, lognormal FSDs emerge consistently in a large variety of model configurations. Care is taken to evaluate the statistical significance of this finding. This result contrasts with the power law FSD behaviour often assumed by modellers. We discuss the properties of these modelled distributions with respect to the ice rheological properties and the forcing waves. The projected output can be used to improve empirical parametrisations used to couple sea ice and ocean wave GCM components.

Novel electrode design having gradually increasing porosity in a vanadium redox flow battery

Vanadium redox flow battery systems safely and efficiently store energy obtained from conventional energy sources (e.g. natural gas, petroleum and coal, and so on.) and renewable energy sources (e.g. solar, and, so on.). These battery systems need to be analyzed in depth since they are one of the most attractive energy storage technology. This study numerically investigates the effects of gradually increasing electrode porosity on potential performance, overpotentials, current densities, species concentrations, and pressure distributions in the VRFB during discharge process. In order to reveal these effects, two-dimensional and steady-state numerical models composed of four different cases are solved. Electrode porosity is constant at 94% for both electrodes in Case 1. While the negative electrode porosity is gradually increased from 64% to 94% (step by 10%) along the electrode thickness towards the membrane in Case 2, the positive electrode porosity is gradually changed from 64% to 94% along the same direction in Case 3. The other electrode porosity is kept constant at 94% in Case 2 and Case 3. Case 4 has electrodes having gradually increasing porosity in both half-cells. The numerical results show that Case 4 has higher potential performance, lower overpotential, and more uniform current distributions than the other cases. Furthermore, there are slight changes in the concentration distributions of consumed and produced species in all cases. Moreover, pressure drop in the region having 64% porosity is too much in cases applied gradually increasing porosity.

Fabricating and determining representative volume elements of two-dimensional random aggregate numerical model for asphalt concrete without damage

In order to optimize the fabricating method of representative volume element (RVE) for two-dimensional (2D) random aggregate numerical model of asphalt concrete, and determine the minimum size of RVE corresponding to the 2D random aggregate numerical model of typical graded asphalt concrete in the nondestructive state. The authors first proposed an RVE fabricating method that includes “aggregate region identification” algorithm and can place “partial aggregate”. Then, by using this method, 9 different sizes of RVEs corresponding to AC-13C, AC-20C and AC-25C asphalt concrete (AC-13C, AC-20C and AC-25C asphalt concrete are dense-graded asphalt mixtures, and their nominal maximum aggregate sizes are 13.2 mm, 19 mm and 26.5 mm respectively) were fabricated. Subsequently, the geometric characteristics (aggregate area proportion, aggregate angle, “density” homogeneity) of these RVEs were analyzed. After that, virtual dynamic modulus tests were performed on these RVEs. Finally, the results of geometric analysis and virtual dynamic modulus test were comprehensively analyzed. The results show that the “aggregate region identification” algorithm can effectively improve the aggregate placement efficiency. By placing “partial aggregate”, the RVEs with small-sizes that meet the requirements of aggregate area proportion in each particle size range can be fabricated. In the nondestructive state, the minimum sizes of 2D random aggregate numerical model’s RVE corresponding to AC-13C, AC-20C and AC-25C asphalt concrete are 20 mm, 30 mm and 40 mm respectively.