923 resultados para mesh-free method
Resumo:
Purpose - In this study we aim to validate a method to assess the impact of reduced visual function and observer performance concurrently with a nodule detection task. Materials and methods - Three consultant radiologists completed a nodule detection task under three conditions: without visual defocus (0.00 Dioptres; D), and with two different magnitudes of visual defocus (−1.00 D and −2.00 D). Defocus was applied with lenses and visual function was assessed prior to each image evaluation. Observers evaluated the same cases on each occasion; this comprised of 50 abnormal cases containing 1–4 simulated nodules (5, 8, 10 and 12 mm spherical diameter, 100 HU) placed within a phantom, and 25 normal cases (images containing no nodules). Data was collected under the free-response paradigm and analysed using Rjafroc. A difference in nodule detection performance would be considered significant at p < 0.05. Results - All observers had acceptable visual function prior to beginning the nodule detection task. Visual acuity was reduced to an unacceptable level for two observers when defocussed to −1.00 D and for one observer when defocussed to −2.00 D. Stereoacuity was unacceptable for one observer when defocussed to −2.00 D. Despite unsatisfactory visual function in the presence of defocus we were unable to find a statistically significant difference in nodule detection performance (F(2,4) = 3.55, p = 0.130). Conclusion - A method to assess visual function and observer performance is proposed. In this pilot evaluation we were unable to detect any difference in nodule detection performance when using lenses to reduce visual function.
Resumo:
This dissertation covers two separate topics in statistical physics. The first part of the dissertation focuses on computational methods of obtaining the free energies (or partition functions) of crystalline solids. We describe a method to compute the Helmholtz free energy of a crystalline solid by direct evaluation of the partition function. In the many-dimensional conformation space of all possible arrangements of N particles inside a periodic box, the energy landscape consists of localized islands corresponding to different solid phases. Calculating the partition function for a specific phase involves integrating over the corresponding island. Introducing a natural order parameter that quantifies the net displacement of particles from lattices sites, we write the partition function in terms of a one-dimensional integral along the order parameter, and evaluate this integral using umbrella sampling. We validate the method by computing free energies of both face-centered cubic (FCC) and hexagonal close-packed (HCP) hard sphere crystals with a precision of $10^{-5}k_BT$ per particle. In developing the numerical method, we find several scaling properties of crystalline solids in the thermodynamic limit. Using these scaling properties, we derive an explicit asymptotic formula for the free energy per particle in the thermodynamic limit. In addition, we describe several changes of coordinates that can be used to separate internal degrees of freedom from external, translational degrees of freedom. The second part of the dissertation focuses on engineering idealized physical devices that work as Maxwell's demon. We describe two autonomous mechanical devices that extract energy from a single heat bath and convert it into work, while writing information onto memory registers. Additionally, both devices can operate as Landauer's eraser, namely they can erase information from a memory register, while energy is dissipated into the heat bath. The phase diagrams and the efficiencies of the two models are solved and analyzed. These two models provide concrete physical illustrations of the thermodynamic consequences of information processing.
Resumo:
Introduction. Endoscopic hernia repair methods have become increasingly popular over the past 15 years. Nonetheless, there is no consensus regarding an optimal fixation method. Transabdominal sutures and titanium tacks or staples are the most traditional ones. Case report. We present a case of mechanic small bowel obstruction due to mesh migration occurring one year and a half after incisional hernia repair with polytetrafluoroethylene mesh fixed by spiral tacks. Discussion. Titanium spiral tacks are dangerous because of their sharp components, which can damage organs such as the small intestine, by causing microperforations. The type of prosthesis used has also contributed to the intraluminal migration, since polytetrafluoroethylene mesh is very flexible and poorly integrates in the abdominal wall. Conclusion. A prosthesis of a different material combined with a different fixation system such as absorbable tacks, biological glue, or mechanical tacks without sharp components, would have obviated mesh migration. .
Resumo:
Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.
Resumo:
Background and objectives Central venous catheterization of the internal jugular vein is a commonly performed invasive procedure associated with a significant morbidity and even mortality. Ultrasound-guided methods have shown to significantly improve the success of the technique and are recommended by various scientific societies, including the American Society of Anesthesiologists. The aim of this report is to describe an innovative ultrasound-guided central line placement of the internal jugular vein. Technique The authors describe an innovative ultrasound-guided central line placement of the internal jugular vein based on an oblique approach – the “Syringe-Free” approach. This technique allows immediate progression of the guide wire in the venous lumen, while maintaining a real-time continuous ultrasound image. Conclusions The described method adds to the traditional oblique technique the possibility of achieving a continuous real-time ultrasound-guided venipuncture and a guide wire insertion that does not need removing the probe from the puncture field, while having a single operator performing the whole procedure.
Resumo:
Khark & Kharko Islands are the last Northern point for fringing coral reefs in Iranian side of the Persian Gulf. These Coralline habitats are the Protected Area and Wildlife Refugees with the total area of 2400 ha which located in the territory of Bushehr Province. This research carried out during 2006-2007 with monthly sampling from 12 stations, which selected around Islands and inshore waters with maximum depth of 20 meter. Sampling was conducted using by Bongo-Net plankton sampler with 500μ of mesh size. Totally, 1808 specimen from 45 family fish larvae was identified in studied area, including: 21 coralline fish larva families and 24 shore fish larvae such as pelagic and demersal fishes which some of them known as indicator, sentinel or endemic species for coral reef ecosystems. The results was shown that coral reef diversity in coral reefs (Khark & Kharko Islands) is more than other habitats such as estuary and river mouth, creeks, mangrove forest sites, and off shore water of the Persian Gulf and Oman Sea Iranian side. Among Identified families, Clupeidae, Blenniidae, Sillaginidae, Atherinidae and Tripterygiidae; with more abundance were dominant families in studied area. The pick of fish larvae abundance family were estimated in spring. There were significant differences between seasonally abundance and sub areas, but there were not significant differences in diversity indexes between Khark and Kharko stations with coastal stations (p< 0.05). The mean abundance of fish larvae were estimated 18.7083 larvae under 10m² of sea surface, and the mean diversity indexes and evenness were estimated 0.7135 and 0.565342 consequently, that was showed the area is under ecological stress for fish larvae, and wasn’t stable. Therefore, from the ecological point of view, only some of the fish larvae groups as like Clupeidae were dominant. Thus, they were the main cause of the fish larvae abundance change in studied area. Due to geographical location of Khark and Kharko Islands and among the environmental parameters, Its seems that the condition of sea current is the main cause for present or absent and distribution patterns of fish larvae in area. Abundance of fish larvae in west of Islands was higher than eastern parts in the spring. But this condition will be reversed in eastern part of Island and several coastal stations, so that the Islands surrounding clock wise current to cause fish larvae distribution patterns.
Resumo:
Introduction: Free tissue transfer using an abdominal tissue flap is a commonly used method of breast reconstruction. However, there are well recognised complications including venous congestion, fat necrosis and flap loss associated with the perfusion of these flaps. Post-operative aesthetic outcome assessment of such breast reconstructions have also proven to be difficult with current methods displaying poor inter-rater reliability and patient correlation. The aim of this research was to investigate potential improvements to the post-operative outcome of free abdominal tissue transfer breast reconstruction by assessing the effects of vascular augmentation interventions on flap perfusion and to assess the use of real-time digital video as a post-operative assessment tool. Methods: An in-vivo pilot study carried out on 12 patients undergoing DIEP flap breast reconstruction assessed the effect on Zone IV perfusion, using LDI and ICG angiography, of vascular augmentation of the flap using the contralateral SIEA and SIEV. A further animal experimental study was carried out on 12 Sprague Dawley rats to assess the effects on main pedicle arterial blood flow and on Zone I and Zone IV perfusion of vascular augmentation of the abdominal flap using the contralateral vascular system. A separate post-operative assessment study was undertaken on 35 breast reconstruction patients who evaluated their own reconstructions via patient questionnaire and underwent photograph and real-time digital video capture of their reconstructions with subsequent panel assessment. Results: Our results showed that combined vascular augmentation of DIEP flaps, using both the SIEA and SIEV together, led to an increase in Zone IV perfusion. Vascular augmentation of the rat abdominal flaps also led to a significant increase in Zone I/IV perfusion, but the augmentation procedure resulted in a decreased main pedicle arterial blood flow. Our post-operative assessment study revealed that real-time digital video footage led to greater inter-rater agreement with regards to cosmesis and shape than photography and also correlated more with patient self-assessment. Conclusion: Vascular augmentation of abdominal free tissue flaps using the contralateral vascular system results in an increase to Zone IV perfusion, however this may lead to decreased main pedicle arterial blood flow. Real-time digital video is a valid post-operative aesthetic assessment method of breast reconstruction outcome and is superior to static photography when coupled with panel assessment.
Resumo:
This article presents the Konstanz Method of Dilemma Discussion ® (KMDD ®) and explains the integration of the KMDD ® in ethics lessons. In this paper, some special learning effects of this inclusive teaching and learning method are shown. Furthermore, it investigates the questions of how to achieve more knowledge in ethics lessons by dialogue and how to realize better moral development, particularly by handling of differentiation. Moral education of all participants who are involved in the learning process (learners and teacher alike) is a crucial task of every true inclusion. True inclusion means building optimal learning conditions in keeping with the free will of all participants. Because our society is transforming constantly in both global and demographic aspects, coping with these challenges is mandatory.
Resumo:
A new type of space debris was recently discovered by Schildknecht in near -geosynchronous orbit (GEO). These objects were later identified as exhibiting properties associated with High Area-to-Mass ratio (HAMR) objects. According to their brightness magnitudes (light curve), high rotation rates and composition properties (albedo, amount of specular and diffuse reflection, colour, etc), it is thought that these objects are multilayer insulation (MLI). Observations have shown that this debris type is very sensitive to environmental disturbances, particularly solar radiation pressure, due to the fact that their shapes are easily deformed leading to changes in the Area-to-Mass ratio (AMR) over time. This thesis proposes a simple effective flexible model of the thin, deformable membrane with two different methods. Firstly, this debris is modelled with Finite Element Analysis (FEA) by using Bernoulli-Euler theory called “Bernoulli model”. The Bernoulli model is constructed with beam elements consisting 2 nodes and each node has six degrees of freedom (DoF). The mass of membrane is distributed in beam elements. Secondly, the debris based on multibody dynamics theory call “Multibody model” is modelled as a series of lump masses, connected through flexible joints, representing the flexibility of the membrane itself. The mass of the membrane, albeit low, is taken into account with lump masses in the joints. The dynamic equations for the masses, including the constraints defined by the connecting rigid rod, are derived using fundamental Newtonian mechanics. The physical properties of both flexible models required by the models (membrane density, reflectivity, composition, etc.), are assumed to be those of multilayer insulation. Both flexible membrane models are then propagated together with classical orbital and attitude equations of motion near GEO region to predict the orbital evolution under the perturbations of solar radiation pressure, Earth’s gravity field, luni-solar gravitational fields and self-shadowing effect. These results are then compared to two rigid body models (cannonball and flat rigid plate). In this investigation, when comparing with a rigid model, the evolutions of orbital elements of the flexible models indicate the difference of inclination and secular eccentricity evolutions, rapid irregular attitude motion and unstable cross-section area due to a deformation over time. Then, the Monte Carlo simulations by varying initial attitude dynamics and deformed angle are investigated and compared with rigid models over 100 days. As the results of the simulations, the different initial conditions provide unique orbital motions, which is significantly different in term of orbital motions of both rigid models. Furthermore, this thesis presents a methodology to determine the material dynamic properties of thin membranes and validates the deformation of the multibody model with real MLI materials. Experiments are performed in a high vacuum chamber (10-4 mbar) replicating space environment. A thin membrane is hinged at one end but free at the other. The free motion experiment, the first experiment, is a free vibration test to determine the damping coefficient and natural frequency of the thin membrane. In this test, the membrane is allowed to fall freely in the chamber with the motion tracked and captured through high velocity video frames. A Kalman filter technique is implemented in the tracking algorithm to reduce noise and increase the tracking accuracy of the oscillating motion. The forced motion experiment, the last test, is performed to determine the deformation characteristics of the object. A high power spotlight (500-2000W) is used to illuminate the MLI and the displacements are measured by means of a high resolution laser sensor. Finite Element Analysis (FEA) and multibody dynamics of the experimental setups are used for the validation of the flexible model by comparing with the experimental results of displacements and natural frequencies.
Resumo:
Measuring the extent to which a piece of structural timber has distorted at a macroscopic scale is fundamental to assessing its viability as a structural component. From the sawmill to the construction site, as structural timber dries, distortion can render it unsuitable for its intended purposes. This rejection of unusable timber is a considerable source of waste to the timber industry and the wider construction sector. As such, ensuring accurate measurement of distortion is a key step in addressing ineffciencies within timber processing. Currently, the FRITS frame method is the established approach used to gain an understanding of timber surface profile. The method, while reliable, is dependent upon relatively few measurements taken across a limited area of the overall surface, with a great deal of interpolation required. Further, the process is unavoidably slow and cumbersome, the immobile scanning equipment limiting where and when measurements can be taken and constricting the process as a whole. This thesis seeks to introduce LiDAR scanning as a new, alternative approach to distortion feature measurement. In its infancy as a measurement technique within timber research, the practicalities of using LiDAR scanning as a measurement method are herein demonstrated, exploiting many of the advantages the technology has over current approaches. LiDAR scanning creates a much more comprehensive image of a timber surface, generating input data multiple magnitudes larger than that of the FRITS frame. Set-up and scanning time for LiDAR is also much quicker and more flexible than existing methods. With LiDAR scanning the measurement process is freed from many of the constraints of the FRITS frame and can be done in almost any environment. For this thesis, surface scans were carried out on seven Sitka spruce samples of dimensions 48.5x102x3000mm using both the FRITS frame and LiDAR scanner. The samples used presented marked levels of distortion and were relatively free from knots. A computational measurement model was created to extract feature measurements from the raw LiDAR data, enabling an assessment of each piece of timber to be carried out in accordance with existing standards. Assessment of distortion features focused primarily on the measurement of twist due to its strong prevalence in spruce and the considerable concern it generates within the construction industry. Additional measurements of surface inclination and bow were also made with each method to further establish LiDAR's credentials as a viable alternative. Overall, feature measurements as generated by the new LiDAR method compared well with those of the established FRITS method. From these investigations recommendations were made to address inadequacies within existing measurement standards, namely their reliance on generalised and interpretative descriptions of distortion. The potential for further uses of LiDAR scanning within timber researches was also discussed.
Resumo:
The transfer coefficients for momentum and heat have been determined for 10 m neutral wind speeds (U-10n) between 0 and 12 m/s using data from the Surface of the Ocean, Fluxes and Interactions with the Atmosphere (SOFIA) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments. The inertial dissipation method was applied to wind and pseudo virtual temperature spectra from a sonic anemometer, mounted on a platform (ship) which was moving through the turbulence held. Under unstable conditions the assumptions concerning the turbulent kinetic energy (TKE) budget appeared incorrect. Using a bulk estimate for the stability parameter, Z/L (where Z is the height and L is the Obukhov length), this resulted in anomalously low drag coefficients compared to neutral conditions. Determining Z/L iteratively, a low rate of convergence was achieved. It was concluded that the divergence of the turbulent transport of TKE was not negligible under unstable conditions. By minimizing the dependence of the calculated neutral drag coefficient on stability, this term was estimated at about -0.65Z/L. The resulting turbulent fluxes were then in close agreement with other studies at moderate wind speed. The drag and exchange coefficients for low wind speeds were found to be C-en x 10(3) = 2.79U(10n)(-1) + 0.66 (U-10n < 5.2 m/s), C-en x 10(3) = C-hn x 10(3) = 1.2 (U-10n greater than or equal to 5.2 m/s), and C-dn x 10(3) = 11.7U(10n)(-2) + 0.668 (U-10n < 5.5 m/s), which imply a rapid increase of the coefficient values as the wind decreased within the smooth flow regime. The frozen turbulence hypothesis and the assumptions of isotropy and an inertial subrange were found to remain valid at these low wind speeds for these shipboard measurements. Incorporation of a free convection parameterization had little effect.
Resumo:
Currently, many consumers search for food with functional characteristics beyond their nutritional properties. Thus, the concept of functional food becomes a hot topic, allowing the obtaining of health benefits, including disease prevention. In this context, plants are recognized as sources of a wide range of bioactives, mainly phenolic compounds. In particular, the Rosmarinus officina/is L., commonly referred as rosemary, has several phenolic compounds with different bioactive properties such as antioxidant, antiinflammatory and antimicrobial activities, among others [!]. Hence, this plant has great potential for incorporation into foods in order to confer bioactivity to the final products. However, it should be highlighted that the bioactive compounds if exposed to adverse environments, for example: light, moisture, extreme pH, storage, food processing conditions, can be degraded leading to the consequent loss of bioactivity [2]. The microencapsulation is an alternative to overcome this problematic of bioactive compounds, as also to ensure controlled release, or target deliver to a specific site [3]. In this work, lyophilized rosemary aqueous extract prepared by in:'usion was used as a functional ingredient for cottage cheeses, after proving that it possesses, both higher content in phenolic compounds and higher antioxidant activity, comparatively with the corresponding hydroethanolic extract. The rosemary aqueous extract revealed, for example, a DPPH scavenging activity with an EC50 value of 73.44±0.54j!g/mL and presented as main phenolic compound the caffeic acid dimer, commonly named as rosmarinic acid. For the functionalized cottage cheeses, a decrease of bioactivity was observed after seven days under storage in fridge, when the extracts were incorporated in its free form. Therefore, to preserve the antioxidant activity, the rosemary aqueous extract was efficiently microencapsulated by using an atomization/coagulation technique and alginate as the matrix material and thereafter incorporated into the cottage cheeses. The final microspheres showed a size, estimated by OM using a magnification of I OOx, ranging between 51.1 and 122.6 J!m and an encapsulation efficiency, estimated through an indirect method, approaching 100%. Overall, the introduction of both free and microencapsulated extracts did not change the nutritional value of cottage cheeses, providing bioactivity that was more preserved with microencapsulated extracts putting in evidence the importance of using microencapsulation to develop effective functional foods.
Resumo:
The research described in this thesis was motivated by the need of a robust model capable of representing 3D data obtained with 3D sensors, which are inherently noisy. In addition, time constraints have to be considered as these sensors are capable of providing a 3D data stream in real time. This thesis proposed the use of Self-Organizing Maps (SOMs) as a 3D representation model. In particular, we proposed the use of the Growing Neural Gas (GNG) network, which has been successfully used for clustering, pattern recognition and topology representation of multi-dimensional data. Until now, Self-Organizing Maps have been primarily computed offline and their application in 3D data has mainly focused on free noise models, without considering time constraints. It is proposed a hardware implementation leveraging the computing power of modern GPUs, which takes advantage of a new paradigm coined as General-Purpose Computing on Graphics Processing Units (GPGPU). The proposed methods were applied to different problem and applications in the area of computer vision such as the recognition and localization of objects, visual surveillance or 3D reconstruction.
Resumo:
We introduce and analyze a discontinuous Galerkin method for the numerical discretization of a stationary incompressible magnetohydrodynamics model problem. The fluid unknowns are discretized with inf-sup stable discontinuous P^3_{k}-P_{k-1} elements whereas the magnetic part of the equations is approximated by discontinuous P^3_{k}-P_{k+1} elements. We carry out a complete a-priori error analysis and prove that the energy norm error is convergent of order O(h^k) in the mesh size h. We also show that the method is able to correctly capture and resolve the strongest magnetic singularities in non-convex polyhedral domains. These results are verified in a series of numerical experiments.