Indirect reduced-order modelling: using nonlinear manifolds to conserve kinetic energy.
Ontology highlight
ABSTRACT:
SUBMITTER: Nicolaidou E
PROVIDER: S-EPMC7735295 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM's prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards' equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM's accuracy while speeding up nonlinear solver by 18-73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM's prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM's accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.
| S-EPMC9684583 | biostudies-literature
Project description:There is rich variety in the activity of single neurons recorded during behaviour. Yet, these diverse single neuron responses can be well described by relatively few patterns of neural co-modulation. The study of such low-dimensional structure of neural population activity has provided important insights into how the brain generates behaviour. Virtually all of these studies have used linear dimensionality reduction techniques to estimate these population-wide co-modulation patterns, constraining them to a flat "neural manifold". Here, we hypothesised that since neurons have nonlinear responses and make thousands of distributed and recurrent connections that likely amplify such nonlinearities, neural manifolds should be intrinsically nonlinear. Combining neural population recordings from monkey motor cortex, mouse motor cortex, mouse striatum, and human motor cortex, we show that: 1) neural manifolds are intrinsically nonlinear; 2) the degree of their nonlinearity varies across architecturally distinct brain regions; and 3) manifold nonlinearity becomes more evident during complex tasks that require more varied activity patterns. Simulations using recurrent neural network models confirmed the proposed relationship between circuit connectivity and manifold nonlinearity, including the differences across architecturally distinct regions. Thus, neural manifolds underlying the generation of behaviour are inherently nonlinear, and properly accounting for such nonlinearities will be critical as neuroscientists move towards studying numerous brain regions involved in increasingly complex and naturalistic behaviours.
| S-EPMC10370078 | biostudies-literature
Project description:Modern science confronts us with massive amounts of data: expression profiles of thousands of human genes, multimedia documents, subjective judgments on consumer products or political candidates, trade indices, global climate patterns, etc. These data are often highly structured, but that structure is hidden in a complex set of relationships or high-dimensional abstractions. Here we present a self-organizing algorithm for embedding a set of related observations into a low-dimensional space that preserves the intrinsic dimensionality and metric structure of the data. The embedding is carried out by using an iterative pairwise refinement strategy that attempts to preserve local geometry while maintaining a minimum separation between distant objects. In effect, the method views the proximities between remote objects as lower bounds of their true geodesic distances and uses them as a means to impose global structure. Unlike previous approaches, our method can reveal the underlying geometry of the manifold without intensive nearest-neighbor or shortest-path computations and can reproduce the true geodesic distances of the data points in the low-dimensional embedding without requiring that these distances be estimated from the data sample. More importantly, the method is found to scale linearly with the number of points and can be applied to very large data sets that are intractable by conventional embedding procedures.
| S-EPMC138530 | biostudies-literature
Project description:In this paper, we extend a recently introduced rigid body minimization algorithm, defined on manifolds, to the problem of minimizing the energy of interacting flexible molecules. The goal is to integrate moving the ligand in six dimensional rotational/translational space with internal rotations around rotatable bonds within the two molecules. We show that adding rotational degrees of freedom to the rigid moves of the ligand results in an overall optimization search space that is a manifold to which our manifold optimization approach can be extended. The effectiveness of the method is shown for three different docking problems of increasing complexity. First, we minimize the energy of fragment-size ligands with a single rotatable bond as part of a protein mapping method developed for the identification of binding hot spots. Second, we consider energy minimization for docking a flexible ligand to a rigid protein receptor, an approach frequently used in existing methods. In the third problem, we account for flexibility in both the ligand and the receptor. Results show that minimization using the manifold optimization algorithm is substantially more efficient than minimization using a traditional all-atom optimization algorithm while producing solutions of comparable quality. In addition to the specific problems considered, the method is general enough to be used in a large class of applications such as docking multidomain proteins with flexible hinges. The code is available under open source license (at http://cluspro.bu.edu/Code/Code_Rigtree.tar) and with minimal effort can be incorporated into any molecular modeling package.
| S-EPMC4606885 | biostudies-literature
Project description:Cable bacteria of the family Desulfobulbaceae couple spatially separated sulfur oxidation and oxygen or nitrate reduction by long-distance electron transfer, which can constitute the dominant sulfur oxidation process in shallow sediments. However, it remains unknown how cells in the anoxic part of the centimeter-long filaments conserve energy. We found 16S rRNA gene sequences similar to groundwater cable bacteria in a 1-methylnaphthalene-degrading culture (1MN). Cultivation with elemental sulfur and thiosulfate with ferrihydrite or nitrate as electron acceptors resulted in a first cable bacteria enrichment culture dominated >90% by 16S rRNA sequences belonging to the Desulfobulbaceae. Desulfobulbaceae-specific fluorescence in situ hybridization (FISH) unveiled single cells and filaments of up to several hundred micrometers length to belong to the same species. The Desulfobulbaceae filaments also showed the distinctive cable bacteria morphology with their continuous ridge pattern as revealed by atomic force microscopy. The cable bacteria grew with nitrate as electron acceptor and elemental sulfur and thiosulfate as electron donor, but also by sulfur disproportionation when Fe(Cl)2 or Fe(OH)3 were present as sulfide scavengers. Metabolic reconstruction based on the first nearly complete genome of groundwater cable bacteria revealed the potential for sulfur disproportionation and a chemo-litho-autotrophic metabolism. The presence of different types of hydrogenases in the genome suggests that they can utilize hydrogen as alternative electron donor. Our results imply that cable bacteria not only use sulfide oxidation coupled to oxygen or nitrate reduction by LDET for energy conservation, but sulfur disproportionation might constitute the energy metabolism for cells in large parts of the cable bacterial filaments.
| S-EPMC6976610 | biostudies-literature
Project description:Biochemical systems embed complex networks and hence development and analysis of their detailed models pose a challenge for computation. Coarse-grained biochemical models, called reduced-order models (ROMs), consisting of essential biochemical mechanisms are more useful for computational analysis and for studying important features of a biochemical network. The authors present a novel method to model-reduction by identifying potentially important parameters using multidimensional sensitivity analysis. A ROM is generated for the GTPase-cycle module of m1 muscarinic acetylcholine receptor, Gq, and regulator of G-protein signalling 4 (a GTPase-activating protein or GAP) starting from a detailed model of 48 reactions. The resulting ROM has only 17 reactions. The ROM suggested that complexes of G-protein coupled receptor (GPCR) and GAP--which were proposed in the detailed model as a hypothesis--are required to fit the experimental data. Models previously published in the literature are also simulated and compared with the ROM. Through this comparison, a minimal ROM, that also requires complexes of GPCR and GAP, with just 15 parameters is generated. The proposed reduced-order modelling methodology is scalable to larger networks and provides a general framework for the reduction of models of biochemical systems.
| S-EPMC3417759 | biostudies-other
Project description:IntroductionFractional flow reserve (FFR) improves assessment of the physiological significance of coronary lesions compared with conventional angiography. However, it is an invasive investigation. We tested the performance of a virtual FFR (1D-vFFR) using routine angiographic images and a rapidly performed reduced order computational model.MethodsQuantitative coronary angiography (QCA) was performed in 102 with coronary lesions assessed by invasive FFR. A 1D-vFFR for each lesion was created using reduced order (one-dimensional) computational flow modelling derived from conventional angiographic images and patient specific estimates of coronary flow. The diagnostic accuracy of 1D-vFFR and QCA derived stenosis was compared against the gold standard of invasive FFR using area under the receiver operator characteristic curve (AUC).ResultsQCA revealed the mean coronary stenosis diameter was 44% ± 12% and lesion length 13 ± 7 mm. Following angiography calculation of the 1DvFFR took less than one minute. Coronary stenosis (QCA) had a significant but weak correlation with FFR (r = -0.2, p = 0.04) and poor diagnostic performance to identify lesions with FFR <0.80 (AUC 0.39, p = 0.09), (sensitivity - 58% and specificity - 26% at a QCA stenosis of 50%). In contrast, 1D-vFFR had a better correlation with FFR (r = 0.32, p = 0.01) and significantly better diagnostic performance (AUC 0.67, p = 0.007), (sensitivity - 92% and specificity - 29% at a 1D-vFFR of 0.7).Conclusions1D-vFFR improves the determination of functionally significant coronary lesions compared with conventional angiography without requiring a pressure-wire or hyperaemia induction. It is fast enough to influence immediate clinical decision-making but requires further clinical evaluation.
| S-EPMC7656870 | biostudies-literature
Project description:Patterning is a familiar approach for imparting novel functionalities to free surfaces. We extend the patterning paradigm to interfaces between crystalline solids. Many interfaces have non-uniform internal structures comprised of misfit dislocations, which in turn govern interface properties. We develop and validate a computational strategy for designing interfaces with controlled misfit dislocation patterns by tailoring interface crystallography and composition. Our approach relies on a novel method for predicting the internal structure of interfaces: rather than obtaining it from resource-intensive atomistic simulations, we compute it using an efficient reduced order model based on anisotropic elasticity theory. Moreover, our strategy incorporates interface synthesis as a constraint on the design process. As an illustration, we apply our approach to the design of interfaces with rapid, 1-D point defect diffusion. Patterned interfaces may be integrated into the microstructure of composite materials, markedly improving performance.
| S-EPMC4148653 | biostudies-literature
Project description:Nonlinear susceptibilities are key to ultrafast lightwave driven optoelectronics, allowing petahertz scaling manipulation of the signal. Recent experiments retrieved a 3rd order nonlinear susceptibility by comparing the nonlinear response induced by a strong laser field to a linear response induced by the otherwise identical weak field. The highly nonlinear nature of high harmonic generation (HHG) has the potential to extract even higher order nonlinear susceptibility terms. However, up till now, such characterization has been elusive due to a lack of direct correspondence between high harmonics and nonlinear susceptibilities. Here, we demonstrate a regime where such correspondence can be clearly made, extracting nonlinear susceptibilities (7th, 9th, and 11th) from sapphire of the same order as the measured high harmonics. The extracted high order susceptibilities show angular-resolved periodicities arising from variation in the band structure with crystal orientation. Our results open a door to multi-channel signal processing, controlled by laser polarization.
| S-EPMC6646338 | biostudies-literature
Project description:Many practical systems can be described by dynamic networks, for which modern technique can measure their outputs, and accumulate extremely rich data. Nevertheless, the network structures producing these data are often deeply hidden in the data. The problem of inferring network structures by analyzing the available data, turns to be of great significance. On one hand, networks are often driven by various unknown facts, such as noises. On the other hand, network structures of practical systems are commonly nonlinear, and different nonlinearities can provide rich dynamic features and meaningful functions of realistic networks. Although many works have considered each fact in studying network reconstructions, much less papers have been found to systematically treat both difficulties together. Here we propose to use high-order correlation computations (HOCC) to treat nonlinear dynamics; use two-time correlations to decorrelate effects of network dynamics and noise driving; and use suitable basis and correlator vectors to unifiedly infer all dynamic nonlinearities, topological interaction links and noise statistical structures. All the above theoretical frameworks are constructed in a closed form and numerical simulations fully verify the validity of theoretical predictions.
| S-EPMC5359559 | biostudies-literature