Patient-specific closed-loop model of the fontan circulation: Calibration and validation.

The Fontan circulation, designed for managing patients with a single functional ventricle, presents challenges in long-term outcomes. Computational methods offer potential solutions, yet their application in cardiology practice remains largely unexplored. Our aim was to assess the ability of a patient-specific, closed-loop, reduced-order blood flow model to simulate pulsatile blood flow in the Fontan circulation. Using one-dimensional models, we simulated the aorta, superior and inferior venae cavae, and right and left pulmonary arteries, while lumping heart chambers and remaining vessels into zero-dimensional models. The model was calibrated with patient-specific haemodynamic data from combined cardiac catheterisation and magnetic resonance exams, using a novel physics-based stepwise methodology involving simpler open-loop models. Testing on a 10-year-old, anesthetised patient, demonstrated the model's capability to replicate pulsatile pressure and flow in the larger vessels and ventricular pressure. Average relative errors in mean pressure and flow were 2.9 % and 3.6 %, with average relative point-to-point errors (RPPE) in pressure and flow at 5.2 % and 16.0 %. Comparing simulation results to measurements, mean aortic pressure and flow values were 50.7 vs. 50.4 mmHg and 41.6 vs. 41.9 ml/s, respectively, while ventricular pressure values were 28.7 vs. 27.4 mmHg. The model accurately described time-varying ventricular volume with a RPPE of 2.9 %, with mean, minimum, and maximum ventricular volume values for simulation results vs. measurements at 59.2 vs. 58.2 ml, 38.0 vs. 37.6 ml, and 76.0 vs. 74.4 ml, respectively. It provided physiologically realistic predictions of haemodynamic changes from pulmonary vasodilation and atrial fenestration opening. The new model and calibration methodology are freely available, offering a platform to virtually investigate the Fontan circulation's response to clinical interventions and explore potential mechanisms of Fontan failure. Future efforts will concentrate on broadening the model's applicability to a wider range of patient populations and clinical scenarios, as well as testing its effectiveness.

A numerical study on the siphonic effect of enhanced external counterpulsation at lower extremities with a coupled 0D-1D closed-loop personalized hemodynamics model.

Enhanced external counterpulsation (EECP) is a treatment and rehabilitation approach for ischemic diseases, including coronary artery disease. Its therapeutic benefits are primarily attributed to the improved blood circulation achieved through sequential mechanical compression of the lower extremities. However, despite the crucial role that hemodynamic effects in the lower extremity arteries play in determining the effectiveness of EECP treatment, most studies have focused on the diastole phase and ignored the systolic phase. In the present study, a novel siphon model (SM) was developed to investigate the interdependence of several hemodynamic parameters, including pulse wave velocity, femoral flow rate, the operation pressure of cuffs, and the mean blood flow changes in the femoral artery throughout EECP therapy. To verify the accuracy of the SM, we coupled the predicted afterload in the lower extremity arteries during deflation using SM with the 0D-1D patient-specific model. Finally, the simulation results were compared with clinical measurements obtained during EECP therapy to verify the applicability and accuracy of the SM, as well as the coupling method. The precision and reliability of the previously developed personalized approach were further affirmed in this study. The average waveform similarity coefficient between the simulation results and the clinical measurements during the rest state exceeded 90%. This work has the potential to enhance our understanding of the hemodynamic mechanisms involved in EECP treatment and provide valuable insights for clinical decision-making.

Tube law parametrization using in vitro data for one-dimensional blood flow in arteries and veins: TUBE LAW PARAMETRIZATION IN ARTERIES AND VEINS.

The deformability of blood vessels in one-dimensional blood flow models is typically described through a pressure-area relation, known as the tube law. The most used tube laws take into account the elastic and viscous components of the tension of the vessel wall. Accurately parametrizing the tube laws is vital for replicating pressure and flow wave propagation phenomena. Here, we present a novel mathematical-property-preserving approach for the estimation of the parameters of the elastic and viscoelastic tube laws. Our goal was to estimate the parameters by using ovine and human in vitro data, while constraining them to meet prescribed mathematical properties. Results show that both elastic and viscoelastic tube laws accurately describe experimental pressure-area data concerning both quantitative and qualitative aspects. Additionally, the viscoelastic tube law can provide a qualitative explanation for the observed hysteresis cycles. The two models were evaluated using two approaches: (i) allowing all parameters to freely vary within their respective ranges and (ii) fixing some of the parameters. The former approach was found to be the most suitable for reproducing pressure-area curves.

A 1-D model to retrieve the vertical profiles of minor atmospheric constituents for cloud microphysical modeling: II. Simulation of diurnal cycle.

The accurate forecast of the diurnal cycle of the number concentration of trace gases is vital due to their influence on precipitation processes by controlling the number concentration of cloud condensation nuclei (CCN). 1-D hybrid Monte Carlo-Gear solver developed to retrieve vertical profiles of the number concentration of CCNs for microphysics modeling has been tested for representation of the diurnal cycle in the present paper. The retrieved profiles of CH4 and SO2 have been tested with the Copernicus Atmosphere Monitoring Service (CAMS) model at 3-hour time intervals for four megacities: Delhi, Kolkata, Chennai, and Mumbai for rainy and non-rainy days. The retrieved profiles have shown diurnal variation up to 18 UTC at all pressure levels with lead or lag with that of the CAMS model. After 18 UTC there was a furious increase in the number concentrations. During non-rainy days, the 1-D model slightly overestimated (underestimated) the maximum (minimum) number concentrations of CH4 over Delhi whereas concentrations are overestimated over Kolkata, Chennai, and Mumbai. Forecasted CH4 has a good (weak) correlation over Chennai (Mumbai) respectively. The 1-D model overestimated (overestimated) the maximum (minimum) number concentrations of SO2 over Delhi but the maximum (minimum) concentrations are underestimated (overestimated) over Kolkata, Chennai, and Mumbai. The number concentrations of SO2 have shown a good correlation for all megacities except Delhi. CH4 number concentration is overestimated during rainy days. Delhi and Kolkata show a good correlation of CH4 during rainy days. SO2 during rainy days is underestimated except over Chennai and both models show a good correlation except over Mumbai. Overall, it can be stated that the 1-D hybrid solver is successful in simulating the monthly mean diurnal variation of vertical profiles of CH4 and SO2, and its implementation in the global model may estimate the number concentrations with better accuracies.

Evaluation and Correlation Analysis of the Rheological Properties of Ground Tire Rubber and Styrene Butadiene Styrene Compound-Modified Asphalt.

With the increase in highway traffic volume, many waste tires are being produced, which puts serious pressure on the global ecological environment. Processing waste tires into powder and adding them to asphalt is an important and effective way to solve this noticeable environmental challenge. In this paper, to produce ground tire rubber (GTR) and styrene-butadiene-styrene (SBS) compound-modified asphalt, GTR was put into SBS-modified asphalt (GTRSA). Subsequently, some ordinary property tests, frequency sweep tests, and multiple stress creep recovery tests were conducted to investigate the conventional properties and rheological properties of GTRSA. Moreover, the 2S2P1D (two springs, two parabolic elements, and one dashpot) model was adopted to analyze the consequences of adding GTR content on the rheological properties of GTRSA. Finally, the Pearson correlation coefficient was employed to reveal the connection between the conventional properties and the rheological properties. The results show that GTR has a great impact on improving the rutting resistance, thermo-sensitive performance, shear resistance capability, stress sensitivity, and creep recovery performance of GTRSA. Adding 20% GTR can improve the creep recovery rate to 80.8%. The 5 °C ductility index suggests that GTR makes a difference to the low-temperature properties. The rheological properties and conventional properties had a strong linear link.

Parameter Determination of the 2S2P1D Model and Havriliak-Negami Model Based on the Genetic Algorithm and Levenberg-Marquardt Optimization Algorithm.

This study utilizes the genetic algorithm (GA) and Levenberg-Marquardt (L-M) algorithm to optimize the parameter acquisition process for two commonly used viscoelastic models: 2S2P1D and Havriliak-Negami (H-N). The effects of the various combinations of the optimization algorithms on the accuracy of the parameter acquisition in these two constitutive equations are investigated. Furthermore, the applicability of the GA among different viscoelastic constitutive models is analyzed and summarized. The results indicate that the GA can ensure a correlation coefficient of 0.99 between the fitting result and the experimental data of the 2S2P1D model parameters, and it is further proved that the fitting accuracy can be achieved through the secondary optimization via the L-M algorithm. Since the H-N model involves fractional power functions, high-precision fitting by directly fitting the parameters to experimental data is challenging. This study proposes an improved semi-analytical method that first fits the Cole-Cole curve of the H-N model, followed by optimizing the parameters of the H-N model using the GA. The correlation coefficient of the fitting result can be improved to over 0.98. This study also reveals a close relationship between the optimization of the H-N model and the discreteness and overlap of experimental data, which may be attributed to the inclusion of fractional power functions in the H-N model.

One-dimensional analysis method of pulsatile blood flow in arterial network for REBOA operations.

Based on the generalized Darcy model, here we develop a linear one-dimensional (1D) composite model to predict the effects of the inserted balloon under REBOA operations on the dynamic characteristics of blood flow in flexible arterial networks. We first consider the effect of the decrease of cardiac output under different degrees of blood loss through employing the fourth-order lumped parameter model of cardiovascular system. Then, the effect of the inserted balloon is included by developing the relation between flow resistance and occlusion ratio with the neural network approach. Finally, the accuracy of the developed 1D composite model for REBOA operations, which can be solved analytically in the frequency domain, is verified by comparing to computational fluid dynamics (CFD) simulations. It is demonstrated that the 1D model is able to reproduce main features of the systemic circulation under balloon occlusion of the aorta during REBOA surgery. The 1D composite model could substantially reduce the computational time, which makes it possible to give the instant prediction of the working parameters during RABOA operations.

Tumor Evolution Models of Phase-Field Type with Nonlocal Effects and Angiogenesis.

In this survey article, a variety of systems modeling tumor growth are discussed. In accordance with the hallmarks of cancer, the described models incorporate the primary characteristics of cancer evolution. Specifically, we focus on diffusive interface models and follow the phase-field approach that describes the tumor as a collection of cells. Such systems are based on a multiphase approach that employs constitutive laws and balance laws for individual constituents. In mathematical oncology, numerous biological phenomena are involved, including temporal and spatial nonlocal effects, complex nonlinearities, stochasticity, and mixed-dimensional couplings. Using the models, for instance, we can express angiogenesis and cell-to-matrix adhesion effects. Finally, we offer some methods for numerically approximating the models and show simulations of the tumor's evolution in response to various biological effects.

A multiscale computational framework to evaluate flow alterations during mechanical thrombectomy for treatment of ischaemic stroke.

The treatment of ischaemic stroke increasingly relies upon endovascular procedures known as mechanical thrombectomy (MT), which consists in capturing and removing the clot with a catheter-guided stent while at the same time applying external aspiration with the aim of reducing haemodynamic loads during retrieval. However, uniform consensus on procedural parameters such as the use of balloon guide catheters (BGC) to provide proximal flow control, or the position of the aspiration catheter is still lacking. Ultimately the decision is left to the clinician performing the operation, and it is difficult to predict how these treatment options might influence clinical outcome. In this study we present a multiscale computational framework to simulate MT procedures. The developed framework can provide quantitative assessment of clinically relevant quantities such as flow in the retrieval path and can be used to find the optimal procedural parameters that are most likely to result in a favorable clinical outcome. The results show the advantage of using BGC during MT and indicate small differences between positioning the aspiration catheter in proximal or distal locations. The framework has significant potential for future expansions and applications to other surgical treatments.

Multimedia fate of sulfamethoxazole (SMX) in a water-scarce city by coupling fugacity model and HYDRUS-1D model.

Increased concentrations of pharmaceutical and personal care products (PPCPs) have raised concerns about their impact on the ecological system and human health. To understand the environmental impact of PPCPs, we evaluated the fate of a typical PPCP of sulfamethoxazole (SMX) in a water-scarce city of Tianjin during 2013-2020 using a coupled model based on the dynamic fugacity model and HYDRUS-1D model. The results showed that the coupled model successfully simulated the reported SMX concentrations in the main fate media of water and soils, which accounted for 46.4 % and 53.0 % with equilibrium concentrations of 135-165 ng/L and 0.4-0.5 ng/g, respectively. The cross-media transfer flux results showed that advection was the prime input path for SMX in water, while degradation was the dominant output path. Wastewater irrigation and degradation were the main transfer processes of SMX in the soil. Moreover, human activities (i.e., emission loads) and climate (i.e., temperature and precipitation) change can significantly affect the concentrations and transfer rate of SMX in the media. These findings provide basic data and methods for the risk assessment of SMX in water-scarce regions.

A 1-D model to retrieve the vertical profiles of minor atmospheric constituents for cloud microphysical modeling: I. Formulation and validation.

Determining the number concentration of minor constituents in the atmosphere is very important as it determines the whole tropospheric chemistry process. These constituents may act as cloud condensation nuclei (CCN) and ice nuclei (IN), impacting heterogeneous nucleation inside the cloud. However, the estimations of the number concentration of CCN/IN in cloud microphysical parameters are associated with uncertainties. In the present work, a hybrid Monte Carlo Gear solver has been developed to retrieve profiles of CH4, N2O, and SO2. The idealized experiments have been carried out using this solver for retrieving vertical profiles of these constituents over four megacities, viz., Delhi, Mumbai, Chennai, and Kolkata. Community Long-term Infrared Microwave Coupled Atmospheric Product System (CLIMCAPS) dataset around 0800 UTC (2000UTC) has been used for initializing the number concentration of CH4, N2O, and SO2 for daytime (nighttime). The daytime (nighttime) retrieved profiles have been validated using 2000 UTC (next day 0800 UTC) CLIMCAPS products. ERA5 temperature dataset has been used to estimate the kinematic rate of reactions with 1000 perturbations determined using Maximum Likelihood Estimation (MLE). The retrieved profiles and CLIMCAPS products are in very good agreement, as evidenced by the percentage difference between them within the range of 1.3 × 10-5-60.8 % and the coefficient of determination mainly within the range between 81 and 97 %. However, during the passage of tropical cyclone and western disturbance, its value became as low as 27 and 65 % over Chennai and Kolkata, respectively. The enactment of synoptic scale systems such as western disturbances, tropical cyclone Amphan, and easterly waves caused disturbed weather over these megacities-the retrieved profiles during disturbed weather cause large deviations of vertical profiles of N2O. However, the profiles of CH4 and SO2 have less deviation. It is inferred that incorporating this methodology in the dynamical model will be useful to simulate the realistic vertical profiles of the minor constituents in the atmosphere.

Assessing groundwater recharge rates, water quality changes, and agricultural impacts of large-scale water recycling.

The over-exploitation and insufficient replenishment of groundwater (GW) have resulted in a pressing need to conserve freshwater and reuse of treated wastewater. To address this issue, the Government of Karnataka launched a large-scale recycling (440 million liters/day) scheme to indirectly recharge GW using secondary treated municipal wastewater (STW) in drought-prone areas of Kolar district in southern India. This recycling employs soil aquifer treatment (SAT) technology, which involves filling surface run-off tanks with STW that intentionally infiltrate and recharge aquifers. This study quantifies the impact of STW recycling on GW recharge rates, levels, and quality in the crystalline aquifers of peninsular India. The study area is characterized by hard rock aquifers with fractured gneiss, granites, schists, and highly fractured weathered rocks. The agricultural impacts of the improved GW table are also quantified by comparing areas receiving STW to those not receiving it, and changes before and after STW recycling were measured. The AMBHAS_1D model was used to estimate the recharge rates and showed a tenfold increase in daily recharge rates, resulting in a significant increase in the GW levels. The results indicate that the surface water in the rejuvenated tanks meets the country's stringent water discharge standards for STW. The GW levels of the studied boreholes increased by 58-73 %, and the GW quality improved significantly, turning hard water into soft water. Land use land cover studies confirmed an increase in the number of water bodies, trees, and cultivated land. The availability of GW significantly improved agricultural productivity (11-42 %), milk productivity (33 %), and fish productivity (341 %). The study's outcomes are expected to serve as a role model for the rest of Indian metro cities and demonstrate the potential of reusing STW to achieve a circular economy and a water-resilient system.

A predictive patient-specific computational model of coronary artery bypass grafts for potential use by cardiac surgeons to guide selection of graft configurations.

Cardiac surgeons face a significant degree of uncertainty when deciding upon coronary artery bypass graft configurations for patients with coronary artery disease. This leads to significant variation in preferred configuration between different surgeons for a particular patient. Additionally, for the majority of cases, there is no consensus regarding the optimal grafting strategy. This situation results in the tendency for individual surgeons to opt for a "one size fits all" approach and use the same grafting configuration for the majority of their patients neglecting the patient-specific nature of the diseased coronary circulation. Quantitative metrics to assess the adequacy of coronary bypass graft flows have recently been advocated for routine intraoperative use by cardiac surgeons. In this work, a novel patient-specific 1D-0D computational model called "COMCAB" is developed to provide the predictive haemodynamic parameters of functional graft performance that can aid surgeons to avoid configurations with grafts that have poor flow and thus poor patency. This model has significant potential for future expanded applications.

Predicting Dynamic Properties of Asphalt Mastic Considering Asphalt-Filler Interaction Based on 2S2P1D Model.

The relationship between the various phases of asphalt materials, from asphalt binder to mastic and mixture, has received great attention over the years, with efforts being made to establish linkages among these phases. Many methods for predicting the rheology properties of asphalt mastics from those of asphalt and filler volume fractions exist. However, most prediction methods are based on an empirical formula and on the micromechanical model. Very few research studies focus on the constitutive model. In addition, relatively little research has explored the influence of asphalt-filler interaction on mastic's rheology properties, which is believed to be an important factor. In this study, the 2S2P1D (two springs, two parabolic elements, and one dashpot) model was applied to link the behavior of asphalt binder, filler volume fraction, asphalt-filler interaction and asphalt mastic. First, the interaction between asphalt and filler was evaluated, and the interaction parameter C of the Palierne model was used as an assessment indicator to calculate the effective filler volume fraction of asphalt mastic. Then, the relation between the 2S2P1D model parameters of asphalt mastic and those of asphalt binder and the effective filler volume fraction was analyzed. Finally, a simple relationship associating the 2S2P1D model parameters h, log(τ0) of mastic and that of asphalt binder and the effective filler volume fraction was developed. The proposed expression was validated, and the result showed that it was an efficient model for the shear complex modulus prediction of virgin asphalt mastic.

Effect of IL4 and IL10 on a human in vitro type 1 diabetes model.

We defined the effect of the anti-inflammatory cytokines IL4 and IL10 on an in vitro model of human T1D. After preincubation with IL4 or IL10, human islet microtissues were co-cultured with PBMC and proinflammatory cytokines for a few hours or for multiple days to assess acute and chronic effects. This resulted in an immune attack with infiltration of T cells into the islet, a loss of beta cell endocrine function, and an upregulation of HLA-I on the beta cells. HLA-I upregulation was associated with infiltration of T cells and both HLA-I expression and infiltration were associated with impaired insulin secretion. Preincubation with IL4 or IL10 did not preserve beta cell function but decreased infiltration of T cells. Our data support the hypothesis that a loss of beta cell endocrine function mediates an increase in beta cell specific antigen presentation. IL4 and IL10 failed to preserve beta cell endocrine function.

Fast-Acquiring High-Quality Prony Series Parameters of Asphalt Concrete through Viscoelastic Continuous Spectral Models.

Prony series representations have been extensively applied to characterizing the time-domain linear viscoelastic (LVE) material functions for asphalt concrete. However, existing methods that can generate high-quality Prony series parameters (i.e., discrete spectra) mostly involve complicated programming algorithms, which poses a challenge for quick access of Prony series parameters. Also, very limited research has been devoted to establishing methods for simultaneously determining both retardation and relaxation spectra. To resolve these issues, this study presented a practical approach to fast acquiring high-quality Prony series parameters for both relaxation modulus and creep compliance of asphalt concrete by using the complex modulus test data. The approach adopts the analytical representations of the continuous relaxation and retardation spectra from the Havriliak-Negami (HN) and 2S2P1D complex modulus models to directly determine the discrete spectra, and the elastic constants, Ee and Dg, for both LVE modulus and compliance functions are further calculated by fitting the corresponding generalized Maxwell model representations to smoothed data from the storage modulus representations of the HN and 2S2P1D complex modulus models. In this way, all the procedures in the proposed method can be easily implemented in Microsoft Excel. The results showed that the HN and 2S2P1D models yielded slightly different continuous spectral patterns at shorter relaxation times and longer retardation times. However, at the region covered by the test data, the continuous spectra of the two complex modulus models were very close to each other. Thus, the two models can generate comparable Prony series parameters within the time or frequency range covered by the test data. Considering that the quality of the resulting Prony series parameters are closely related to the master curve models used for presmoothing, the HN and 2S2P1D models were compared with the conventional Sigmoidal model. Additionally, the Black diagram was recommended for examining the quality of the complex modulus test data before constructing the master curves.

Uranium isotopes in groundwater in Ho Chi Minh City and related issues: Health risks, environmental effects, and mitigation methods.

Groundwater is regularly used for many purposes, such as drinking and agricultural irrigation systems. Still, it contains high levels of radionuclides (e.g., 238U, 232Th, and 226Ra) that are potentially hazardous to humans and the environment. In this study, activity concentrations of uranium isotopes were analyzed in 15 groundwater samples taken from 15 bored wells in Thu Duc district, Ho Chi Minh City, Vietnam. Environmental effects of the irrigation system with groundwater on agricultural soil in the study area were assessed by models. It was found that the activity concentrations of 238U and 234U in groundwater samples were in the ranges of (13.5-268.7) mBq l-1 and (20.2-438.3) mBq l-1, respectively. The ratio 234U/238U values were ranged from 1.12 to 2, with an average value of 1.44. Based on the model prediction, 25 years irrigation with the groundwater can inject 94.8 Bq both uranium isotopes in 1 kg topsoil. For investigated groundwater samples, the proposed removal method using K2FeO4 removed 74.28% and 81.04% for 234U and 238U, respectively.

T1D -weighted ihMT imaging - Part I. Isolation of long- and short-T1D components by T1D -filtering.

PURPOSE: To identify T1D -filtering methods, which can specifically isolate various ranges of T1D components as they may be sensitive to different microstructural properties. METHODS: Modified Bloch-Provotorov equations describing a bi-T1D component biophysical model were used to simulate the inhomogeneous magnetization transfer (ihMT) signal from ihMTRAGE sequences at high RF power and low duty-cycle with different switching time values for the dual saturation experiment: Δt = 0.0, 0.8, 1.6, and 3.2 ms. Simulations were compared with experimental signals on the brain gray and white matter tissues of healthy mice at 7T. RESULTS: The lengthening of Δt created ihMT high-pass T1D -filters, which efficiently eliminated the signal from T1D components shorter than 1 ms, while partially attenuating that of longer components (≥ 1 ms). Subtraction of ihMTR images obtained with Δt = 0.0 ms and Δt = 0.8 ms generated a new ihMT band-pass T1D -filter isolating short-T1D components in the 100-µs to 1-ms range. Simulated ihMTR values in central nervous system tissues were confirmed experimentally. CONCLUSION: Long- and short-T1D components were successfully isolated with high RF power and low duty-cycle ihMT filters in the healthy mouse brain. Future studies should investigate the various T1D -range microstructural correlations in in vivo tissues.

Effects of the Water-Based Foaming Process on the Basic and Rheological Properties of Bitumen 70/100.

The process of water-based foaming of bitumen produces binders that can be incorporated in cold recycled asphalt mixes and pavement upper structural layers made of half-warm mix asphalt prepared at 100-130 °C. During the foaming process, cold water and air act on hot bitumen (160-170 °C), which results in the explosive vaporization of water leading to changes in the binder structure. The impact of foaming on the properties of bitumen 70/100 was evaluated by investigating the binder characteristics before and after foaming. Determination of two foaming parameters, maximum expansion and half-life, was followed by measurements of penetration at 25 °C, softening point, Fraass breaking point, and dynamic viscosity at 60, 90, and 135 °C. Rheological and low-temperature tests were also performed before and after foaming bitumen 70/100. The Bending Beam Rheometer method was applied to determine the low temperature stiffness modulus. A DHR-2 rheometer was used to determine the dynamic modulus and phase angle of the tested binder. The Black and master curves before and after foaming were plotted in the 2S2P1D model and the model parameters were analysed. Analysis of the test results confirmed the effects of the foaming process on the basic, low-temperature, and rheological characteristics of the bitumen.

Influence of cross-sectional velocity profile on flow characteristics of arterial wall modeled as elastic and viscoelastic material.

In the present work, blood flow behavior in a single artery and in arterial network is studied using time domain based one-dimensional wave propagation model retaining the nonlinear convective force. 1-D Navier-Stokes equation is used to model the flow behavior of the blood, using three unknown parameters: flow rate (q), cross-sectional area of artery (A) and pressure (p) based formulation. Three different approximate velocity profile functions across the cross-section namely modified flat, parabolic and the one proposed by Bessems are used to calculate the nonlinear convective force and the frictional force. Two different constitutive models, linear elastic model and standard linear solid model (Zener model) are used to model the arterial wall mechanical behavior. The system of partial differential equations is discretized using finite element and Crank Nicolson methods in space and time domains, respectively. Based on the formulation, an in house finite element code is developed to simulate flow behavior in both a single artery as well as in arterial network consisting of 20 small and large size arteries. Simulations are performed by enforcing a flow rate at the inlet and Windkessel model at the outlet. The results for elastic arterial wall model are found to be in good agreement with the results available in the literature. The flow rate/pressure predictions using different velocity profile functions are found to be nearly the same, however the Bessems velocity profile predicts more closer to 3D results compared to modified flat and parabolic profiles. Whereas, significant difference is found in the results predicted using elastic and viscoelastic artery wall models.