Empirical Parameter Values Research Articles

Feedback-driven real-time forecasting method for the arrival times of electric vehicles

Affected by weather conditions, traffic conditions, and driver behavior, the arrival characteristics of electric vehicles (EVs) vary significantly from day to day. This study proposes a feedback-driven real-time forecasting approach that combines historical data to improve the forecasting accuracy of arrival times of EVs. For model-based forecasting methods that sample from probability density functions (PDFs), the related parameter values are dynamically optimized. Compared with sampling from PDFs with empirical parameter values, the dynamic optimal parameter values can track the characteristics of EV arrivals by fully using the continuously updated EV feedback. Considering robustness, a historical data-based support vector clustering technology is utilized to obtain the optimization range of optimal parameter values. As a key of this study, the conservativeness of the optimization range is dynamically adjusted with the periodic updates of EV feedback. The experimental results indicate that, by making full utilization of EV feedback, the proposed method can effectively reduce the forecasting errors of EV arrival times.

Chemotaxis-driven stationary and oscillatory patterns in a diffusive HIV-1 model with CTL immune response and general sensitivity.

In this paper, a reaction-diffusion-chemotaxis HIV-1 model with a cytotoxic T lymphocyte (CTL) immune response and general sensitivity is investigated. We first prove the global classical solvability and L∞-boundedness for the considered model in a bounded domain with arbitrary spatial dimensions, which extends the previous existing results. Then, we apply the global existence result to the case with a linear proliferation immune response and an incidence rate. We study the spatiotemporal dynamics about the three types of spatially homogeneous steady states: infection-free steady state S0, CTL-inactivated infection steady state S1, and CTL-activated infection steady state S∗. Our analyses indicate that S0 is globally asymptotically stable if the basic reproduction number R0 is less than 1; if R0 is between 1 and a threshold, then S1 is globally asymptotically stable. However, if R0 is larger than the threshold, then the chemoattraction and chemorepulsion can destabilize S∗, and thus, a spatiotemporal pattern forms as the chemotactic sensitivity crosses certain critical values. We obtain two kinds of important patterns, which are induced by chemotaxis: stationary Turing pattern and irregular oscillatory pattern. We also find that different chemotactic response functions can affect system's dynamics. Based on some empirical parameter values, numerical simulations are given to illustrate the effectiveness of the theoretical predications.

ГЕОДЕЗІЯ, КАРТОГРАФІЯ І АЕРОФОТОЗНІМАННЯ

Nowadays, GIS Technologies are used in many areas of human life, both in everyday life and scientific research. The purpose of the presented study is to identify the relationship between the scale of topographic maps and the main hydrographic characteristics of the river based on the data of observations in the Siverskyi Donets basin. The study is based on the results of the identification of the hydrographic network, which was performed based on the Open Street Map in the GIS environment of the QGIS program using the method of A. N. Straller and I. N. Hartzman. The process of identifying, describing, and analyzing subcontracted connections consists in assigning its identification order to each element of the river network, which makes it possible to compare and standardize streams. Operating with a hierarchical "tree" of the channel network, the main characteristic of which is the number of elementary, unbranched streams, it is possible to identify and analytically describe the dependencies between the detailing of the map and the main characteristics of the river network structures, such as water discharges, network density, drainage basin area, and river length. The basis for describing these relationships was the method of B.V. Kindiuk, who introduced the concept of the coefficient of the river network structure or the fractional order of the stream as a basis for approximating the above-mentioned dependencies, which allows mathematically describing the obtained functions and obtaining numerical values of empirical parameters. Using QGIS made it possible to create maps of the Siverskyi Donets hydrographic network within Ukraine based on maps of scales 1:50 000 and 1:200, 000. With their help, as well as data from a 1:100,000 map, the number of elementary unbranched watercourses was calculated, and each element of the system was identified, where the order of the main river changes depending on the map scale. The change in these indicators shows a tendency to increase the density and complexity of the river network with increasing map detail, and, as a result, potential changes in indicators of the catchment area, water runoff, and river length. The identified dependencies were expressed mathematically in the form of functions, and are also characterized by high values of the approximation reliability coefficient, which made it possible to construct a general transitional graph from the order of the water flow to the scale of the map with the corresponding values of the calculation parameters. The novelty and practical significance lie in the fact that the use of modern GIS technologies in hydrological science significantly increases the quality of cartographic data and concerning the studied object - the Siverskyi Donets River creates a database in the form of digital maps for further use in hydrographic studies. This sub‑basin has not been previously studied using the methodology proposed by B. V. Kindiuk about the influence of map scales on the characteristics of the river network structure. Such study from a practical point of view, can significantly help the work of engineers, researchers, and designers with cartographic data. This study is designed to explain the peculiarities in the scaling of river networks and propose a mechanism for a scientifically based transition from the existing map scale to the desired one within the Siverskyi Donets sub-basin.

Hydroelectricity Modeling for Low‐Carbon and No‐Carbon Grids: Empirical Operational Parameters for Optimization and Dispatch Models

AbstractHydroelectric power has unusual technical characteristics that could become more valuable as the penetration of variable generation renewables grows, but its use for electricity generation is constrained by complex physical, safety, and socioenvironmental considerations. Hydroelectricity can therefore be difficult to represent in national‐scale energy models, and is frequently presumed to be either overly flexible or inflexible. While a few grid models address this complexity via detailed hydraulic process models, more simplified optimization and dispatch models could benefit from the use of empirical parameter values. In this study, we combine a new data set comprising 7.8 million flow‐hours of data from 2011 to 2016 at 158 dams across the United States with monthly and hourly generation data from the U.S. Energy Information Administration (EIA) to elucidate such empirical constraints for the continental United States. We introduce an approach for estimating power generation from hourly water discharge, then present regionally resolved interannual seasonal and diurnal generation patterns; frequency analyses of ramp rates; minimum and maximum generation rates; daily reversals; and load duration curves, all available interactively with a new data visualization tool. We suggest that due largely to hydropower's role as power generation that also serves non‐energy purposes, it acts as a predictable variable generator with constrained dispatchability, more like a supply side analog to demand response resources than like a battery. This observation is particularly relevant for high penetration renewable energy scenarios, given hydroelectric generators' expected value for grid stabilization and load balancing.

Toward a Global Model for Soil Inorganic Phosphorus Dynamics: Dependence of Exchange Kinetics and Soil Bioavailability on Soil Physicochemical Properties.

The representation of phosphorus (P) cycling in global land models remains quite simplistic, particularly on soil inorganic phosphorus. For example, sorption and desorption remain unresolved and their dependence on soil physical and chemical properties is ignored. Empirical parameter values are usually based on expert knowledge or data from few sites with debatable global representativeness in most global land models. To overcome these issues, we compiled from data of inorganic soil P fractions and calculated the fraction of added P remaining in soil solution over time of 147 soil samples to optimize three parameters in a model of soil inorganic P dynamics. The calibrated model performed well (r 2 > 0.7 for 122 soil samples). Model parameters vary by several orders of magnitude, and correlate with soil P fractions of different inorganic pools, soil organic carbon and oxalate extractable metal oxide concentrations among the soil samples. The modeled bioavailability of soil P depends on, not only, the desorption rates of labile and sorbed pool, inorganic phosphorus fractions, the slope of P sorbed against solution P concentration, but also on the ability of biological uptake to deplete solution P concentration and the time scale. The model together with the empirical relationships of model parameters on soil properties can be used to quantify bioavailability of soil inorganic P on various timescale especially when coupled within global land models.

Financial wellbeing of households in instability

In instability and economic turbulence, the wellbeing of households as market economy entities constitutes the financial-investment capacity of a region, the level of which is determined by the conditions of the competitive socio-economic environment. The paper aims to estimate the financial wellbeing of households on the example of the oblasts of the Carpathian region of Ukraine in instability. The study is based on a system-integral estimation method, which includes the implementation of three stages: (1) development of a system of indicators, (2) determination and substantiation of weight significance, and (3) construction of time series of empirical parameters of households’ wellbeing based on temporal and spatial approaches. The analysis reveals that the financial wellbeing of households differentiates in a competitive economic environment and with the spread of behavioral factors (COVID-19, consumer reflections). Among the oblasts of the Carpathian region of Ukraine, the highest values of empirical parameters of financial wellbeing were in Zakarpatska (0.537) and Chernivetska (0.459) oblasts (2019). Meanwhile, the level of the financial wellbeing of households is higher in Lvivska oblast by several indicators. The divergence of the Carpathian region from Ukraine by the level of the financial wellbeing of households was mostly observed in 2018–2019. Zakarpatska oblast was the leader by the level of the financial wellbeing of households in 2010–2019. The study is of the practical nature for framing the regional economic policy in terms of detecting the critical “pressure” of financial wellbeing on the economic growth of the region and economic ability to increase investment capacity. AcknowledgmentsThe study has been conducted within the framework of the Applied Research “Financial determinants of the provision of economic growth in the regions and territorial communities based on behavioural economy” with the support of the National Research Foundation of Ukraine (M. Dolishniy Institute of Regional Research of the National Academy of Sciences of Ukraine, the grant Reg. No. 2020.02/0215, 2020–2022).

Construction of a photosynthetic rate prediction model for greenhouse strawberries with distributed regulation of light environment.

In winter and spring, for greenhouses with larger areas and stereoscopic cultivation, distributed light environment regulation based on photosynthetic rate prediction model can better ensure good crop growth. In this paper, strawberries at flowering-fruit stage were used as the test crop, and the LI-6800 portable photosynthesis system was used to control the leaf chamber environment and obtain sample data by nested photosynthetic rate combination experiments under temperature, light and CO2 concentration conditions to study the photosynthetic rate prediction model construction method. For a small-sample, nonlinear real experimental data set validated by grey relational analysis, a photosynthetic rate prediction model was developed based on Support vector regression (SVR), and the particle swarm algorithm (PSO) was used to search the influence of the empirical values of parameters, such as the penalty parameter C, accuracy ε and kernel constant g, on the model prediction performance. The modeling and prediction results show that the PSO-SVR method outperforms the commonly used algorithms such as MLR, BP, SVR and RF in terms of prediction performance and generalization on a small sample data set. The research in this paper achieves accurate prediction of photosynthetic rate of strawberry and lays the foundation for subsequent distributed regulation of greenhouse strawberry light environment.

On the identification and mapping of three distinct stages of single vapor bubble growth with the corresponding microlayer dynamics

An experimental study to establish a one-on-one correspondence between microlayer dynamics and growth process of single vapor bubble under nucleate pool boiling regime has been presented. The transient evolution of microlayer and the corresponding vapor bubble growth have been simultaneously recorded using a thin film interferometer and one of the gradients-based imaging techniques, namely rainbow schlieren deflectometry. Based on interferometry measurements, the microlayer is observed to undergo three distinct stages. Each of these stages is identified by distinct microlayer structure namely, inverted conical, flattened peripheral region and truncated conical structure during the initial, transitional and diffusion-controlled growth phases of the vapor bubble. Contribution of microlayer towards bubble growth has been established by quantifying the time rate of change of net mass of the superheated fluid entrapped in the microlayer region and the corresponding rate of change of mass of the growing bubble. The analysis showed that the ratio of mass of the superheated fluid in the microlayer region to the mass of the vapor bubble does not exceed a maximum of ~ 15%. For any given heat flux level, the growing vapor bubble is seen to acquire its maximum size in the transitional regime of the growth process during which the microlayer exhibits a flattened peripheral region. Possible impact of varying heat flux level on microlayer dynamics and movement of the apparent contact line of the vapor bubble has also been discussed. A reasonably good agreement of the experimental data with some of the well-established bubble growth models for initial microlayer thickness has been observed. Based on the experimental results, critical value of empirical parameter ‘Cslope’ (parameter correlating the bubble base radius and the maximum microlayer thickness) for transition from inertia-driven to diffusion dominated growth has been identified to be ~ 0.0042.

Comparisons of Two Approaches for Geotechnical Model Calibration with Scarce Data

Geotechnical models are usually built upon assumptions and simplifications, inevitably resulting in discrepancies between model predictions and measurements. To enhance prediction accuracy, geotechnical models are typically calibrated against measurements by bringing in additional empirical or semiempirical correction terms. Different approaches have been used in the literature to determine the optimal values of empirical parameters in the correction terms. When measured data are abundant, calibration outcomes using different approaches can be expected to be practically the same. However, if measurements are scarce or limited, calibration outcomes could differ significantly, depending largely on the adopted calibration approach. In this study, we examine two most commonly used approaches for geotechnical model calibration in the literature, namely, (1) purely data‐catering (PDC) approach, and (2) root mean squared error (RMSE) method. Here, the purely data‐catering approach refers to selection of empirical parameter values that minimize coefficient of variation of model factor while maintains its mean value of one, based solely on measured data. A real case of calibrating the Federal Highway Administration (FHWA) simplified facing load model for design of soil nail walls is illustrated to thoroughly elaborate the differences in practical calibration and design outcomes using the two approaches under scarce data conditions.

Study on Evaluation Model of Chinese P2P Online Lending Platform Based on Hybrid Kernel Support Vector Machine

Accurate evaluation of the risk level and operation performances of P2P online lending platforms is not only conducive to better functioning of information intermediaries but also effective protection of investors’ interests. This paper proposes a genetic algorithm (GA) improved hybrid kernel support vector machine (SVM) with an index system to construct such an evaluation model. A hybrid kernel consisting of polynomial function and radial basis function is improved, specifically kernel parameters and the weight of two kernels, by GA method with excellent global optimization and rapid convergence. Empirical testing based on cross-sectional data from Chinese P2P lending market demonstrates the superiority of the improved hybrid kernel SVM model. The classification accuracy of credit risk level and operation quality is higher than the single kernel SVM model as well as the hybrid kernel model with empirical parameter values.

Compartmental approach for modelling twin-screw granulation using population balances.

In this study, a compartmental population balance model (CPBM) is developed as a predictive tool of particle size distribution (PSD) for wet granulation in co-rotating twin-screw granulator (TSG). This model is derived in terms of liquid to solid ratio (L/S) and screw speed representing the main process parameters of the TSG. The mathematical model accounts for aggregation and breakage of the particles occurring in five compartments of the TSG with inhomogeneous screw configurations (3 conveying zones and 2 kneading zones). Kapur's aggregation kernel is implemented in granulation and finite volume numerical method is adapted for solving the mathematical model. The results show a dramatic improvement in solution accuracy compared to the cell average numerical method. Moreover, Kriging interpolation is used to interpolate for new values of empirical parameters at different L/S and screw speeds. Finally, the CPBM model is calibrated and validated using the experimental data.

New approaches for the determination of electrochemical parameters in the model of proton-conducting solid oxide fuel cell

Proton-conducting solid oxide fuel cell is an electrical energy generation device with great potentials in the industrial and civil applications. The values of empirical parameters adopted in the corresponding simulations significantly determines the quality and reliability of the modeling results. Two new and practicable approaches are proposed in this paper, in order to accurately evaluate several critical electrochemical parameters in the relevant simulations. According to our models, effective hole conductivity of the electrolyte under working conditions, along with cathode electrochemical reaction order, are directly derived from the experiment. To examine the proposed approaches, single cells with BaCe0.8Sm0.2O3-δ electrolyte and Sm0.5Sr0.5CoO3-δ–BaCe0.8Sm0.2O3-δ cathode are manufactured, further the related experiments are performed as case studies. To determine the cathode activation overpotentials for cells with different particle sizes, a transferable parameter named TPB-specific exchange current density is introduced, and an expression which is much more convenient for numerical simulation of the cathode activation overpotential than the previously reported one is put forward in our model. These proposed expressions and parameters possess desirable transferability, and accurate numerical predictions of the cell output performances are achieved for both our and reported experiments. Then the influences of the involved polarizations are explored. Emphasis is drawn on the electrolyte hole conductivity. Its impact on the cell performance is non-negligible particularly under relatively high temperature and low output current, and hence has to be considered seriously in the simulation.

Tephra Fallout Hazard Assessment for VEI5 Plinian Eruption at Kuju Volcano, Japan, Using TEPHRA2

Tephra fallout has a potential impact on engineered structures and systems at nuclear power plants. We provide the first report estimating potential accumulations of tephra fallout as big as VEI5 eruption from Kuju Volcano and calculated hazard curves at the Ikata Power Plant, using the TEPHRA2 computer program. We reconstructed the eruptive parameters of Kj-P1 tephra fallout deposit based on geological survey and literature review. A series of parameter studies were carried out to determine the best values of empirical parameters, such as diffusion coefficient and the fall time threshold. Based on such a reconstruction, we represent probabilistic analyses which assess the variation in meteorological condition, using wind profiles extracted from a 22 year long wind dataset. The obtained hazard curves and probability maps of tephra fallout associated to a Plinian eruption were used to discuss the exceeding probability at the site and the implications of such a severe eruption scenario.

Refractive index, speed of sound, FT-IR and computational study of intermolecular interaction between binary mixtures of 1,8-cineole with o-, m- and p-cresol at 303.15, 308.15 and 313.15 K

Refractive indexes (nD) and speed of sound (u) have been experimentally measured for binary mixtures of 1,8-cineole+o-, m- and p-cresol at entire composition range at 303.15, 308.15 and 313.15K at atmospheric pressure. The deviation in refractive index (∆nD), molar refraction (Rm) and deviation in molar refraction (∆Rm) has been calculated using experimentally measured values of the refractive index. Similarly, from the measured values of speed of sound (u), some of its related properties like deviation in speed of sound (∆u), isentropic compressibility (ks), deviation in isentropic compressibility (∆ks), acoustic impedance (z), deviation in acoustic impedance (∆z) have been calculated. The partial molar isentropic compressibility K−s,m,io, excess partial molar isentropic compressibility K−s,m,io,E, apparent molar isentropic compressibility (Ks,m,Φ,i) have been also calculated. The values of apparent molar isentropic compressibility (Ks,m,Φ,i) have been fitted with the Redlich-Rosenberg-Mayer equation to get limiting apparent molar isentropic compressibility (Kos,m,Φ,i) and Sk and Bk empirical parameter values. All these calculated results were fitted with the Redlich-Kister polynomial equation. The FT-IR spectra of pure and their binary mixtures are also studied. The geometrically optimization, harmonic vibrational wave numbers, bond characteristics and single point energy has been theoretically calculated using density functional theory (DFT-B3LYP) method with 6-31+G(d,p) basis sets for pure components and their binary mixtures. Positive values of ∆nD, ∆u, ∆z and negative values of ∆Rm, ∆ks and K−s,m,io,E support for the presence of strong intermolecular interactions between components of these binary mixtures.

Field evidence of beach profile evolution toward equilibrium

AbstractAn equilibrium framework is used to describe the evolution of the cross‐shore profile of five beaches (medium grain size sand) in southern California. Elevations were observed quarterly on cross‐shore transects extending from the back beach to 8 m depth, for 3–10 years. Transects spaced 100 m in the alongshore direction are alongshore averaged into nineteen 700–900 m long sections. Consistent with previous observations, changes about the time average profile in many sections are captured by the first mode empirical orthogonal function (EOF). The first EOF poorly describes sections with hard substrate (less than roughly 80% sandy bottom) and also fails near the head of a submarine canyon and adjacent to an inlet. At the 12 well‐described sections, the time‐varying amplitude of the first EOF, the beach state A, describes the well‐known seasonal sand exchange between the shoreline and offshore (roughly between 4 and 7 m depth). We show that the beach state change rate dA/dt depends on the disequilibrium between the present state A and wave conditions, consistent with the equilibrium concepts of Wright and Short (1984) and Wright et al. (1985). Empirically determined, optimal model coefficients using the framework of Yates et al. (2009a, 2011) vary between sections, but a single set of globally optimized values performs almost as well. The model implements equilibrium concepts using ad hoc assumptions and empirical parameter values. The similarity with observed profile change at five southern California beaches supports the underlying model equilibrium hypotheses, but for unknown reasons the model fails at Duck, NC.

Air density effects on aeolian sand movement: Implications for sediment transport and sand control in regions with extreme altitudes or temperatures

AbstractAeolian sand transport results from interactions between the ground surface and airflow. Previous research has focused on the effects on sand entrainment and mass transport of surface features and wind velocity, but the influence of air density, which strongly constrains airflow characteristics and the resulting sand flow, has not been widely considered. In the present study, entrainment, saltation characteristics and transport rates were examined at nine experimental sites ranging in elevation from −154 m below sea‐level (Aiding Lake) to 5076 m above sea‐level (Tanggula Mountain pass on the Qinghai–Tibetan plateau). At each site, a portable wind tunnel and high‐speed camera system were set up, and the friction wind velocity, threshold friction velocity and sand flow structure were observed systematically. For a given volumetric airflow, lower air density increases the wind velocity. Low air density also creates a high threshold friction velocity. The Bagnold wind erosion threshold model remains valid, but the value of empirical parameter A decreased with decreasing air density and ranged from 0·10 to 0·07, the smallest values reported in the literature. For a given wind velocity, increased altitude reduced total sand transport and creeping, but the saltation rate and saltation height increased. The present results provide insights into the fundamental mechanisms of the initiation and transport of sand by wind in regions with an extreme temperature or altitude (for example, alpine deserts and low‐lying lake basins) or on other planets, including Mars. These results also provide theoretical support for improved sand‐control engineering measures. The data and empirical equations provided in this paper improve the ability to estimate threshold and transport conditions for wind‐blown sand.

New Results for an Approximate Method for Calculating Two-Dimensional Furrow Infiltration

AbstractThis study expands the analysis of a proposed furrow infiltration formulation based on an approximate solution to the two-dimensional Richards equation. The approach calculates two-dimensional infiltration flux as the sum of one-dimensional infiltration and a second term labeled the edge effect. The edge effect varies linearly with time when the applied water pressure is constant. It is a function of sorptivity, soil water content, and the empirical parameters γ and W* that require calibration for the specific soil and furrow geometry. The primary objectives of the analysis were to better understand how furrow geometry affects the edge effect and the resulting empirical parameter values and to provide additional guidance for calibration. For a given flow depth, furrow geometry defines the wetted perimeter and the water pressure applied along that wetted perimeter. A subobjective was to evaluate the effect of the assumed soil water retention model on the calibration results. The results show that t...

Determination of correlation parameters for evaluation of mechanical properties by Small Punch Test and Automated Ball Indentation Test for Zr–2.5% Nb pressure tube material

Empirical correlation parameters for evaluation of key mechanical properties of anisotropic Zr–2.5% Nb pressure tube (PT) material have been determined for Small Punch Test (SPT) and Automated Ball Indentation (ABI) Test. Parameters have been developed as a function of loading direction due to highly anisotropic nature of the material. Slopes of the linear fit correlations established for the SPT and ABI tests, showed good sensitivity of these techniques to wide range of microstructures in the alloy that resulted due to different heat treatment conditions. In the case of ABI, the difference in yield strength offset constant between the two directions of loading studied, represented the effect of texture. Similarly in the case of SPT the difference in intercept values of the, linear fit correlations established for evaluation of yield strength (YS) and the ultimate tensile strength (UTS) with respect to two directions of loading studied, clearly indicated the effect of anisotropy in the material. Comparison of tensile strength data, the ABI and SPT derived σy and σUTS data obtained for heat treated specimens used in the study, through currently established correlation parameters showed good agreement with each other for both the directions of loading, indicating proper standardisation of ABI and SPT techniques and proper values of empirical parameters have been arrived at for evaluation of key mechanical properties of Zr–2.5% Nb material. The overall study established that, the direction of loading bears significance while dealing with anisotropic material like Zr–2.5% Nb alloy and hence, direction specific correlations need to be used for evaluation of true mechanical properties for a given miniature test.

Impact of capillary pressure on permeability impairment during CO2 injection into deep saline aquifers

During CO2 injection in deep saline aquifers, salt precipitation happens around the injection well because of capillary driven back flow, inducing permeability impairment. The permeability impairment affects CO2 injectivity and migration. Different values of three characteristic parameters for capillary pressure function (air entry pressure, empirical parameter m and liquid residual saturation) as well as input absolute value of maximum capillary are chosen in numerical simulation to figure out their effects on salt precipitation. Verma & Pruess model is then used for quantifying permeability impairment. Results show that permeability decreases with higher air entry pressure, larger liquid residual saturation, and especially smaller value of empirical parameter m. To enhance CO2 injectivity and avoid blocking of CO2 migration, a homogenous formation with large pore size should be chosen before CO2 injection into deep saline aquifer.

Optimal Nonlinear Control in MWD assisted Directional Drilling Process

In the present work an algorithm is proposed for automatic guidance of bit in directional drilling process. In the manual mode control of bit's movement some problems arise due to operational and human inaccuracies which make the actual path differ from the desired one. In this article an automatic control system based on the bang-bang control strategy equipped with MWD (Measuring While Drilling) has been designed to lead the bit to the target. The bang-bang control strategy is one of the optimal nonlinear control methods. According to the kinematic constraints and system's input variables of the present problem this control strategy is suitable and simple to implement. In the first step the kinematic model of the drillstring movement has been derived. In the second step using Matlab-Simulink software, the system's kinematic equations besides some effects such as mechanical behavior of the well, dead zone between the formation and the bit, changing the geological formation during the drilling, time delay in sending and receiving mechanical and electrical signals, and environmental noises have been taken into account and are modeled. Directional drilling has been simulated with empirical parameter values of the well formation in Nargesi field. To study the control system abilities and limitations, results were obtained and studied. It is shown that the control system can generate an appropriate path when the result is compared with an actual drilling data. In the final part the concluding remarks are presented.