Two different methods for numerical solution of the modified Burgers' equation.
Ontology highlight
ABSTRACT:
SUBMITTER: Karakoc SB
PROVIDER: S-EPMC4076646 | biostudies-other | 2014
REPOSITORIES: biostudies-other
Similar Datasets
Project description:To better describe its constitutive relation, we need a new constitutive equation for an important nonlinear elastic material, Mn-Cu damping alloy. In this work, we studied the nonlinear and hysteretic characteristics of the stress-strain curve of the M2052 alloy with the uniaxial cyclic tensile test with constant strain rate. The strain rate and amplitude correlations of M2052 resembled those of nonlinear viscoelastic material. Therefore, we created a new constitutive equation for the M2052 damping alloy by modifying the fractional Maxwell model, and we used the genetic algorithm to carry out numerical fitting with MATLAB. By comparing with the experimental data, we confirmed that the new constitutive equation could accurately depict the nonlinear constitutive relation and hysteretic property of the damping alloy. Taken together, this new constitutive equation for Mn-Cu damping alloy based on the fractional Maxwell model can serve as an effective tool for further studies of the constitutive relation of the Mn-Cu damping alloys.
| S-EPMC7254341 | biostudies-literature
Project description:BackgroundThe area of the left ventricular outflow tract (ALVOT) represents a major component of the continuity equation (CE), which is, i.a., crucial to calculate the aortic valve (AV) area (AAV). The ALVOT is typically calculated using 2D echo assessments as the measured anterior-posterior (a/p) extension, assuming a round LVOT base. Anatomically, however, usually an elliptical shape of the LVOT base is present, with the long diameter extending from the medial-lateral axis (m/l), which is not recognized by two-dimensional (2D) echocardiography.ObjectiveWe aimed to compare standard and three-dimensional (3D)-echocardiography-derived ALVOT calculation and its use in a standard CE (CEstd) and a modified CE (CEmod) to calculate the AAV vs. computed tomography (CT) multi-planar reconstruction (MPR) measurements of the anatomical ALVOT, and AAV, respectively.MethodsPatients were selected if 3D transthoracic echocardiography (TTE), 3D transesophageal echocardiography (TEE), and cardiac CT were all performed, and imaging quality was adequate. The ALVOT was assessed using 2D calculation, (a/p only), 3D-volume MPR, and 3D-biplane calculation (a/p and m/l). AAV was measured using both CEstd and CEmod, and 3D-volume MPR. Data were compared to corresponding CT analyses.ResultsFrom 2017 to 2018, 107 consecutive patients with complete and adequate imaging data were included. The calculated ALVOT was smaller when assessed by 2D- compared to both 3D-volume MPR and 3D-biplane calculation. Calculated AAV was correspondingly smaller in CEstd compared to CEmod or 3D-volume MPR. The ALVOT and AAV, using data from 3D echocardiography, highly correlated and were congruent with corresponding measurements in CT.ConclusionDue to the elliptic shape of the LVOT, use of measurements and calculations based on 2D echocardiography systematically underestimates the ALVOT and dependent areas, such as the AAV. Anatomically correct assessment can be achieved using 3D echocardiography and adapted calculations, such as CEmod.
| S-EPMC9321790 | biostudies-literature
Project description:We present a simple, one-dimensional model of an atom exposed to a time-dependent intense, short-pulse EM field with the objective of teaching undergraduates how to apply various numerical methods to study the behavior of this system as it evolves in time using several time propagation schemes.In this model, the exact Coulomb potential is replaced by a soft-core interaction to avoid the singularity at the origin. While the model has some drawbacks, it has been shown to be a reasonable representation of what occurs in the fully three-dimensional hydrogen atom.The model can be used as a tool to train undergraduate physics majors in the art of computation and software development.Program summaryProgram Title:: 1d hydrogen light interactionProgram Files doi:: http://dx.doi.org/10.17632/2275fmvdzc.1Code Ocean Capsule:: https://doi.org/10.24433/CO.1476487.v1Licensing provisions:: MIT licenseProgramming language:: FORTRAN90Nature of problem:: The one dimensional time dependent Schrödinger equation has been shown to be quite useful as a model to study the Hydrogen atom exposed to an intense, short pulse, electromagnetic field. We use a model potential that is cut-off near x = 0 and avoids the singularity of the true 1-D potential, but retains the characteristic Rydberg series and continuum to study excitation and ionization of the true H atom. The code employs a number of numerical methods to understand and compare the efficacy and accuracy when applied to this model problem.Solution method:: The program uses and contrasts a number of approaches; the Crank-Nicolson, Short Iterative Lanczos, various incarnations of the split-operator and the Chebychev method have been programmed. These methods have been compared using a 3-point finite difference (FD) discretization of the space coordinate. For completeness, some attention has also been given to using 5-9 FD formulas in order to show how higher order discretization affects the accuracy and efficiency of the methods but the primary focus of the method is the time propagation.Additional comments including restrictions and unusual features:: The main purpose of this code is as a teaching tool for undergraduates interested in acquiring knowledge of numerical methods and programming skills useful to a practicing computational physicist.
| S-EPMC7594019 | biostudies-literature
Project description:BackgroundThe calculation of diffusion-controlled ligand binding rates is important for understanding enzyme mechanisms as well as designing enzyme inhibitors.MethodsWe demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) BC, is considered on the reactive boundaries. This new BC treatment allows for the analysis of enzymes with "imperfect" reaction rates.ResultsThe numerical method is first verified in simple systems and then applied to the calculation of ligand binding to a mouse acetylcholinesterase (mAChE) monomer. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates.ConclusionsAlthough this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.
| S-EPMC4438506 | biostudies-literature
Project description:The purpose of this study is to compare the effectiveness (quality of bowel cleansing) and ease of use (patient satisfaction) of two methods of administration of a precolonoscopy PEG solution in hospitalized patients.
| 2051444 | ecrin-mdr-crc
Project description:Differential equations are extensively used for modeling dynamics of physical processes in many scientific fields such as engineering, physics, and biomedical sciences. Parameter estimation of differential equation models is a challenging problem because of high computational cost and high-dimensional parameter space. In this article, we propose a novel class of methods for estimating parameters in ordinary differential equation (ODE) models, which is motivated by HIV dynamics modeling. The new methods exploit the form of numerical discretization algorithms for an ODE solver to formulate estimating equations. First, a penalized-spline approach is employed to estimate the state variables and the estimated state variables are then plugged in a discretization formula of an ODE solver to obtain the ODE parameter estimates via a regression approach. We consider three different order of discretization methods, Euler's method, trapezoidal rule, and Runge-Kutta method. A higher-order numerical algorithm reduces numerical error in the approximation of the derivative, which produces a more accurate estimate, but its computational cost is higher. To balance the computational cost and estimation accuracy, we demonstrate, via simulation studies, that the trapezoidal discretization-based estimate is the best and is recommended for practical use. The asymptotic properties for the proposed numerical discretization-based estimators are established. Comparisons between the proposed methods and existing methods show a clear benefit of the proposed methods in regards to the trade-off between computational cost and estimation accuracy. We apply the proposed methods t an HIV study to further illustrate the usefulness of the proposed approaches.
| S-EPMC3496749 | biostudies-literature
Project description:The Two-Way Ranging (TWR) method is commonly used for measuring the distance between two wireless transceiver nodes, especially when clock synchronization between the two nodes is not available. For modeling the time-of-flight (TOF) error between two wireless transceiver nodes in TWR, the existing error model, described in the IEEE 802.15.4-2011 standard, is solely based on clock drift. However, it is inadequate for in-depth comparative analysis between different TWR methods. In this paper, we propose a novel TOF Error Estimation Model (TEEM) for TWR methods. Using the proposed model, we evaluate the comparative analysis between different TWR methods. The analytical results were validated with both numerical simulation and experimental results. Moreover, we demonstrate the pitfalls of the symmetric double-sided TWR (SDS-TWR) method, which is the most highlighted TWR method in the literature because of its highly accurate performance on clock-drift error reduction when reply times are symmetric. We argue that alternative double-sided TWR (AltDS-TWR) outperforms SDS-TWR. The argument was verified with both numerical simulation and experimental evaluation results.
| S-EPMC6387385 | biostudies-literature
Project description:The methods of Atkins & Nimmo (1973) and Fernley (1974) for fitting the integrated Michaelis-Menten equation were compared by using the same sets of simulated experimental data. The method of Fernley (1974) is to be preferred because it gives precise and unbiased estimates of the Michaelis-Menten parameters over a wide range of substrate concentrations. However, the estimates may not be symmetrically distributed, especially at low substrate concentrations.
| S-EPMC1168197 | biostudies-other
Project description:This study employs the principles of computer science and statistics to evaluate the efficacy of the linear random effect model, utilizing Lasso variable selection techniques (including Lasso, Elastic-Net, Adaptive-Lasso, and SCAD) through numerical simulation and empirical research. The analysis focuses on the model's consistency in variable selection, prediction accuracy, stability, and efficiency. This study employs a novel approach to assess the consistency of variable selection across models. Specifically, the angle between the actual coefficient vector β and the estimated coefficient vector βˆ is computed to determine the degree of consistency. Additionally, the boxplot tool of statistical analysis is utilized to visually represent the distribution of model prediction accuracy data and variable selection consistency. The comparative stability of each model is assessed based on the frequency of outliers. This study conducts comparative experiments of numerical simulation to evaluate a proposed model evaluation method against commonly used analysis methods. The results demonstrate the effectiveness and correctness of the proposed method, highlighting its ability to conveniently analyze the stability and efficiency of each fitting model.
| S-EPMC10495967 | biostudies-literature
Project description:The present work is concerned with the study of a generalized Langevin equation and its link to the physical theories of statistical mechanics and scale relativity. It is demonstrated that the form of the coefficients of the Langevin equation depends critically on the assumption of continuity of the reconstructed trajectory. This in turn demands for the fluctuations of the diffusion term to be discontinuous in time. This paper further investigates the connection between the scale-relativistic and stochastic mechanics approaches, respectively, with the study of the Burgers equation, which in this case appears as a stochastic geodesic equation for the drift. By further demanding time reversibility of the drift, the Langevin equation can also describe equivalent quantum-mechanical systems in a path-wise manner. The resulting statistical description obeys the Fokker-Planck equation of the probability density of the differential system, which can be readily estimated from simulations of the random paths. Based on the Fokker-Planck formalism, a new derivation of the transient probability densities is presented. Finally, stochastic simulations are compared to the theoretical results.
| S-EPMC7513018 | biostudies-literature