To solve aircraft arrival sequencing and scheduling problems,and improve the typical predatory search algorithm(PSA),an innovative PSA is developed.The new PSA uses variable constraints of local search and global se...To solve aircraft arrival sequencing and scheduling problems,and improve the typical predatory search algorithm(PSA),an innovative PSA is developed.The new PSA uses variable constraints of local search and global search to avoid falling into local optimal solutions and the degeneration of solutions.To test the performance of new PSA,a case study with ten arriving flights and two runways is performed.Test results show that the new PSA performs much better than typical PSA and genetic algorithm(GA)in the aspects of the rate of gaining optimal solutions and the computational time.展开更多
The effect of supercooled melt forced laminar flow at low Reynolds Number on dendritic growth perpendicular to melt flow direction was investigated with the phase-field method by incorporating melt convection and ther...The effect of supercooled melt forced laminar flow at low Reynolds Number on dendritic growth perpendicular to melt flow direction was investigated with the phase-field method by incorporating melt convection and thermal noise under non-isothermal condition. By taking the dendritic growth of high pure succinonitrile (SCN) supercooled melt as an example, side-branching shape difference of melts with flow and without flow was analyzed. Relationships among supercooled melt inflow velocity, deflexion angle of dendritic arm and dendritic tip growth velocity were studied. Results show that the melt inflow velocity has few effects on the dendritic tip growth velocity. A formula of relationship between the velocity of the melt in front of primary dendritic tip and the dendritic growth time was deduced, and the calculated result was in quantitative agreement with the simulation result.展开更多
By utilizing phase field method combined with analysis on free energy and interatomic potentials, pre-precipitation phase formation and transformation process of Ni0.75Al0.05Fe0.2 alloy in early precipitation stage du...By utilizing phase field method combined with analysis on free energy and interatomic potentials, pre-precipitation phase formation and transformation process of Ni0.75Al0.05Fe0.2 alloy in early precipitation stage during the ageing process under 1 000 K were studied. And free energy, microstructures, compositions and volume fractions of pre-precipitation phase and equilibrium phase were analyzed. The simulation results indicate that nonstoichiometric Llo pre-precipitation phase formed first, and then would gradually transform into L12 equilibrium phase. It is discovered that the phase transformation process was closely related to free energy and interatomic potentials. Additionally, it is revealed that free energy of Llo pre-precipitation phase was higher and interatomic potential was smaller than that of L12 equilibrium phase. Therefore, it is concluded that Llo phase was unstable, and phase transformation would occur to L12 which was more stable.展开更多
Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field con...Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field controlled equations. And the alternating direction implicit(ADI) algorithm for solving temperature field controlled equation was also employed to avoid the restriction of time step. Some characteristics of the Ni-Cu alloy were captured in the process of non-isothermal solidification, and the comparative analysis of the isothermal and the non-isothermal solidification was investigated. The simulation results indicate that the non-isothermal model is favorable to simulate the real solidification process of binary alloys, and when the thermal diffusivity decreases, the non-isothermal phase-field model is gradually consistent with the isothermal phase-field model.展开更多
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton...Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation(capillary pressure).So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.展开更多
Small multi-turn coil devices are used with the transient electromagnetic method (TEM) in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and f...Small multi-turn coil devices are used with the transient electromagnetic method (TEM) in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.展开更多
The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled wi...The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.展开更多
With microscopic phase-field kinetic model, atomic-scale computer simulation program for the precipitation sequence and microstructure evolution of the ordered intermetallic compound γ' and θ in ternary Ni75AlxV25-...With microscopic phase-field kinetic model, atomic-scale computer simulation program for the precipitation sequence and microstructure evolution of the ordered intermetallic compound γ' and θ in ternary Ni75AlxV25-x alloy were studied. The simulation results show that Al concentration has important effects on the precipitation sequence. When Al concentration in Ni75AlxV25-x alloy is low, 0(Ni3V) ordered phase will be firstly precipitated, followed by γ'(Ni3Al) ordered phase. With Al concentration increasing, θ and γ' ordered phases are simultaneously precipitated. With A1 concentration further increasing, γ' ordered phase is firstly precipitated, followed by θ ordered phase. There is a competition relationship between θ and γ' ordered phases during growth and coarsening process. No matter which first precipitates, θ ordered phase always occupies advantage in the competition process of coarsening, thus, the microstructure with preferred orientation is formed.展开更多
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidificati...The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidification.Taking high pure succinonitrile (SCN) undercooled melt as an example,the forced laminar flow effect on multiple grains was studied by phase-field model of single grain which coupled with flow equations at non-isothermal condition.The simulation results show that the random grid method can reasonably settle the problem of random distribution and is more effective.When the solid fraction is relatively low,melt particles flow around the downstream side of dendrite,and the flow velocity between two dendrite arms becomes high.At the stage of solidification time less than 1800Δt,every dendrite grows freely;the upstream dendrites are stronger than the downstream ones.The higher the melt flow rate,the higher the solid fraction.However,when the solid fraction is relatively high,the dendrite arm intertwins and only a little residual melt which is not encapsulated can flow;the solid fraction will gradually tend to equal to solid fraction of melt without flow.展开更多
文摘To solve aircraft arrival sequencing and scheduling problems,and improve the typical predatory search algorithm(PSA),an innovative PSA is developed.The new PSA uses variable constraints of local search and global search to avoid falling into local optimal solutions and the degeneration of solutions.To test the performance of new PSA,a case study with ten arriving flights and two runways is performed.Test results show that the new PSA performs much better than typical PSA and genetic algorithm(GA)in the aspects of the rate of gaining optimal solutions and the computational time.
基金Project (10964004) supported by the National Natural Science Foundation of ChinaProject (096RJZA104) supported by the Natural Science Foundation of Gansu Province, China
文摘The effect of supercooled melt forced laminar flow at low Reynolds Number on dendritic growth perpendicular to melt flow direction was investigated with the phase-field method by incorporating melt convection and thermal noise under non-isothermal condition. By taking the dendritic growth of high pure succinonitrile (SCN) supercooled melt as an example, side-branching shape difference of melts with flow and without flow was analyzed. Relationships among supercooled melt inflow velocity, deflexion angle of dendritic arm and dendritic tip growth velocity were studied. Results show that the melt inflow velocity has few effects on the dendritic tip growth velocity. A formula of relationship between the velocity of the melt in front of primary dendritic tip and the dendritic growth time was deduced, and the calculated result was in quantitative agreement with the simulation result.
基金Projects(10902086, 50941020, 50875217) supported by the National Natural Science Foundation of China Projects(JC201005) supported by Basic Research Fund of Northwestern Polytechnical University, China Project supported by Graduate Starting Seed Fund and Doctoral Foundation of Northwestern Polytechnical University, China
文摘By utilizing phase field method combined with analysis on free energy and interatomic potentials, pre-precipitation phase formation and transformation process of Ni0.75Al0.05Fe0.2 alloy in early precipitation stage during the ageing process under 1 000 K were studied. And free energy, microstructures, compositions and volume fractions of pre-precipitation phase and equilibrium phase were analyzed. The simulation results indicate that nonstoichiometric Llo pre-precipitation phase formed first, and then would gradually transform into L12 equilibrium phase. It is discovered that the phase transformation process was closely related to free energy and interatomic potentials. Additionally, it is revealed that free energy of Llo pre-precipitation phase was higher and interatomic potential was smaller than that of L12 equilibrium phase. Therefore, it is concluded that Llo phase was unstable, and phase transformation would occur to L12 which was more stable.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of China
文摘Based on the entropy function, a two-dimensional phase field model of binary alloys was established. Meanwhile, an explicit difference method with uniform grid was adopted to solve the phase field and solute field controlled equations. And the alternating direction implicit(ADI) algorithm for solving temperature field controlled equation was also employed to avoid the restriction of time step. Some characteristics of the Ni-Cu alloy were captured in the process of non-isothermal solidification, and the comparative analysis of the isothermal and the non-isothermal solidification was investigated. The simulation results indicate that the non-isothermal model is favorable to simulate the real solidification process of binary alloys, and when the thermal diffusivity decreases, the non-isothermal phase-field model is gradually consistent with the isothermal phase-field model.
基金supported by the 973 Program (Grant No.2007CB209505)the National Natural Science Foundation of China (Grant No.40674061,40704019)
文摘Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation(capillary pressure).So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.
基金supported by the National Natural Science Foundation of China(Nos.41564001 and 41572185)the Natural Science Foundation of Jiangxi Province(No.20151BAB203045)
文摘Small multi-turn coil devices are used with the transient electromagnetic method (TEM) in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.
基金Project (10964004) supported by the National Natural Science Foundation of ChinaProject (20070731001) supported by Research Fund for the Doctoral Program of ChinaProject (096RJZA104) supported by the Natural Science Foundation of Gansu Province,China
文摘The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.
基金Projects(51174168,51274167)supported by the National Natural Science Foundation of ChinaProject(2014JM7261)supported by the Natural Science Basic Research Plan in Shaanxi Province of ChinaProject(B08040)supported by "111" Project,China
文摘With microscopic phase-field kinetic model, atomic-scale computer simulation program for the precipitation sequence and microstructure evolution of the ordered intermetallic compound γ' and θ in ternary Ni75AlxV25-x alloy were studied. The simulation results show that Al concentration has important effects on the precipitation sequence. When Al concentration in Ni75AlxV25-x alloy is low, 0(Ni3V) ordered phase will be firstly precipitated, followed by γ'(Ni3Al) ordered phase. With Al concentration increasing, θ and γ' ordered phases are simultaneously precipitated. With A1 concentration further increasing, γ' ordered phase is firstly precipitated, followed by θ ordered phase. There is a competition relationship between θ and γ' ordered phases during growth and coarsening process. No matter which first precipitates, θ ordered phase always occupies advantage in the competition process of coarsening, thus, the microstructure with preferred orientation is formed.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
基金Project(10964004) supported by the National Natural Science Foundation of ChinaProject(20070731001) supported by Research Fund for the Doctoral Program of China+1 种基金 Project(096RJZA104) supported by the Natural Science Foundation of Gansu Province,ChinaProject(SB14200801) supported by the Doctoral Fund of Lanzhou University of Technology,China
文摘The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidification.Taking high pure succinonitrile (SCN) undercooled melt as an example,the forced laminar flow effect on multiple grains was studied by phase-field model of single grain which coupled with flow equations at non-isothermal condition.The simulation results show that the random grid method can reasonably settle the problem of random distribution and is more effective.When the solid fraction is relatively low,melt particles flow around the downstream side of dendrite,and the flow velocity between two dendrite arms becomes high.At the stage of solidification time less than 1800Δt,every dendrite grows freely;the upstream dendrites are stronger than the downstream ones.The higher the melt flow rate,the higher the solid fraction.However,when the solid fraction is relatively high,the dendrite arm intertwins and only a little residual melt which is not encapsulated can flow;the solid fraction will gradually tend to equal to solid fraction of melt without flow.