THEORETICAL MAGNETOTELLURIC AND TURAM RESPONSE FROM TWO‐DIMENSIONAL INHOMOGENEITIES

Geophysics ◽  
1971 ◽  
Vol 36 (1) ◽  
pp. 38-52 ◽  
Author(s):  
Charles M. Swift
Keyword(s):  
Numerical Solutions ◽  
Deep Structure ◽  
Solution Technique ◽  
Two Dimensional ◽  
Near Surface ◽  
Anomalous Field ◽  
Strike Direction ◽  
Conductivity Structure ◽  
Representative Model ◽  
Infinite Line

The “network solution” technique for obtaining numerical solutions to Maxwell’s equations in two‐dimensional inhomogeneous media, originally developed by Madden, is presented in detail in this paper. By using plane‐wave and infinite line‐current sources respectively, theoretical magnetotelluric and Turam responses can be calculated above and within an arbitrarily complex two‐dimensional earth conductivity structure. Analysis of a few representative model results emphasizes the importance of considering all the electromagnetic field components, particularly those of the anomalous field. In magnetotellurics, the vertical magnetic field is useful in determining strike direction, and only the E parallel apparent resistivities are representative of deep structure in the presence of near‐surface conductivity inhomogeneities. In Turam, a finite resistivity background causes buried conductors to appear deeper and more conductive than they would appear if the background had infinite resistivity; and strong conductors can affect the total field sufficiently to cause screening.

Two-dimensional elastic pseudo-spectral modeling of wide-aperture seismic array data with application to the Wichita Uplift-Anadarko Basin region of southwestern Oklahoma

1990 ◽  
Vol 80 (6A) ◽  
pp. 1677-1695 ◽  
Author(s):  
Ik Bum Kang ◽  
George A. McMechan
Keyword(s):  
Large Scale ◽  
Deep Structure ◽  
Sedimentary Basins ◽  
P Wave ◽  
Two Dimensional ◽  
Anadarko Basin ◽  
Near Surface ◽  
University Of Texas ◽  
Wide Aperture ◽  
The University

Abstract Full wave field modeling of wide-aperture data is performed with a pseudospectral implementation of the elastic wave equation. This approach naturally produces three-component stress and two-component particle displacement, velocity, and acceleration seismograms for compressional, shear, and Rayleigh waves. It also has distinct advantages in terms of computational requirements over finite-differencing when data from large-scale structures are to be modeled at high frequencies. The algorithm is applied to iterative two-dimensional modeling of seismograms from a survey performed in 1985 by The University of Texas at El Paso and The University of Texas at Dallas across the Anadarko basin and the Wichita Mountains in southwestern Oklahoma. The results provide an independent look at details of near-surface structure and reflector configurations. Near-surface (<3 km deep) structure and scattering effects account for a large percentage (>70 per cent) of the energy in the observed seismograms. The interpretation of the data is consistent with the results of previous studies of these data, but provides considerably more detail. Overall, the P-wave velocities in the Wichita Uplift are more typical of the middle crust than the upper crust (5.3 to 7.1 km/sec). At the surface, the uplift is either exposed as weathered outcrop (5.0 to 5.3 km/sec) or is overlain with sediments of up to 0.4 km in thickness, ranging in velocity from 2.7 to 3.4 km/sec, generally increasing with depth. The core of the uplift is relatively seismically transparent. A very clear, coherent reflection is observed from the Mountain View fault, which dips at ≈40° to the southwest, to at least 12 km depth. Velocities in the Anadarko Basin are typical of sedimentary basins; there is a general increase from ≈2.7 km/sec at the surface to ≈5.9 km/sec at ≈16 km depth, with discontinuous reflections at depths of ≈8, 10, 12, and 16 km.

scholarly journals Modeling Transient Flows in Heterogeneous Layered Porous Media Using the Space–Time Trefftz Method

2021 ◽  
Vol 11 (8) ◽  
pp. 3421
Author(s):  
Cheng-Yu Ku ◽  
Li-Dan Hong ◽  
Chih-Yu Liu ◽  
Jing-En Xiao ◽  
Wei-Po Huang
Keyword(s):  
Porous Media ◽  
Time Domain ◽  
Numerical Solutions ◽  
Space Time ◽  
Two Dimensional ◽  
Transient Flows ◽  
Layered Porous Media ◽  
Time Marching ◽  
Time Basis ◽  
Marching Scheme

In this study, we developed a novel boundary-type meshless approach for dealing with two-dimensional transient flows in heterogeneous layered porous media. The novelty of the proposed method is that we derived the Trefftz space–time basis function for the two-dimensional diffusion equation in layered porous media in the space–time domain. The continuity conditions at the interface of the subdomains were satisfied in terms of the domain decomposition method. Numerical solutions were approximated based on the superposition principle utilizing the space–time basis functions of the governing equation. Using the space–time collocation scheme, the numerical solutions of the problem were solved with boundary and initial data assigned on the space–time boundaries, which combined spatial and temporal discretizations in the space–time manifold. Accordingly, the transient flows through the heterogeneous layered porous media in the space–time domain could be solved without using a time-marching scheme. Numerical examples and a convergence analysis were carried out to validate the accuracy and the stability of the method. The results illustrate that an excellent agreement with the analytical solution was obtained. Additionally, the proposed method was relatively simple because we only needed to deal with the boundary data, even for the problems in the heterogeneous layered porous media. Finally, when compared with the conventional time-marching scheme, highly accurate solutions were obtained and the error accumulation from the time-marching scheme was avoided.

Transient Two-Dimensional Heat Conduction Problem With Partial Heating Near Corners

2016 ◽  
Author(s):  
Robert L. McMasters ◽  
Filippo de Monte ◽  
James V. Beck ◽  
Donald E. Amos
Keyword(s):  
Heat Conduction ◽  
Numerical Integration ◽  
Analytical Solutions ◽  
Numerical Solutions ◽  
Thermal Protection ◽  
Heat Conduction Problem ◽  
Separation Of Variables ◽  
Two Dimensional ◽  
Rectangular Region ◽  
Long Distance

This paper provides a solution for two-dimensional heating over a rectangular region on a homogeneous plate. It has application to verification of numerical conduction codes as well as direct application for heating and cooling of electronic equipment. Additionally, it can be applied as a direct solution for the inverse heat conduction problem, most notably used in thermal protection systems for re-entry vehicles. The solutions used in this work are generated using Green’s functions. Two approaches are used which provide solutions for either semi-infinite plates or finite plates with isothermal conditions which are located a long distance from the heating. The methods are both efficient numerically and have extreme accuracy, which can be used to provide additional solution verification. The solutions have components that are shown to have physical significance. The extremely precise nature of analytical solutions allows them to be used as prime standards for their respective transient conduction cases. This extreme precision also allows an accurate calculation of heat flux by finite differences between two points of very close proximity which would not be possible with numerical solutions. This is particularly useful near heated surfaces and near corners. Similarly, sensitivity coefficients for parameter estimation problems can be calculated with extreme precision using this same technique. Another contribution of these solutions is the insight that they can bring. Important dimensionless groups are identified and their influence can be more readily seen than with numerical results. For linear problems, basic heating elements on plates, for example, can be solved to aid in understanding more complex cases. Furthermore these basic solutions can be superimposed both in time and space to obtain solutions for numerous other problems. This paper provides an analytical two-dimensional, transient solution for heating over a rectangular region on a homogeneous square plate. Several methods are available for the solution of such problems. One of the most common is the separation of variables (SOV) method. In the standard implementation of the SOV method, convergence can be slow and accuracy lacking. Another method of generating a solution to this problem makes use of time-partitioning which can produce accurate results. However, numerical integration may be required in these cases, which, in some ways, negates the advantages offered by the analytical solutions. The method given herein requires no numerical integration; it also exhibits exponential series convergence and can provide excellent accuracy. The procedure involves the derivation of previously-unknown simpler forms for the summations, in some cases by virtue of the use of algebraic components. Also, a mathematical identity given in this paper can be used for a variety of related problems.

scholarly journals Perturbation of eigenvalues due to gaps in two-dimensional boundaries

2006 ◽  
Vol 463 (2079) ◽  
pp. 759-786 ◽  
Author(s):  
Anthony M.J Davis ◽  
Stefan G Llewellyn Smith
Keyword(s):  
Eigenvalue Problems ◽  
Numerical Solutions ◽  
Neumann Boundary ◽  
Diffusion Problem ◽  
Length Scale ◽  
Two Dimensional ◽  
Constructive Approach ◽  
Term Outcome ◽  
Long Term Outcome

Motivated by problems involving diffusion through small gaps, we revisit two-dimensional eigenvalue problems with localized perturbations to Neumann boundary conditions. We recover the known result that the gravest eigenvalue is O (|ln  ϵ | −1 ), where ϵ is the ratio of the size of the hole to the length-scale of the domain, and provide a simple and constructive approach for summing the inverse logarithm terms and obtaining further corrections. Comparisons with numerical solutions obtained for special geometries, both for the Dirichlet ‘patch problem’ where the perturbation to the boundary consists of a different boundary condition and for the gap problem, confirm that this approach is a simple way of obtaining an accurate value for the gravest eigenvalue and hence the long-term outcome of the underlying diffusion problem.

Similar Solutions for Laminar Film Condensation with Adverse Pressure Gradients

1973 ◽  
Vol 95 (2) ◽  
pp. 268-270 ◽  
Author(s):  
P. M. Beckett
Keyword(s):  
Numerical Solutions ◽  
Saturated Vapor ◽  
Film Condensation ◽  
Pressure Gradients ◽  
Two Dimensional ◽  
Laminar Film ◽  
Approximate Models ◽  
Laminar Film Condensation ◽  
Similar Solutions

Steady two-dimensional laminar film condensation is investigated when the saturated vapor has the Falkner–Skan mainstream. Numerical solutions and approximate models are discussed with reference to other published work.

STACKING OF REFLECTIONS FROM COMPLEX STRUCTURES

Geophysics ◽  
1974 ◽  
Vol 39 (4) ◽  
pp. 427-440 ◽  
Author(s):  
Max K. Miller
Keyword(s):  
Deep Structure ◽  
Control Model ◽  
Layer Model ◽  
Low Velocity ◽  
Seismic Reflection Data ◽  
Near Surface ◽  
Static Correction ◽  
Reflection Data ◽  
Actual Field ◽  
Interval Velocities

Common‐depth‐point seismic reflection data were generated on a computer using simple ray tracing and analyzed with processing techniques currently used on actual field recordings. Constant velocity layers with curved interfaces were used to simulate complex geologic shapes. Two models were chosen to illustrate problems caused by curved geologic interfaces, i.e., interfaces at depths which vary laterally in a nonlinear fashion and produce large spatial variations in the apparent stacking velocity. A three‐layer model with a deep structure and no weathering was used as a control model. For comparison, a low velocity weathering layer also of variable thickness was inserted near the surface of the control model. The low velocity layer was thicker than the ordinary thin weathering layers where state‐of‐the‐art static correction methods work well. Traveltime, moveout, apparent rms velocities, and interval velocities were calculated for both models. The weathering introduces errors into the rms velocities and traveltimes. A method is described to compensate for these errors. A static correction applied to the traveltimes reduced the fluctuation of apparent rms velocities. Values for the thick weathering layer model were “over corrected” so that synclines (anticlines) replaced false anticlines (synclines) for both near‐surface and deep zones. It is concluded that computer modeling is a useful tool for analyzing specific problems of processing CDP seismic data such as errors in velocity estimates produced by large lateral variations in overburden.

Multigrid Numerical Solutions of Non-Isothermal Laminar Recirculating Flows

1999 ◽  
Author(s):  
Marcelo J. S. de Lemos ◽  
Maximilian S. Mesquita
Keyword(s):  
Numerical Solutions ◽  
Rectangular Domain ◽  
Optimal Number ◽  
Computational Effort ◽  
Temperature Fields ◽  
Two Dimensional ◽  
Discretization Scheme ◽  
Multigrid Algorithm ◽  
Finite Volume Discretization ◽  
Velocity And Temperature Fields

Abstract The present work investigates the efficiency of the multigrid numerical method applied to solve two-dimensional laminar velocity and temperature fields inside a rectangular domain. Numerical analysis is based on the finite volume discretization scheme applied to structured orthogonal regular meshes. Performance of the correction storage (CS) multigrid algorithm is compared for different inlet Reynolds number (Rein) and number of grids. Up to four grids were used for both V- and W-cycles. Simultaneous and uncoupled temperature-velocity solution schemes were also applied. Advantages in using more than one grid is discussed. Results further indicate an increase in the computational effort for higher Rein and an optimal number of relaxation sweeps for both V- and W-cycles.

Sign in / Sign up

Export Citation Format

Share Document