Browsing by Author "Yang, SY"

Now showing 1 - 20 of 20

A closed form solution for constant flux pumping in a well under partial penetration condition
(2006-05-19)
Yang, SY
;
Yeh, HD
;
Chiu, PY
;
環境工程研究所
;
Institute of Environmental Engineering
[ 1] An analytical model for the constant flux pumping test is developed in a radial confined aquifer system with a partially penetrating well. The Laplace domain solution is derived by the application of the Laplace transforms with respect to time and the finite Fourier cosine transforms with respect to the vertical coordinates. A time domain solution is obtained using the inverse Laplace transforms, convolution theorem, and Bromwich integral method. The effect of partial penetration is apparent if the test well is completed with a short screen. An aquifer thickness 100 times larger than the screen length of the well can be considered as infinite. This solution can be used to investigate the effects of screen length and location on the drawdown distribution in a radial confined aquifer system and to produce type curves for the estimation of aquifer parameters with field pumping drawdown data.
A least-squares finite element method for incompressible flow in stress-velocity-pressure version
(1995-12-01)
Chang, CL
;
Yang, SY
;
Hsu, CH
;
交大名義發表
;
National Chiao Tung University
In this paper we are concerned with the incompressible flow in 2-D. Introducing additional variables of derivatives of velocity, which are called stresses here, the second-order dynamic equations are reduced into a first-order system with variables of stress, velocity and pressure. Combining the compatibility condition and the divergence free condition, we have a system with six first-order equations and six unknowns. Least-squares method is performed over this extended system. The analysis shows that this method achieves optimal rates of convergence in the H-1-norm as the h approaches to zero. Numerical experiences are also available.
A least-squares finite element method for incompressible flow in stress-velocity-pressure version
(1995-12-01)
Chang, CL
;
Yang, SY
;
Hsu, CH
;
交大名義發表
;
National Chiao Tung University
In this paper we are concerned with the incompressible flow in 2-D. Introducing additional variables of derivatives of velocity, which are called stresses here, the second-order dynamic equations are reduced into a first-order system with variables of stress, velocity and pressure. Combining the compatibility condition and the divergence free condition, we have a system with six first-order equations and six unknowns. Least-squares method is performed over this extended system. The analysis shows that this method achieves optimal rates of convergence in the H-1-norm as the h approaches to zero. Numerical experiences are also available.
A new closed-form solution for a radial two-layer drawdown equation for groundwater under constant-flux pumping in a finite-radius well
(2003-07-01)
Yeh, HD
;
Yang, SY
;
Peng, HY
;
環境工程研究所
;
Institute of Environmental Engineering
A mathematical model is presented for describing the groundwater flow in a radial two-layer confined aquifer system with a constant-flux pumping well that has a wellbore skin and finite well radius. The Laplace-domain solution for the model is first derived by the Laplace transforms; and the time-domain solution in terms of the aquifer drawdown is then obtained form the Laplace inversion using the Bromwich integral method. When neglecting the well radius, our Laplace-domain solution is shown to reduce to a Laplace-domain solution given by Butler [J. Hydrol. 101 (1988) 15]. A unified numerical approach including a root search approach, the Gaussian quadrature, and the Shanks method is employed for evaluating this time-domain solution. The evaluated results of the solution agree well with those of the Laplace-domain solution estimated by the modified Crump algorithm. This new solution can be used either to predict the spatial and temporal drawdown distributions in both the skin and formation zones or to investigate the effects of the skin type, skin thickness and well radius on the drawdown distribution. (C) 2003 Elsevier Science Ltd. All rights reserved.
A simple approach using Bouwer and Rice's method for slug test data analysis
(2004-09-01)
Yang, SY
;
Yeh, HD
;
環境工程研究所
;
Institute of Environmental Engineering
Slug test data obtained from tests performed in an unconfined aquifer are commonly analyzed by graphical or numerical approaches to determine the aquifer parameters. This paper derives three fourth-degree polynomials to represent the relationship between Bouwer and Rice's coefficients and the ratio of the screen length to the radius of the gravel envelope. A numerical approach using the nonlinear least squares and Newton's method is used to determine hydraulic conductivity from the best fit of the slug test data. The method of nonlinear least squares minimizes the sum of the squares of the differences between the predicted and observed water levels inside the well. With the polynomials, the hydraulic conductivity can be obtained by simply solving the nonlinear least squares equation by Newton's method. A computer code, SLUGBR, was developed from the derived polynomials using the proposed numerical approach. The results of analyzing two slug test datasets show that SLUGBR can determine hydraulic conductivity with very good accuracy.
A two-stage least-squares finite element method for the stress-pressure-displacement elasticity equations
(1998-05-01)
Yang, SY
;
Chang, CL
;
應用數學系
;
Department of Applied Mathematics
A new stress-pressure-displacement formulation for the planar elasticity equations is proposed by introducing the auxiliary variables, stresses, and pressure. The resulting first-order system involves a nonnegative parameter that measures the material compressibility for the elastic body. A two-stage least-squares finite element procedure is introduced for approximating the solution to this system with appropriate boundary conditions. It is shown that the two-stage least-squares scheme is stable and, with respect to the order of approximation for smooth exact solutions, the rates of convergence of the approximations for all the unknowns are optimal both in the H-1-norm and in the L-2-norm. Numerical experiments with various values of the parameter are examined, which demonstrate the theoretical estimates. Among other things, computational results indicate that the behavior of convergence is uniform in the nonnegative parameter. (C) 1998 John Wiley & Sons, Inc.
A unified analysis of a weighted least squares method for first-order systems
(1998-06-01)
Yang, SY
;
Liu, JL
;
應用數學系
;
Department of Applied Mathematics
A unified analysis of a weighted least squares finite element method (WLSFEM) for approximating solutions of a large class of first-order differential systems is proposed. The method exhibits several advantageous features. For example, the trial and test functions are not required to satisfy the boundary conditions. Its discretization results in symmetric and positive definite algebraic systems with condition number O(h(-2) + w(2)). And a single piecewise polynomial finite element space may be used for all test and trial functions. Asymptotic convergence of the least squares approximations with suitable weights is established in a natural norm without requiring extra smoothness of the solutions. If instead, the solutions are sufficiently regular, a priori error estimates can be derived under two suitable assumptions which are related respectively to the symmetric positive systems of Friedrichs and first-order Agmon-Douglis-Nirenberg (ADN) elliptic systems. Numerous model problems fit into these two important systems. Some selective examples are examined and verified in the unified framework. (C) 1998 Published by Elsevier Science Inc. All rights reserved.
Analysis of least squares finite element methods for a parameter-dependent first-order system
(1998-01-01)
Yang, SY
;
Liu, JL
;
應用數學系
;
Department of Applied Mathematics
A parameter-dependent first-order system arising from elasticity problems is introduced. The system corresponds to the 2D isotropic elasticity equations under a stress-pressure-displacement formulation in which the nonnegative parameter measures the material compressibility for the elastic body. Standard and weighted least squares finite element methods are applied to this system, and analyses for different values of the parameter are performed in a unified manner. The methods offer certain advantages such as they need not satisfy the Babuska-Brezzi condition, a single continuous piecewise polynomial space can be used for the approximation of all the unknowns, the resulting algebraic system is symmetric and positive definite, accurate approximations of all the unknowns can be obtained simultaneously, and, especially, computational results do not exhibit any significant numerical locking as the parameter tends to zero which corresponds to the incompressible elasticity problem (or equivalently, the Stokes problem). With suitable boundary conditions, it is shown that both methods achieve optimal rates of convergence in the H-1-norm and in the L-2-norm for all the unknowns. Numerical experiments with various values of the parameter are given to demonstrate the theoretical estimates.
Analysis of least squares finite element methods for a parameter-dependent first-order system
(1998-01-01)
Yang, SY
;
Liu, JL
;
應用數學系
;
Department of Applied Mathematics
A parameter-dependent first-order system arising from elasticity problems is introduced. The system corresponds to the 2D isotropic elasticity equations under a stress-pressure-displacement formulation in which the nonnegative parameter measures the material compressibility for the elastic body. Standard and weighted least squares finite element methods are applied to this system, and analyses for different values of the parameter are performed in a unified manner. The methods offer certain advantages such as they need not satisfy the Babuska-Brezzi condition, a single continuous piecewise polynomial space can be used for the approximation of all the unknowns, the resulting algebraic system is symmetric and positive definite, accurate approximations of all the unknowns can be obtained simultaneously, and, especially, computational results do not exhibit any significant numerical locking as the parameter tends to zero which corresponds to the incompressible elasticity problem (or equivalently, the Stokes problem). With suitable boundary conditions, it is shown that both methods achieve optimal rates of convergence in the H-1-norm and in the L-2-norm for all the unknowns. Numerical experiments with various values of the parameter are given to demonstrate the theoretical estimates.
As-quenched microstructures of Cu3-XMnxAl alloys
(2006-08-01)
Yang, SY
;
Liu, TF
;
材料科學與工程學系
;
Department of Materials Science and Engineering
In the as-quenched condition, the microstructure of the Cu2.9Mn0.1 Al alloy was of DO3 phase containing gamma(')(1) martensite. This is similar to that reported by other workers in the Cu3Al alloy. However, the as-quenched microstructure of the Cu2.8Mn0.2Al or Cu2.7Mn0.3Al alloy were found to be of DO3 phase containing extremely fine L-J precipitates, whereas for Cu2.6Mn0.4Al alloy, it was a mixture of (DO3 + L2(1) + L-J) phases. These results are different from those proposed by Bouchard et al. (c) 2005 Elsevier B.V. All rights reserved.
Improved numerical evaluation of the radial groundwater flow equation
(2002-06-01)
Peng, HY
;
Yeh, HD
;
Yang, SY
;
環境工程研究所
;
Institute of Environmental Engineering
The closed-form solution for the hydraulic head, derived for the radial groundwater flow equation subject to the constant-head boundary condition at the wellbore, is in an integral form that covers an integration range from zero to infinity. The integral is difficult to evaluate due to the integrand not only consisting of the product and the square of the Bessel functions but also having a singularity at the origin. A unified numerical method is proposed to evaluate the solution with accuracy to five decimal places and for very wide ranges of dimensionless distances and times. This approach includes a singularity removal scheme, Newton's method, the Gaussian quadrature, and Shanks' method. It gives the dimensionless heads in tabular forms with better accuracy while comparing to those given by other approaches. A formula describing the flow rate across the wellbore, derived by the authors based on Darcy's law, is proved to equal those presented by Jaeger [Proc Royal Soc Edinburgh 61 (1942) 223] and Jacob and Lohman [Am Geo Union 33 (4) (1952) 559]. The same singularity removal scheme and the Gaussian quadrature are also employed to evaluate the wellbore flow rate. Computed values of dimensionless flow rate versus dimensionless time expressed in tabular forms are also correct to five decimal places. Both the tabular results of dimensionless hydraulic head and dimensionless flow rate may be useful in engineering applications. (C) 2002 Elsevier Science Ltd. All rights reserved.
Laplace-domain solutions for radial two-zone flow equations under the conditions of constant-head and partially penetrating well
(2005-03-01)
Yang, SY
;
Yeh, HD
;
環境工程研究所
;
Institute of Environmental Engineering
A mathematical model is presented for a constant-head test performed in a partially penetrating well with a finite-thickness skin. The model uses a no-flow boundary condition for the casing and a constant-head boundary condition for the screen to represent the partially penetrating well. The Laplace-domain solutions for the dimensionless flow rate at the wellbore and the hydraulic heads in the skin and formation zones are derived using the Laplace and finite Fourier cosine transforms. The solutions of hydraulic heads have been shown to satisfy the governing equations, related boundary conditions, and continuity requirements for the pressure head and flow rate at the interface of the skin zone and undisturbed formation. In addition, an efficient algorithm for evaluating those solutions is also presented. The dimensionless flow rates obtained from new solutions have been shown to be better than those of Novakowski's solutions, especially when the penetration ratio is large.
Least-squares finite element approximations to the Timoshenko beam problem
(2000-10-06)
Jou, J
;
Yang, SY
;
應用數學系
;
Department of Applied Mathematics
In this paper a least-squares finite element method for the Timoshenko beam problem is proposed and analyzed. The method is shown to be convergent and stable without requiring extra smoothness of the exact solutions. For sufficiently regular exact solutions, the method achieves optimal order of convergence in the H-1-norm for all the unknowns (displacement, rotation, shear, moment), uniformly in the small parameter which is generally proportional to the ratio of thickness to length. Thus the locking phenomenon disappears as the parameter tends to zero, A sharp a posteriori error estimator which is exact in the energy norm and equivalent in the H-1-norm is also briefly discussed. (C) 2000 Published by Elsevier Science Inc. All rights reserved. AMS classification: 65N15; 65N30.
Least-squares finite element methods for the elasticity problem
(1997-12-18)
Yang, SY
;
Liu, JL
;
應用數學系
;
Department of Applied Mathematics
A new first-order system formulation for the linear elasticity problem in displacement-stress form is proposed. The formulation is derived by introducing additional variables of derivatives of the displacements, whose combinations represent the usual stresses. Standard and weighted least-squares finite element methods are then applied to this extended system. These methods offer certain advantages such as that they need not satisfy the inf-sup condition which is required in the mixed finite element formulation, that a single continuous piecewise polynomial space can be used for the approximation of all the unknowns, that the resulting algebraic systems are symmetric and positive definite, and that accurate approximations of the displacements and the stresses can be obtained simultaneously. With displacement boundary conditions, it is shown that both methods achieve optimal rates of convergence in the H-1-norm and in the L-2-norm for all the unknowns. Numerical experiments with various Poisson ratios are given to demonstrate the theoretical error estimates.
Least-squares finite element methods for the elasticity problem
(1997-12-18)
Yang, SY
;
Liu, JL
;
應用數學系
;
Department of Applied Mathematics
A new first-order system formulation for the linear elasticity problem in displacement-stress form is proposed. The formulation is derived by introducing additional variables of derivatives of the displacements, whose combinations represent the usual stresses. Standard and weighted least-squares finite element methods are then applied to this extended system. These methods offer certain advantages such as that they need not satisfy the inf-sup condition which is required in the mixed finite element formulation, that a single continuous piecewise polynomial space can be used for the approximation of all the unknowns, that the resulting algebraic systems are symmetric and positive definite, and that accurate approximations of the displacements and the stresses can be obtained simultaneously. With displacement boundary conditions, it is shown that both methods achieve optimal rates of convergence in the H-1-norm and in the L-2-norm for all the unknowns. Numerical experiments with various Poisson ratios are given to demonstrate the theoretical error estimates.
On the convergence and stability of the standard least squares finite element method for first-order elliptic systems
(1998-07-01)
Yang, SY
;
應用數學系
;
Department of Applied Mathematics
A general framework of the theoretical analysis for the convergence and stability of the standard least squares finite element approximations to boundary value problems of first-order linear elliptic systems is established in a natural norm. With a suitable density assumption, the standard least squares method is proved to be convergent without requiring extra smoothness of the exact solutions. The method is also shown to be stable with respect to the natural norm. Some representative problems such as the grad-div type problems and the Stokes problem are demonstrated. (C) 1998 Published by Elsevier Science Inc. All rights reserved.
Phase transformations in a Cu-25at.%Al-7.5at.%Mn alloy
(2006-03-01)
Yang, SY
;
Liu, TF
;
材料科學與工程學系
;
Department of Materials Science and Engineering
The as-quenched microstructure of the Cu-25at.%Al-7.5at.%Mn alloy was D0(3) phase containing extremely fine L-J precipitates (the L-J phase is a new type of precipitate, which was firstly observed and identified by Liu and Jena (designated L-J phase) in a Cu-2.2-Mn0.8Al alloy.). When the as-quenched alloy was aged at temperatures ranging from 500 degrees C to 700 degrees C, the phase transition sequence was found to be (gamma-brass + L-J + D0(3)) (gamma-brass + L-J + B2) -> beta This result is different from that reported by previous workers in Cu-25at.%Al-(0-8)at.%Mn alloys. (c) 2005 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.
Phase transformations in a Cu-35 at.% Mn-25 at.% Al alloy
(2006-06-29)
Yang, SY
;
Liu, TF
;
材料科學與工程學系
;
Department of Materials Science and Engineering
The phase transformations in the Cu-35 at.% Mn-25 at.% Al alloy have been investigated by means of transmission electron microscopy (TEM) and energy-dispersive X-ray spectrometry (EDS). In as-quenched condition, the microstructure of the alloy was a mixture of (L2(1)+B2+L-J) phases. This is different from that observed by previous workers in the Cu3-xMnxAl alloys with x <= 1.0. When the as-quenched alloy was aged at 460 degrees C for short times, gamma-brass precipitates started to occur at anti-phase boundaries (APBs). After prolonged aging at 460 degrees C, the gamma-brass precipitates grew and beta-Mn precipitates were formed at the regions contiguous to the gamma-brass precipitates. The coexistence of (gamma-brass + beta-Mn) has never been observed by previous workers in Cu-Mn-Al alloy systems before. (c) 2005 Elsevier B.V. All rights reserved.
Solution for flow rates across the wellbore in a two-zone confined aquifer
(2002-02-01)
Yang, SY
;
Yeh, HD
;
環境工程研究所
;
Institute of Environmental Engineering
A closed-form solution for transient flow rates across the wellbore in a confined aquifer is derived from a two-zone radial ground-water flow equation subject to the boundary condition of keeping a constant head at the well radius. An aquifer may be considered as a two-zone system if the formation properties near the wellbore are significantly changed due to the well construction and/or well development. An efficient numerical approach is used to evaluate this newly derived solution. Values of the transient flow rate are provided in a tabular form and compared with those obtained by numerical inversion for the Laplace-domain solution. The results show that the two solutions are in good agreement. This newly derived solution can be used not only for predicting the transient flow rate across the wellbore but also for identifying the effects of a skin with a finite thickness on the estimation of transient flow rates in a ground-water system with two different formation properties.
Stability investigation of single-wafer process by using a spin etcher
(2003-07-01)
Kang, TK
;
Wang, CC
;
Tsui, BY
;
Yang, WL
;
Chien, FT
Etch rate instability was found in a spin etcher. When chemical was pumped to the spin etcher, the flow rate became unstable; possibly resulting in a chemical deficiency due to low mass transport at the wafer edge. Consequently, the etch rate and etch uniformity are unstable. The detailed etch mechanism has been investigated and the mass transport deficiency can be avoided. Less than 2% of etch nonuniformity and a stable process can be achieved by achieving a stable flow rate and control of the surface reaction over the whole wafer in spin etching. (C) 2003 The Electrochemical Society.