Analysis of one-dimensional Helmholtz equation with PML boundary

AbstractIn this paper, the linear conforming finite element method for the one-dimensional Bérenger's PML boundary is investigated and well-posedness of the given equation is discussed. Furthermore, optimal error estimates and stability in the L2 or H1-norm are derived under the assumption that h, h2ω2 and h2ω3 are sufficiently small, where h is the mesh size and ω denotes a fixed frequency. Numerical examples are presented to validate the theoretical error bounds

Refraction traveltime tomography using damped monochromatic wavefield

For complicated earth models, wave-equation–based refraction-traveltime tomography is more accurate than ray-based tomography but requires more computational effort. Most of the computational effort in traveltime tomography comes from computing traveltimes and their Fr´echet derivatives, which for ray-based methods can be computed directly. However, in most wave-equation traveltime-tomography algorithms, the steepest descent direction of the objective function is computed by the backprojection algorithm, without computing a Fr ´echet derivative directly. We propose a new wave-based refraction-traveltime– tomography procedure that computes Fr´echet derivatives directly and efficiently. Our method involves solving a damped-wave equation using a frequency-domain, finite-element modeling algorithm at a single frequency and invoking the reciprocity theorem. A damping factor, which is commonly used to suppress wraparound effects in frequency-domain modeling, plays the role of suppressing multievent wavefields. By limiting the wavefield to a single first arrival, we are able to extract the first-arrival traveltime from the phase term without applying a time window. Computing the partial derivative of the damped wave-equation solution using the reciprocity theorem enables us to compute the Fr ´echet derivative of amplitude, as well as that of traveltime, with respect to subsurface parameters. Using the Marmousi-2 model, we demonstrate numerically that refraction traveltime tomography with large-offset data can be used to provide the smooth initial velocity model necessary for prestack depth migration.This work was financially supported by the National Laboratory Project of the Ministry of Science and Technology and the Brain Korea 21 project of the Ministry of Education. We are also grateful to Prof. K. J. Marfurt of the University of Houston and Dr. M. Schoenberger for editing our manuscript

Evolution of the Stethoscope: Advances with the Adoption of Machine Learning and Development of Wearable Devices

The stethoscope has long been used for the examination of patients, but the importance of auscultation has declined due to its several limitations and the development of other diagnostic tools. However, auscultation is still recognized as a primary diagnostic device because it is non-invasive and provides valuable information in real-time. To supplement the limitations of existing stethoscopes, digital stethoscopes with machine learning (ML) algorithms have been developed. Thus, now we can record and share respiratory sounds and artificial intelligence (AI)-assisted auscultation using ML algorithms distinguishes the type of sounds. Recently, the demands for remote care and non-face-to-face treatment diseases requiring isolation such as coronavirus disease 2019 (COVID-19) infection increased. To address these problems, wireless and wearable stethoscopes are being developed with the advances in battery technology and integrated sensors. This review provides the history of the stethoscope and classification of respiratory sounds, describes ML algorithms, and introduces new auscultation methods based on AI-assisted analysis and wireless or wearable stethoscopes

Traveltime and amplitude calculations using the damped wave solution

Because of its computational efficiency, prestack Kirchhoff depth migration remains the method of choice for all but the most complicated geological depth structures. Further improvement in computational speed and amplitude estimation will allow us to use such technology more routinely and generate better images. To this end, we developed a new, accurate, and economical algorithm to calculate first-arrival traveltimes and amplitudes for an arbitrarily complex earth model. Our method is based on numerical solutions of the wave equation obtained by using well-established finite-difference or finite-element modeling algorithms in the Laplace domain, where a damping term is naturally incorporated in the wave equation. We show that solving the strongly damped wave equation is equivalent to solving the eikonal and transport equations simultaneously at a fixed reference frequency, which properly accounts for caustics and other problems encountered in ray theory. Using our algorithm, we can easily calculate first-arrival traveltimes for given models. We present numerical examples for 2-D acoustic models having irregular topography and complex geological structure using a finite-element modeling code.This work was financially supported by National Research Laboratory Project of the Korea Ministry of Science and Technology, Brain Korea 21 project of the Korea Ministry of Education, grant No. R05-2000-00003 from the Basic Research Program of the Korea Science&Engineering Foundation, and grant No. PM10300 from Korea Ocean Research & Development Institute

Traveltime and amplitude calculation using a perturbation approach

Accurate amplitudes and correct traveltimes are critical factors that govern the quality of prestack migration images. Because we never know the correct velocity initially, recomputing traveltimes and amplitudes of updated velocity models can dominate the iterative prestack migration procedure. Most tomographic velocity updating techniques require the calculation of the change of traveltime due to local changes in velocity. For such locally updated velocity models, perturbation techniques can be a significantly more economic way of calculating traveltimes and amplitudes than recalculating the entire solutions from scratch. In this paper, we implement an iterative Born perturbation theory applied to the damped wave equation algorithm. Our iterative Born perturbation algorithm yields stable solutions for models having velocity contrasts of 30% about the initial velocity estimate, which is significantly more economic than recalculating the entire solution.This work was financially supported by National Research Laboratory Project of the Korea Ministry of Science and Technology, Brain Korea 21 project of the Korea Ministry of Education, grant No. R05-2000-00003 from the Basic Research Program of the Korea Science&Engineering Foundation, and grant No. PM10300 from Korea Ocean Research & Development Institute

PINK1 deficiency impairs osteoblast differentiation through aberrant mitochondrial homeostasis

Background PTEN-induced kinase 1 (PINK1) is a serine/threonine-protein kinase in mitochondria that is critical for mitochondrial quality control. PINK1 triggers mitophagy, a selective autophagy of mitochondria, and is involved in mitochondrial regeneration. Although increments of mitochondrial biogenesis and activity are known to be crucial during differentiation, data regarding the specific role of PINK1 in osteogenic maturation and bone remodeling are limited. Methods We adopted an ovariectomy model in female wildtype and Pink1−/− mice. Ovariectomized mice were analyzed using micro-CT, H&E staining, Masson’s trichrome staining. RT-PCR, western blot, immunofluorescence, alkaline phosphatase, and alizarin red staining were performed to assess the expression of PINK1 and osteogenic markers in silencing of PINK1 MC3T3-E1 cells. Clinical relevance of PINK1 expression levels was determined via qRT-PCR analysis in normal and osteoporosis patients. Results A significant decrease in bone mass and collagen deposition was observed in the femurs of Pink1−/− mice after ovariectomy. Ex vivo, differentiation of osteoblasts was inhibited upon Pink1 downregulation, accompanied by impaired mitochondrial homeostasis, increased mitochondrial reactive oxygen species production, and defects in mitochondrial calcium handling. Furthermore, PINK1 expression was reduced in bones from patients with osteoporosis, which supports the practical role of PINK1 in human bone disease. Conclusions In this study, we demonstrated that activation of PINK1 is a requisite in osteoblasts during differentiation, which is related to mitochondrial quality control and low reactive oxygen species production. Enhancing PINK1 activity might be a possible treatment target in bone diseases as it can promote a healthy pool of functional mitochondria in osteoblasts.So-Young Lee received National Research Foundation Grant of Korea (NRF2019R1A2C4070492), funded by the Korean government (https://www.nrf.re.kr) for this work. Soonchul Lee received National Research Foundation Grant of Korea (NRF-2019R1C1C1004017), funded by the Korean government (https://www.nrf.re.kr) for this work

Magnetotelluric inversion via reverse time migration algorithm of seismic data

We propose a new algorithm for two-dimensional magnetotelluric (MT) inversion. Our algorithm is an MT inversion based on the steepest descent method, borrowed from the backpropagation technique of seismic inversion or reverse time migration, introduced in the middle 1980s by Lailly and Tarantola. The steepest descent direction can be calculated efficiently by using the symmetry of numerical Green's function derived from a mixed finite element method proposed by Nedelec for Maxwell's equation, without calculating the Jacobian matrix explicitly. We construct three different objective functions by taking the logarithm of the complex apparent resistivity as introduced in the recent waveform inversion algorithm by Shin and Min. These objective functions can be naturally separated into amplitude inversion, phase inversion and simultaneous inversion. We demonstrate our algorithm by showing three inversion results for synthetic data.The work of T. Ha was supported by the Korea Research Foundation Grant (KRF-2006-C00014) and the work of C. Shin was financially supported by the Brain Korea 21 Project of the Korea Ministry of Education

Efficient electric resistivity inversion using adjoint state of mixed finite-element method for Poissons equation

We propose an electric resistivity inversion method that is similar to the reverse time migration technique applied to seismic data. For calculating model responses and inversion, we use the mixed finite-element method with the standard P1 P0 pair for triangular decompositions, which makes it possible to compute both the electric potential and the electric field vector economically. In order to apply the adjoint state of the Poisson equation in the resistivity inverse problem, we introduce an apparent electric field defined as the dot product between the computed electric field vector and a weighting factor and then defining a virtual source to compute the partial derivative of the electric field vector. We exploit the adjoint state (the symmetry of Green s function) of matrix equations derived from solving the Poisson equation by the mixed finiteelement method, for the calculation of the steepest descent direction of our objective function. By computing the steepest descent direction by a dot product of backpropagated residual and virtual source, we can avoid the cumbersome and expensive process of computing the Jacobian matrix directly. We calibrate our algorithm on a synthetic of a buried conductive block and obtain an image that is compatible with the limits of the resistivity method.The work of Ha was supported by the Korea Research Foundation Grant (KRF-2004-C00007) and the works of other people were financially supported by the National Research Laboratory Project of the Korea Ministry of Science and Technology, the Brain Korea 21 Project of the Korea Ministry of Education

Detection of large sequence insertions by a hybrid approach that combine de novo assembly and resequencing of medium-coverage genome sequences

Large sequence insertion (LSI) is one of the structural variations (SVs) that may cause phenotypic differences in plants. To identify the LSIs using medium-coverage sequencing data of four wild soybean (Glycine soja) genotypes, we designed a hybrid approach combining de novo assembly and read mapping. Total reads and reads with both ends unmapped were independently assembled into “ordinary contigs” and “orphan contigs”, respectively, and subjected to pairwise alignment and stringent filtering. This approach predicted 24 LSIs averaging 2,682 bp in size, with no overlap with SVs detected by Pindel, BreakDancer, or ScanIndel, and they were validated by PCR. Compared with the soybean (Glycine max) reference genome, 20 LSIs were located outside genic regions. One of the four LSIs within a genic region, LSI05, is located in the coding DNA sequence region of a protein kinase superfamily gene (Glyma.08G123500). It caused delayed translation initiation and loss of 24 amino acids in the wild soybean genotype CW12. LSI05 was more frequently observed in 29 G. soja accessions than in 34 G. max accessions. Identified LSIs would be genomic resources harboring novel gene contents for studying SVs and improving crops. Moreover, our cost-efficient approach may be applicable to other plant species.The accepted manuscript in pdf format is listed with the files at the bottom of this page. The presentation of the authors' names and (or) special characters in the title of the manuscript may differ slightly between what is listed on this page and what is listed in the pdf file of the accepted manuscript; that in the pdf file of the accepted manuscript is what was submitted by the author