• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.259-267
• /
• 2007

An efficient algorithm for inverting static-shifted magnetotelluric (MT) data has been proposed to produce a three-dimensional (3D) resistivity model. In the Gauss-Newton approach, computational costs associated with construction of a full sensitivity matrix usually make 3D MT inversion impractical. This computational difficulty may be overcome by using approximate sensitivities. We use four kinds of sensitivities in particular orders in the inversion process. These sensitivities are computed 1) analytically for an initial, homogeneous earth, 2) exactly for a current model, 3) approximately by the Broyden method, and 4) approximately using the previous adjoint fields. Inversion experiments with static-shifted synthetic and field MT data indicate that inversion results are highly dependent on characteristics of data and thus applying various combinations of sensitivities is helpful in obtaining a good image of the subsurface structure with reasonable computation time.

• Journal of Educational Research in Mathematics
• /
• v.19 no.4
• /
• pp.479-491
• /
• 2009

This paper researches the possibility of introducing Descartes' theorem to mathematically gifted students. Not only is Descartes' theorem logically equivalent to Euler's theorem but is hierarchically connected with Gauss-Bonnet theorem which is the core concept on differential geometry. It is possible to teach mathematically gifted students Descartes' theorem by generalizing mathematical property in solid geometry through analogical reasoning, that is, so in a polyhedrons the sum of the deficient angles is $720^\circ$ as in an polygon the sum of the exterior angles is $360^\circ$. This study introduces an alternative method of instruction that we enable mathematically gifted students to reinvent Descartes' theorem through analogical reasoning instead of deductive reasoning.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.291-309
• /
• 2023

This paper addresses the finite element modeling of functionally graded porous (FGP) beams for free vibration and buckling behaviour cases. The formulated finite element is based on simple and efficient higher order shear deformation theory. The key feature of this formulation is that it deals with Euler-Bernoulli beam theory with only three unknowns without requiring any shear correction factor. In fact, the presented two-noded beam element has three degrees of freedom per node, and the discrete model guarantees the interelement continuity by using both C⁰ and C¹ continuities for the displacement field and its first derivative shape functions, respectively. The weak form of the governing equations is obtained from the Hamilton principle of FGP beams to generate the elementary stiffness, geometric, and mass matrices. By deploying the isoparametric coordinate system, the derived elementary matrices are computed using the Gauss quadrature rule. To overcome the shear-locking phenomenon, the reduced integration technique is used for the shear strain energy. Furthermore, the effect of porosity distribution patterns on the free vibration and buckling behaviours of porous functionally graded beams in various parameters is investigated. The obtained results extend and improve those predicted previously by alternative existing theories, in which significant parameters such as material distribution, geometrical configuration, boundary conditions, and porosity distributions are considered and discussed in detailed numerical comparisons. Determining the impacts of these parameters on natural frequencies and critical buckling loads play an essential role in the manufacturing process of such materials and their related mechanical modeling in aerospace, nuclear, civil, and other structures.

• Progress in Superconductivity
• /
• v.3 no.1
• /
• pp.130-133
• /
• 2001

A multifilamentary Bi-2223 HTS tape for superconducting power applications was studied through the fabrication of 250-meter long tapes by the PIT(powder in tube) process. To fabricate continuous long wire, a drawing machine, a two-drum bull block and a rolled tape winding machine were developed. Especially, 250-meter long tapes were heat treated in the shape of pancake coil to reduce the heat affect zone and to achieve the high critical current. Engineering critical current density was improved through both the enhancements of critical current density by control of thermal process and the increase of filling factor by using thin Ag alloy sheath tubes less than 1.5 mm in thickness. We have made successfully 250-meter long 37 filamentary tapes with high filling factor up to 31 % employing the modified drawing and rolling technique. The critical current of 250-meter long tapes with pancake coil type was measured by transport method at self-field up to 250 gauss of center field. The measured values, based on the transport critical current at self-field, $I_{c}$ -B characteristics and magnetic field analysis, are 34 A of I$_{c}$ and 4.0 $kA/\textrm{cm}^2$ of $J_{e}$ at 250 m, 77 K, and 0 T. We also have achieved the 56 A of I$_{c}$ and 7.0 $0 kA/\textrm{cm}^2$ of$ J_{e}$ in short tapes at 77K, self-field, and 1$mutextrm{V}$/cm.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.413-420
• /
• 1994

Dynamics of the river plume is a very complicated non-linear problem with the free boundary changing in time and space. Mixing with the ambient water through the boundary makes the problem more complicated. In this paper we reduced 3-dimensional problem into 1-dimensional one by using the integral analysis method. Basic equations have been integrated over the lateral and vertical variations. For these integrations we adopted the well-established assumption that the flow-axis component of plume velocity and the density difference of the plume with the ambient water have Gaussian distributions in directions which are perpendicular to the flow-axis of the plume. We also used the result of our previous study on the lateral spreading velocity of the plume derived under the same assumption. And entrainment was included as a mixing process. The resultant 1-dimensional equations were solved by Runge-Kutta numerical method. Consequently, comparatively easy method of numerical analysis is presented for the 3-dimensional river plume. The method can also be used for the analysis of the thermal plume of cooling water of power plants.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.2
• /
• pp.167-180
• /
• 2008

The resolution of Earth images taken from a satellite has close relation with satellite's altitude. If a satellite has lower altitude, it gets a picture having better resolution. However the satellite will be exposed to heavier air drag and will spend more fuel to maintain its altitude for a desired mission. Therefore, in this study, the required fuel to maintain very low earth orbit(LEO) with severe air drag is analyzed using optimization method such as collocation method. The required fuel to maintain the low altitude has significantly increased as the mission altitude is lowered and the solar activity is maximized. This study also shows that the fuel reduced by increasing the period of the satellite maneuver is very small, and that slightly increasing the satellite's mission altitude is much effective in reducing the amount of fuel to maintain its altitude. The calculated fuel to maintain very low earth orbit in this study would give useful information in planning the budget of fuel and cost for LEO satellites.

• Applied Chemistry for Engineering
• /
• v.29 no.2
• /
• pp.215-224
• /
• 2018

Influences of the magnetic field on 5, 10 and 15 wt% (1.5, 3 and 4.5 M) HCl solution behaviour, which has widespread applications in petroleum well acidizing, were investigated in various conditions. Differences in the pH of magnetized hydrochloric acid compared to that of normal hydrochloric acid were measured. Taguchi design of experimental (DoE) method were used to model effects of the magnetic field intensity, concentration, velocity and temperature of acid in addition to the elapsed time. The experimental results showed that the magnetic field decreases [$H^+$] concentration of hydrochloric acid up to 42% after magnetization. Increasing the magnetic field intensity (with 28% contribution), concentration (with 42% contribution), and velocity of acid increases the effect of magnetic treatment. The results also demonstrated that the acid magnetization was-not influenced by the fluid velocity and heating. It was also displayed that the acid preserves its magnetic memory during time. The optimum combination of factors with respect to the highest change of [$H^+$] concentration was obtained as an acid concentration of 10% and an applied magnetic field of 4,300 Gauss. Due to the reduction of HCl reaction rate under the magnetization process, it can be proposed that the magnetized HCl is a cost effective and reliable alternative retarder in the matrix acidizing of hydrocarbon (crude oil and natural gas) wells.

• Journal of Radiation Protection and Research
• /
• v.17 no.1
• /
• pp.43-56
• /
• 1992

As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.49-56
• /
• 2024

An automated dynamic structural analysis module stands as a crucial element within a structural integrated mitigation system. This module must deliver prompt real-time responses to enable timely actions, such as evacuation or warnings, in response to the severity posed by the structural system. The finite element method, a widely adopted approximate structural analysis approach globally, owes its popularity in part to its user-friendly nature. However, the computational efficiency and accuracy of results depend on the user-provided finite element mesh, with the number of elements and their quality playing pivotal roles. This paper introduces a computationally efficient adaptive mesh generation scheme that optimally combines the h-method of node movement and the r-method of element division for mesh refinement. Adaptive mesh generation schemes automatically create finite element meshes, and in this case, representative strain values for a given mesh are employed for error estimates. When applied to dynamic problems analyzed in the time domain, meshes need to be modified at each time step, considering a few hundred or thousand steps. The algorithm's specifics are demonstrated through a standard cantilever beam example subjected to a concentrated load at the free end. Additionally, a portal frame example showcases the generation of various robust meshes. These examples illustrate the adaptive algorithm's capability to produce robust meshes, ensuring reasonable accuracy and efficient computing time. Moreover, the study highlights the potential for the scheme's effective application in complex structural dynamic problems, such as those subjected to seismic or erratic wind loads. It also emphasizes its suitability for general nonlinear analysis problems, establishing the versatility and reliability of the proposed adaptive mesh generation scheme.