Two-Dimensional Static Numerical Modeling and Simulation of AlGaN/GaN HEMT

AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered in this model. Then the two-dimensional conduction band and electron distribution, electron temperature characteristics, Id versus Vd and Id versus Vg, transfer characteristics and transconductance curves are obtained. Corresponding analysis and discussion based on the simulation results are subsequently given.

Simulation of Horizontal-Two-Dimension Focused Waves Using A Two-Layer Boussinesq-Type Model

Accurate simulation of the horizontal-two-dimension(H2D)focused wave group in deep water requires high accuracy of a numerical model.The two-layer Boussinesq-type model(Liu and Fang,2016;Liu et al.,2018)with the highest spatial derivative of 2 has high accuracy in both linear and nonlinear properties.Based on the further development of the velocity equations(Liu et al.,2023),the H2D numerical model for water waves is established with the prediction-correction-iteration model in the finite difference method,and a composite fourth-order Adams-Bashforth-Moulton scheme is used for time integration.The wave generation method proposed by Hsiao et al.(2005)is applied and calibrated in this H2D model.The numerical calculations lead to the following three main conclusions:First,compared with the analytical solution of Stokes linear waves,the calculated velocity profiles show higher accuracy by using the improved velocity formulas.Second,the simulations of the focused multidirectional wave group are carried out,and good agreements are found,demonstrating that the present H2D numerical model shows high accuracy in simulating focused multidirectional wave groups,and the effectiveness of the improved velocity formulas is also validated.Furthermore,the velocity profiles throughout the computational domain at the time of maximum wave crest are given.Finally,the FFT method is used to obtain the amplitude with different frequencies for several locations,and the changes of the wavelet energy spectrum at different locations are presented for several cases.

Development and application of a two-dimensional water quality model for the Daqinghe River Mouth of the Dianchi Lake

Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Performance assessment of two-dimensional hydraulic models for generation of flood inundation maps in mountain river basins

Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.

Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition

A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.

Comprehensive two-dimensional river ice model based on boundary-fitted coordinate transformation method

River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.

Numerical Simulation for Two-Phase Water Hammer Flows in Pipe by Quasi-Two-Dimensional Model

The features of a quasi-two-dimensional( quasi-2D) model for simulating two-phase water hammer flows with vaporous cavity in a pipe are investigated. The quasi-2D model with discrete vaporous cavity in the pipe is proposed in this paper. This model uses the quasi-2D model for pure liquid zone and one-dimensional( 1D) discrete vapor cavity model for vaporous cavity zone. The quasi-2D model solves two-dimensional equations for both axial and radial velocities and 1D equations for both pressure head and discharge by the method of characteristics. The 1D discrete vapor cavity model is used to simulate the vaporous cavity occurred when the pressure in the local pipe is lower than the vapor pressure of the liquid. The proposed model is used to simulate two-phase water flows caused by the rapid downstream valve closure in a reservoir-pipe-valve system.The results obtained by the proposed model are compared with those by the corresponding 1D model and the experimental ones provided by the literature,respectively. The comparison shows that the maximum pressure heads simulated by the proposed model are more accurate than those by the corresponding 1D model.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

HIGH RESOLUTION PARAMETRIC MODELLING FOR TWO-DIMENSIONAL RADAR TARGET USING PRONY ALGORITHM

On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.

Water Vapor,Cloud,and Surface Rainfall Budgets Associated with the Landfall of Typhoon Krosa(2007):A Two-Dimensional Cloud-Resolving Modeling Study

Water vapor, cloud, and surface rainfall budgets associated with the landfall of Typhoon Krosa on 6-8 October 2007 are analyzed based on a two-dimensional cloud-resolving model simulation. The model is integrated with imposed zonally-uniform vertical velocity, zonal wind, horizontal temperature, and vapor advection from NCEP/Global Data Assimilation System （GDAS） data. The simulation data that are validated with observations are examined to study physical causes associated with surface rainfall processes during the landfall. The time- and domain-mean analysis shows that when Krosa approached the eastern coast of China on 6 October, the water vapor convergence over land caused a local atmospheric moistening and a net condensation that further produced surface rainfall and an increase of cloud hydrometeor concentration. Meanwhile, latent heating was balanced by advective cooling and a local atmospheric warming. One day later, the enhancement of net condensation led to an increase of surface rainfall and a local atmospheric drying, while the water vapor convergence weakened as a result of the landfall-induced deprivation of water vapor flux. At the same time, the latent heating is mainly compensated the advective cooling. Further weakening of vapor convergence on 8 October enhanced the local atmospheric drying while the net condensation and associated surface rainfall was maintained. The latent heating is balanced by advective cooling and a local atmospheric cooling.

Modeling of air quality with a modified two-dimensional Eulerian model: A case study in the Pearl River Delta (PRD) region of China

A modified two-dimensional Eulerian air quality model was used to simulate both the gaseous and particulate pollutant concentrations during October 21-24, 2004 in the Pearl River Delta （PRD） region, China. The most significant improvement to the model is the added capability to predict the secondary organic aerosols （SOA） concentrations because of the inclusion of the SOA formation chemistry. The meteorological input data were prepared using the CALMET meteorological model. The concentrations of aerosol-bound species such as NO3^-, NH4^＋, SO4^2-, and SOA were calculated in the fine particle size range （〈2.5 μm）. The results of the two-dimensional model were compared to the measurements at the ground level during the PRD Intensive Monitoring Campaign （IMC）. Overall, there were good agreements between the measured and modeled concentrations of inorganic aerosol components and O3. Both the measured and the modeled results indicated that the maximum hourly O3 concentrations exceeded the China National Air Quality Standard. The predicted 24-h average SOA concentrations were in reasonable agreement with those predicted by the method of minimum OC/EC ratio.

Precision of meshfree methods and application to forward modeling of two-dimensional electromagnetic sources

Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.

Two-Dimensional Mathematical Model of Tidal Current and Sediment for Oujiang Estuary and Wenzhou Bay

A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the Wenzhou Bay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport are presented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with their initial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of the critical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested. The model verification computation shows that, the tidal levels computed from the model are in good agreement with the field data at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measured along the totally 52 synchronously observed verticals distributed over 8 cross sections. The coraputed tidal water throughputs through the Huangda'ao cross section are close to the measured data. And the computed values of bed deformation from Yangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992. The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of the reclamation project on the Wenzhou shoal are predicted by means of this model.

Two-dimensional Numerical Modeling Research on Continent Subduction Dynamics

Continent subduction is one of the hot research problems in geoscience. New models presented here have been set up and two-dimensional numerical modeling research on the possibility of continental subduction has been made with the finite element software, ANSYS, based on documentary evidence and reasonable assumptions that the subduction of oceanic crust has occurred, the subduction of continental crust can take place and the process can be simplified to a discontinuous plane strain theory model. The modeling results show that it is completely possible for continental crust to be subducted to a depth of 120 km under certain circumstances and conditions. At the same time, the simulations of continental subduction under a single dynamical factor have also been made, including the pull force of the subducted oceanic lithosphere, the drag force connected with mantle convection and the push force of the mid-ocean ridge. These experiments show that the drag force connected with mantle convection is critical for continent subduction.

A Vertical Two-Dimensional Model to Simulate TidalHydrodynamics in A Branched Estuary

A vertical (laterally averaged) two-dimensional hydrodynamic model is developed for tides, tidal current, and salinity in a branched estuarine system. The governing equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. An explicit scheme is employed to solve the continuity equations. The momentum and mass balance equations are solved implicitly in the Cartesian coordinate system. The tributaries are governed by the same dynamic equations. A control volume at the junctions is designed to conserve mass and volume transport in the finite difference schemes, based on the physical principle of continuum medium of fluid. Predictions by the developed model are compared with the analytic solutions of steady wind-driven circulatory flow and tidal flow. The model results for the velocities and water surface elevations coincide with analytic results. The model is then applied to the Tanshui River estuarine system. Detailed model calibration and verification have been conducted with measured water surface elevations, tidal current, and salinity distributions. The overall performance of the model is in qualitative agreement with the available field data. The calibrated and verified numerical model has been used to quantify the tidal prism and flushing rate in the Tanshui River-Tahan Stream, Hsintien Stream, and Keelung River.

Effects of Doubled Carbon Dioxide on Rainfall Responses to Large-Scale Forcing:A Two-Dimensional Cloud-Resolving Modeling Study

ABSTRACT Rainfall responses to doubled atmospheric carbon dioxide concentration were investigated through the analysis of two pairs of two-dimensional cloud-resolving model sensitivity experiments. One pair of experiments simulated pre-summer heavy rainfall over southern China around the summer solstice, whereas the other pair of experiments simulated tropical rainfall around the winter solstice. The analysis of the time and model domain mean heat budget revealed that the enhanced local atmospheric warming was associated with doubled carbon dioxide through the weakened infrared radiative cooling during the summer solstice. The weakened mean pre-summer rainfall corresponded to the weakened mean infrared radiative cooling. Doubled carbon dioxide increased the mean tropical atmospheric warming via the enhanced mean latent heat in correspondence with the strengthened mean infrared radiative cooling during the winter solstice. The enhanced mean tropical rainfall was associated with the increased mean latent heat.

Vertical two-dimensional non-hydrostatic pressure model with single layer

The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and unstable, which makes them unsuitable for wider application. In this study, an efficient model with a single layer is developed. Decomposing the pressure into the hydrostatic and dynamic components and integrating the x-momentum equation from the bottom to the free surface can yield a horizontal momentum equation, in which the terms relevant to the dynamic pressure are discretized semi-implicitly. The convective terms in the vertical momentum equation are ignored, and the rest of the equation is approximated with the Keller-box scheme. The velocities expressed as the unknown dynamic pressure are substituted into the continuity equation, resulting in a tri-diagonal linear system solved by the Thomas algorithm. The validation of solitary and sinusoidal waves indicates that the present model can provide comparable results to the models with multiple layers but at much lower computation cost.

Two-dimensional hydrodynamic robust numerical model of soil erosion based on slopes and river basins

Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.

DETERMINATION OF PARAMETERS IN A TWO-DIMENSIONAL MODEL FOR FIXED-BED ADSORBERS

A new method is proposed for the determination of the parameters in a two-dimensionalmodel which characterizes the properties of axial and radial mixing and mass transport in afixed-bed adsorber.Parameter estimation for the model is carried out with methane-air-5A molecularsieve in a bed under the condition of step injection of tracer from a point on the main axis of thebed by the curve fitting method in the time domain.