We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point. 相似文献
Plane ideal incompressible flow in a rectangular channel partitioned by a thin permeable barrier (lattice) is considered. In flowing through the lattice the stream suddenly (jumpwise) changes direction and loses energy. The flow is assumed to be vortical; the vorticity is discontinuous on the lattice. A mathematical formulation of the problem for the stream function is proposed in the form of a nonlinear elliptic equation with coefficients discontinuous on the lattice line. A numerical solution is constructed using the finite-element iteration method. The results of the numerical simulation show how the flow velocity profile in the channel can be controlled by means of permeable barriers. 相似文献
The local and the terminal velocities, the size and the degree of bubbles’ shape deformations were determined as a function of distance from the position of the bubble formation (capillary orifice) in solutions of n-octyltrimethylammonium bromide, n-octyldimethylphosphine oxide, n-octyl-β-D-glucopyranoside and n-octanoic acid.
These surface-active compounds have different polar groups but an identical hydrocarbon chain (C8) in the molecule. The motion of the bubbles was monitored and recorded using a stroboscopic illumination, a CCD camera, and a JVC professional video. The recorded bubble images were analyzed by the image analysis software. The bubbles accelerated rapidly and their shape was deformed immediately after detachment from the capillary. The extent of the bubbles’ shape deformation (ratio of horizontal and vertical diameters) was 1.5 in distilled water and dropped rapidly down to a level of ca. 1.05–1.03 with increasing surfactant concentration. After the acceleration period the bubbles either attained a constant value of the terminal velocity (distilled water and high concentrations of the solutions), or a maximum in the velocity profiles was observed (low concentrations). The values of the terminal velocity diminished drastically with increasing concentration, from the value of 35 cm/s in water down to about 15 cm/s, while the bubble diameter decreased by ca. 10% only. The surfactant adsorption at the surface of the bubbles was evaluated and the minimum adsorption coverages required to immobilize the bubbles’ surface were determined. It was found that this minimum adsorption coverage was ca. 4% for n-octyldimethylphosphine oxide, n-octyl-β-D-glucopyranoside, n-octanoic acid and 25% for n-octyltrimethylammonium bromide. The difference in the adsorption coverage together with the surfactants’ surface activities indicate that it is mainly the adsorption kinetics of the surfactants that governs the fluidity of interfaces of the rising bubbles. 相似文献